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SUMMARY 


. This report is a user’s guide for a computer code which calculates the 
-nmmar and turbulent hypersonic flows about blunt axisymmetric bodies such a 
»|*»r,e»u y blunted nones, hyperboloids, etc., at zero C small angles of a- 1 
‘.A' C ° de rtril0 ^ en 111 orPAR FORTRAIi language for the CDC-STAR-100 computer, 
ime ependent viscous-shock-layer- type equations are used to describe the ' fl o- 
„^ e 7. r * J Tn ® se equat - ons are solved by an explicit, two-step, time-asymptotic 
f tnite-dliference met nod. For the turbulent flow, a two-layer, eddv-vinccsitv 
model is used. The code provides complete flow-field properties including shoi 
location, surface pressure _ distribution, surface heating rates, and skin-*’' ' 
r j c * ion . cos t f ici ent..; . T.ais report contains descriptions of the input and out r 
■ :ie listing of trie program, and a sample flow-field solution. 


INTRODUCTION 

This report is a user's guide for a computer code which calculates the 
ammar and turbulent hypersonic flows about blunt axisymmetric bodies such a 
-.[■her ic ally blunted cones, hyperboloids, etc., at zero and small angles of aid 
iiMo dependent viscous-shock-layer-type equations are used to describe the flo 
lieid bounded by the body, shock wave, and outflow boundaries as shown in 
ir:ure . . Die equations are put in the conservative form as; 
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iVo independent variable transformations are applied to the governing equate 
! nt ! r,jt brans formation maps the computational domain into a rectangular r 
ir; wnich the body and the shock wave are made boundary mesh lines. The sec, 
-runs Jormation further maps the rectangular region into another plane to all 
;iigher resolution near the surface, which is desirable for resolving the vis 
° W ‘ . time-asymptotic, two-step finite-difference method is used to sol ve 
governing equations. Solutions are obtained only in the planes of symmetry 
the flow field, i.e. the windward and leeward planes. The code given in thi 
report considers the flow of a perfect gas only. For the turbulent flo? a 
two-layer, eddy-yiscosity model Is used. The details of the governing enuati 
erty viscosity model, and the method of solution are given in references 1 

The code is written in STAR FORTRAN language for the CDC-STAR-3 00 comp- Me 
The present method of solution works very efficiently on this computer The “ ' 
results of the code include shock location, surface pressure distribution. 
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surface heating rate and skin friction, and detailed flow properties at all the 
grid points of the flow field. Since it is not necessary to solve the flow in 
tv/o planes at zero angle of attack, a separate program listing is given for the 
zero unrle-of-attack case for solving the flow only in one plane of the flow 

VI eld . 

The program listing for the ang'l e-of-attack case is gi ven in appendix A, 
the program .listing for the zero angle-of- attack case is given in appendix B, 
and a sample input .and output are given in appendixes " and D, respectively. 


CODE STRUCTURE 

The code has a main program which simply calls either BOON or BHYP depending 
upon whether the body is a spherically blunted cone or a hyperboloid. All the 
inputs a^e made in ROOM or nod. These subroutines then calculate the body geome- 
try and starting :■ jlution at all the grid points. Calculations for the time step 
and the nite d ; fferonc in. *. for the predictor and corrector steps are also 
performed in these subroutines. 

Subroutine CHOCK calculates the flow properties immediately behind the shock 
by using the shock relations and updates the shock shape. Subroutine DERV calcu- 
lates the normal derivatives of the velocities and enthalpy. These are required 
for calculating the transport terms in the governing equations. Subroutine VEC2 
calculates the flow vectors U, M, N, and Q of equation (l) at the mesh 
points on the axis of symmetry, of the body and subroutine VEC1 calculates at the 
remaining mesh points in the flow field. Eddy viscosity is calculated in sub- 
routine EDDY which is called after every 25 time steps instead of every time step 
to save on computing time. Subroutine VECU calculates the physical flow quantities 
such as velocities, density, etc., at the mesh points on the axis of symmetry of the 
body and subroutine VEC3 performs these calculations at the remaining mesh points. 
Eor the turbulent flow, subroutine VE05 is called in place of VEC3., In the case 
of zero angle of attack, VECU is eliminated and the physical flow quantities at 
the body axis are also calculated in either VEC3 or VEC5 depending upon whether 
the flow is laminar or turbulent. The boundary conditions on the body surface- 
are applied in subroutine BOUND. 

Finally, subroutine BOON or BHYP prints the surface quantities, shock shape, 
and drag coefficients after every 500 time steps. At this time, the convergence 
of the solution is also checked. The convergence criterion used in these codes 
is that the surface heating rates should not change by more than 1 percent over 
500 time steps. The detailed flow-field properties are printed at all the mesh 
points after every 5000 time steps and after the solution is converged. 


PROGRAM INPUT 

The program Is set up for 15 mesh points in the tangential direction in each- 
plane and 101 mesh points in the normal direction. It uses the properties of air 
under perfect gas assumption. Sutherland’s viscosity law is used to calculate 
the viscosity. The symbolic program name and a brief description of each input 
parameter are given. 
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NBODY Parameter used to choose the type of body. NBODY = 1 for spherical!.' 

blunted cone and NBODY = 2 for hyperboloid. Rest of the input 
parameters are prescribed in subroutine BOON for NBODY = 1 and in 
subroutine BHYP for NBODY = 2. 


cc 

Artificial damping coefficient 

of the order of 0.001. 

FM 

Free- stream Mach number . 


TF 

Free-stream temperature in °K, 


PF 

p 

Free-stream pressure in N/m . 


THC 

Half-tody angle in radians. 


TW 

Wall temperature in °K. 


RN 

Body- nose radius in m. 


ALPHA 

Angle of attack in radians. 


DY 

Mesh spacing in the tangential 

direction nondimensionalized with nose 


radius. For spherically blunted cones, it should be so chosen that 
the juncture point lies approximately in the middle of two mesh points 

BETA Stretching factor greater than 1. The value of BETA has to be chosen 

such that the viscous region near the surface is properly resolved. 

To help the user in choosing the value of BETA, the distances of the 
first mesh point away from the body are given below as a percentage of 
shock-layer thickness for various BETA. With no stretching, this 
distance will be 1 percent for 101 mesh points in the normal direction 

BETA Distance of first mesh point away from the body as 

percentage of the shock-layer thickness. 

2.0 0.8262 

1-5 0.6742 

1.2 0.1+4U 

0.295 

0.1845 

0. 093 

0.05418 

0.0308 

0.01429 

0.0079 


1.1 

1.05 
1.02 
1.01 

1.005 
1.002 
1.001 
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FDT 


THESH 


LTURB 


LMAX 


Time-step parameter equal to 1 or less. It determines what fraction 
of CFL time step be used to march the solution in time. For low 
Reynolds numbers and turbulent flows, FDT may have to be much less th^n 

Initial shock angle in radians between the body surface and shock. 

For large angle bodies, it can be taken as zero but for bodies with 
half angle less than 20°, where the free-stream Mach number 
normal to shock may approach 1 for THESH = 0, a suitable small value 
should be prescribed for THESH. 

Number of time steps up to which the flow remains laminar. Since the 
starting solution is obtained by linear interpolation between the body 
and the "shock values, LTURB should be large enough (4000 or so) so 
that reasonable flow profiles are established before the turbulence 
is turned on. 

Maximum number of time steps. The code will stop after LMAX time steps 
even though the solution may not have converged. 


CRIT Convergence criterion. 

NT Number of mesh points in the tangential direction up to which the 

flow is laminar. The flow becomes fully turbulent for NT + 1. HI 
is an input only for the zero angle-of-attack code. For the angle o ■ 
attack code, the flow is assumed fully turbulent at all mesh points. 


PROGRAM OUTPUT 


The first quantities to be output are the geometry parameters for the body. 
These quantities are: 


Y 

THE 

R 

CUR 


Distance measured along the body surface nondimensionalized with nose 


radius, R . 

n 


Local body angle in radians- 


Local body radius normal to axis nondimensionalized with R n - 
Local body curvature nondimensionalized with R r . 


The complete starting solution is then printed at each grid point in the 
flow field. The details of the quantities printed m it will be given later. 

The code prints gross flow quantities after every 500 time steps. These 
quantities are: 


Y 

S 

US 


Distance along the body nondimensionalized with R Q . 

Shock stand off distance nondimensionalized with R q . 

Shock speed nondimensionalized with free-stream velocity, V o 



CII 


w 1 1 f 1 r>- 


CF 


PW 


3 

Surface heating rate nondimensionalized with p V /2 

,J ^00 oo 

the free -stream density. 

Skin friction coefficient nondimensionalized with p V' 1 /' 1 . 

r oo On' 

2 

Surface pressure nondimensionalized with p V . 

oo oo 


CDP Local pressure drag coefficient. 

CDF Local ski u friction drag coefficient (only for zero angle of attack). 


The code prints detailed flow quantities at each grid point after every ^00 ) 
time steps and after the solution is converged. The same quantities are printed 
in the starting solution also. These quantities are: 


ZN 


U 

V 

WP 

p 

RO 


T 


H 


V i ST 


Distance normal to body nondimensionalized with R . 

Tangential velocity nondimensionalized with V . 

Formal velocity nondimensionalized with V . 

Cross flow velocity derivative with respect to crossflow direct':-".: 
nondimensionalized with V (only for angle of attack). 

oo 

Pressure nondimensionalized with p . 

^oo oo 

Density nondimensionalized with p . 

17 ^00 

Temperature nondimensionalized with free-stream temperature, T . 

o 

Total enthalpy nondimensi onalized with V"' . 

1 ^ oo 

Eddy viscosity nondimensionalized with laminar viscosity. 


GENERAL DISCUSSION 

The code presented in this report calculates the hypersonic flow of a poriV-'t 
gas past axi symmetric bodies at zero and small angles of attack. No-slip and 
no-mass flux boundary conditions are used at the surface. The flow-field solution 
for the sample Input conditions, given in appendix C, required 7500 time stops lor 
the zero angle of attack and 10 000 time steps for 5° angle of 'attack. !i took 
about 5 minutes for the zero angle-of-attack case and about 15 minutes for . do V 
angle-of-attack case. 

The code works well over a wide range of input conditions and body angles. 
With suitable modifications, it has been used for the flow with low Reynold;, 
numbers (ref. 3) as well as for the flow with massive surface blowing, (rot. 1). 
Recently, the code has been made operational for calculating the aerothornu..] 
environment of the Jovian entry probe (ref. d) in which the solutions are obtained 
for radiating and reacting flow under chemical equilibrium and massive ablation 
injection. The code also works well for bodies with ablated nose shapes. 
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APPENDIX A 


LISTING FOR ANGLE- OF- ATTACK CODE 


PFOGRA m RROPjY ( INPUT, OUTPUT, TAPF5= INPUT, TAPEfceOtlTPlJT) 

TIME-DEPfmOFNT PIWJTE-filFPEOENCE W E THOO FOR CALCULATING THE LAMINAR 
AND TURBULENT HYPERSONIC FL.OS ABOUT BL'JNT AXISYmkeTRIC BOOIES AT 

small angle oe attack, local TIME STEP TS USFO TO march the 
solution in TIME, 

REFERENCE A T A a PAPER no, 77.17? (KIJMaR ANO GRAVES), 

NHOO Ys 1 FOR SPHERE. CONE, CALL OCON, 

NROOYc? for hyperroloto, call hhyp. 

MFODYsI 

IFIMIOOy' EO.1 ) C A 1 l rc on 
TFfNBOOY.EO.?) CALL RHYP 
STOP 
F no 
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COROUTINE SCON 

DIMFMSIOn CPI (IS) ,CP2(l«>),C0Pl MS),Cr>P2(l5),COP(tS) 

DIMENSION v (?9) ,THE f?0),P(?q) , CM (?Q). TEMPI ( 29 ) , TEMP2 ( 1 0 1) , Z d 0 1) , 
172(101), TZ1 (!On,GM(«j) 

DIMENSION TaU(? 9,5) ,TAM(29 > S),TANi(2g,5),TAN2(?9,5),TAQ(29,5) # 

1TAU1 (29,5) ,CF(29) 

DIMENSION TEA (29 I ,TF9(?9) ,TF10(29) 

COMMON /F /THM (1 0 1 ,29),THE2d (U , 29), VAIRft 01,29), 07(101,29) 

COMMON /F1/N1, Nil, M2, N1h,M J5, Ml, MU, M50.N19 

COMMON/F^ /G M 1 ,0M2,OM3,r,M<i,r,M5,r.Mf 5 ,GM7,RMfl,CM9,e2,RE,TW,GAMA,SIG,MM 

CPMM0M/F2/N¥M,NXMt,NXM3,MXM<|,NXM8,NXM5,NXM61 

COMMON/F/J/XX, ZZ,FM,META,SIRT 

C0MM0N/F5/Sf?91 ,SS(291 ,0(291 ,CUR(29) ,U3(2 a )*DS(?9),VS(29) ,VN(29) 
C0MM0N/F6 / Int ( ?9 ) , T I md ( ? gi , temp ( 29 ) , TE (29) , TE l { 29 ) , TE 2 (29) # TE 3 (29) 
1,TE9(29),TF5(29),TF(,(?9),TE7(29) 

C0MM0N/F7/U( 1 01 , 29) , V (101 ,29) , T ( loi ,29) ,P(101 ,29) ,Rfl( 101 ,29) , 

1 WP( t 01 ,29) , VTS(IO) ,291 ,H( 101 ,293 ,SM(t 01 ,29),PW{29) , VIST (10 t ,29) 
COMMON/FP/BT ( 1 01 ,20) ,EM(5A ( 1 01 ,29) ,ZN( 1 01 ,29) 

COMMON/Fo/SI ( 101 ,?0) ,SS1 ( 101 ,29) ,DS1 f 1 01 ,29) , CUR 1 ( 1 0 1 , 29 ) , R1 (1 0 1 , 
129) 

fOMMON/Pl O/THEC (lot, 29), THFSd 01,29), TV1 f 101,29), TV2H01 ,29) 

COMMON /FI t/P!j( 101 ,29),riSH( 101 ,29),DWP( 101, 29), *6(10 1,29), 

1*7(101, 20), *«(101, 20) 

COMMON /F 1 2/ T 7 2 ( 1 0 i , ?9 1 , 2 1 ( i o ! , 29 ) , T 73 ( 1 0 1 , 29) , TZO ( 1 0 1 , 29 ) 

COMMON/Fl J/4II1 (101 ,29,5), AU2(1 01 ,?9,5) 

COMMON /F 10/ A 1 1 ( 10T ,29,5) , * m d o i ,29,5) ,**1(101 ,29,5) ,*0(101 ,29,5) 
COMMON/Fl ?/4l'0( 101 ,«) , *M0( 101 ,2,9) , ANfldOl ,R),*^0i 1 01 ,9) 
rOMM0N/Flfc/HTNF,PEl,RE2 

V I S i ! , 1 | 2920 ) s 0 , 

vn;T(l,l}2929)sO. 

r. ( 1 , S 1 1 1 u s ) a n . 

**■ ( 1 , t , 1 1 l«(*a5)so, 

A \' ( 1 , 1 , 1 | 1 0 > '1 5 ) 3 n , 

ion# i , i » i ubusioo. 

*01(1*1,1)1 0*05)30, 

*1 5 { 1 , i , 1 ) 1 1=0. 

CH( 1 )2Q)=0. 

C CONSTANTS Ain r RFEST9f:*M CONDITIONS 

P T s u , *■ a T a N ; : , ) 

G A M A 3 1 U 
R G A S = 2 H 7 

5 T G = , 7 2 

CFsG*’1**RP*S/(r.AMA-i . 1 

SIGTs.9 

CCo.OOl 

FMti 0,3 

TFsAJ6,2* 
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PFet00.77 
ROFaPF/TF/RGAS 
ALPHA=5.*PI/180. 

THC=US.*P1/1R0. 

TWaJSO.N 
ThaTK/TF 
RMs,0?1 75 
OYa.ieS 
RFT As l , 1 
FOTal, 

THfSH=?,*PI/lPO, 

'/FoFM*(Gama*RGAS*TF)**,« 5 
VISFsl ..'i-.T/l :',**fc*TF**l ,5/fTF+l 10,) 
nFoVF*RnF*«N/viSF 

|,M fi s is THF N'JMRFR OF TJMp STEPS AFTER whJCH THE PROGRAM 

will. STOP FVEN THOUGH the SOLUTION may NOT BE CCMVERGFD, 

L *' A X = 2 0 < 5 ft T 
LTItRB = ?Snno 

FOR LA«tNAR FLO* r ALL VEC3 AND FOR ’LRmn.ENT PLOW, CALL VFC5 
LTIJPH IS THF N‘ '*m}£R 0 F ITERATIONS UP TO WHICH THE FLOW REMAINS 
I. ahtuar. 

CONVERGENCE CRITERION, 
r « I T r , fit 
00 5 Is ' ,5 
5 G M (I)3f-1 . 1**T 

G M * p1 . /GAM A /fm**? 

GM?=2,*GAMA/f \ ,+GAHA ) 

GM^3CP*TF/VF**a 
G u fisSlN ( ALPHA 5 
G M SsUft,/TF 
r,M(,ai ,*rms 

‘ G m 73(GAMA-1 . W(GAMA + 1 , ) 

G M Re ( 1 ,+gAMAW(?,*GAMa*Fm**P) 

RMq a (GAMA««l,l/GAHA 
GMlfi = Tw*r.M3/GAMA 
HINFc.QRS* ( ,5+GMi) 

RF 1 aSORT ( RF ) 

RF?30 , n i pfl*RE 

C MFSH SI7FS ANO VFCTOP LENGTHS 

N 1 a 1 5 
N 1 1 a2*N 1 -1 
N 1 ?a N 1 + 1 
NIUsNi-l 
N15SN] 1-1 
N503N1-2 
NSlaUl 1-2 
NIJaNUj 
N173N1+U 
NSpaNl-3 
NS3SN1 J-3 

NlficNU2 . _ . . . . "L ' 
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Mia 101 
MUaHi.i 
MSOaM *? 

MSl»Ml-5 

M*52»Mi.a 

NXMaNt 1 ★ m 1 

NXHlaMVM.Mjj 

NYM^oNXM.i 

NXH3SNXM.? 

NXMaaKiX w -Ml 

NXM5afgta*Ml 

NXMfcBN»51 *MJ 
NYM61 8Ki5t **1-1 
WXH7»NXm i ,3*m^ + 1 
WXMRtKlXM.Hi.^ 

NXM<3eNXM«£i*M 1 ♦! 

MXM10c5*n/XH 

NX** 1 a q*KjX*fl 
WXMM?aNXM.?^M j »u 

NXM*3aMX*»4**t»i? 

c step sizes 

nZa»l , /Mil 
XXa?,*OY 
ZZ=?.*o I 

v c 1 )aQBVTKiTL( 0 , ,r>Y| Y(1 |Ni ) ) 

ZM |Ml)aQAV!MTL(t,»DZlZ(l|M|)} 

AlePI/?,.THC 

YMcnY*Ml« 

RR»( (Ym -4 j )«SIN(THr.)+COS(THC) )*RN 
C WJ! IS the *ESH PRINT UPSTREAM OF JUNTllBE POINT. 

NjsAl/OY 
NJ1«NJ+J 
NJ?*NJ+N1 

0 

NJJlONJt.l 

NJi2aNjj.2 

N.»2laNJ2»l 

NJ2?rNJ?-2 

c initial Shock Shape ano shock slopf. 
sn»Nin«,i7 
on 1 5 nb3,ni 

15 S(N)«S(N-1 > ♦r>Y*SIN(THFSM) 

ssn )«cs(?).s(Ni 2 n/xx 

00 16 N*?,NJ 

16 ssfN)*(8fNTi >-s(N«m/xx 
SS(Nl2)afS(Nt 6)-S(!))/XX 
00 17 NaNJl.M 

17 SS(N)a.(a,*S(N-l )-3,*S(N)-8(N«2) J/XX 

00 18 1 6, nJ3 

18 SS(N)«(S(NTM-S(N-1))/XX 
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DO 1*> NaNJ2.Nl i 

10 sS(N)a.(«,*S(N«n-3.*S(N)-$(N-2))/XX 

C BODY GEOMETRY. 

00 10 Nal.N) 

IFfY(N) .(JT, A1 )00 TO 11 
THE(N)sPT/2.-Y(N) 

CUR(N)a1 . 

R(N)cCOS(THE(N)) 

go to io 

11 THF(N)*THe 
CUR(N)aO’. 

R(N)sC03(THC)*(Y(N)-4t )*SIN(TMC) 

10 CONTINUE 

Y(N12|N1 <0*Y(2»N1 BJ 
R(Nl2|NlB)o«f2|MU) 

THE(N12rNltt)«THE (2|NU) 

CUR(N12iN1B1«CUR(2jN14> 

00 20 Nr 1 f N i 
VS(N)sC0S(Thf(N1-4LPM4) 

20 VN(N)s9IN(ALPH4-THEfN)) 

00 25 Nbn 1 2 * N 1 1 
VS(N)«CPS{4LPW4+TmE(N)) 

25 vn(n)3-$tn ( alpha+Thf (N) ) 

00 35 Not, Nil 
81 ( 1 » N 1 Ml ) *S (N) 

cur i n,N,MnscuR(N> 

R 1 (l,N|MUsPfN) 

35 THE1 ( 1 ,N|M1)bThe (n> 

THECt 1 , 1 |NXM)«VC0S(THE1 (1 , 1 fNXM) iTHECtl , 1 |NXM)) 

THES(1, 1 |MXM)nVS!NfTHEl (l,l |NXM) jTHfSil , l |NXM)) 
Blr(HETA*t.)/(BETA-l.) 

B2b2#*BETA/AL0G(B1 ) 

XL 1 » -B2/(nFT A**?-l , ) 

XLs-2.*B?/( BETA ** 2-1 .J/SIG/RE 
TZ(iiMi)a«i**(zn »Min , 

TZ1(tl M l)aHE TA*(TZ(ll M l)"l,)/(TZ(ltMl)*1.) 

00 30 Nai ,Ni i 
30 Z! (l»N,Mt)8TZ1 (1 tMD 

ZNC1 , 1 INXMISSI (1 , 1 JN*M)*(1 ,-ZI (1, 1 |N*M) 1 
E«0A(1,1 »NXM)*1 , + ZN(l ,1 |NXM)*CUR 1(1,1 JNXM) 

BTC1 » 1 |NXH)»P1 (1 , I |NXN)*ZN( 1 , 1 »NXM) *THEC ( 1 , 1 |NX*) 

TZ2( 1 » 1 »NXB)oRETA**2-Zl ( 1 , 1 |NXM)**2 
TZ3(1#1 |NXM)e2.*Z1 (1 ,1 »NXM)/TZ2(1,1|NXM) 

TZfltl. 1 |NXM)oB2/TZ?(1 ,1 »MXM) 

TNT d »N1 1 JsOBVlNTL d,BljINT(l|Nin) 

C CALCULATION OP The STARTING SOLUTION'. 

00 00 Na 1 , N 1 

00 TpMP(N)a3IN(THE (N)-ALPHA)**2*{ ,R+, 1*(M-1)/N1 ) 

on 05 N«N12,Nll 

05 TEMP(N)aSIN (THE (N1* ALPHA) **2* (.«»♦, 1*(N-Nn/Nl) 

Pd#l|NXM)»QBV8CATR(TEHPd»Nin,lNT(HNin»Pd,l»NXH) J 
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TEMPdtNinsO. ----- - - — - 

UCl* l|MXM)snRWSCATR(TEMPM ,N11 ),INT(1 |N1 i ) l ,NVMn 

v(tii>NXM)aQgvscATecTE w 0fitnin,iNTfi|NJin>vn.ifiMXMn 

WPM # 1 |NXM)sO«VSCATR(TEMP(l ,N V. > , IN* ' S « Kl 1 1 ) HP(1 , 1 |MX “) ) 
TE M P(!iM!)aT« 

Tn«ljNXM)aQBVSCAT»fTE«P(lr‘Jlt',lNTntMM)r' r n,!jNXM)) 
us!i|Nt! )sss(i»»»n VM,+S(iiwii)*ruP(i»Nin) 

T 1 * iU I N 1 \ ) si *TE5( 1 |M1J )*TF5< 1 (Ml) 

TF6 ( ’ » M nsM-jp - (TF*( ) f M 1 ) j TFM1 }M 1) ) 

*ft T i : i i vi n-vs ' '.^ i n*TF5n j^i i ) )/TF6c i ini i ) 

TF wp ( 1 , m i 1 ) s ( T E 7 ( ( I M t ) ) *TF7 f 1 | 11 ) /f,Mfl-r,-7 ) ♦r.v \ 

P(Mli1^XHl) 3 QHVSC4TF'fTtHPf|.NH),IvT'ljVn)rf , l K ‘l,l^XMl)) 

T EHPf 1 f VI •vm/fjM; 

TF ( 1 J N 1 1 )S (TPMPfl |H| ( ' «-GM ) / ( l". ♦TEMPfl t K < J ) 

TEi c i » m « i )a n ./tf n (»." i ) )hf’( • jM' 1 * /■'E'j'i «ni n 
TF?f 1 IM 1 =V.SM (M n-TFI n |N1 1 )*TES( » ,M n 

■’’F?. ' 1 I M 1 f i ;M n+TFI ( t » M 1 ) 

Vr a ( l ini narpwpn j*-i 1 1 /rr d |N* •. ) 

H'Ml # 1 | MX^n*I3flVSC4TPf TF2( 1 |N1 n # INTO |N1 n |U(Mt , 1 |NXM1 ) ) 

V(‘l M |NXMnsOHVSC4TRfTF3( 1 ,N1 t ),INT ( 1 ,MJ 1 ) I V(M1 , , |NXM1 ) 

Tt«l , 1 »NX^neO0VSr'ATR(Tf«n f N 1 1 ), INTO |M 1 ) |T (Ml , 1 ,NXM ) ) 

DO 61 N=2,N1 

DDF l a -f S f N )»8 fN + wj 0)1/2, 

M WP(Ml #N)a“D“o»Tfi (N)*DDEL/DTf M l , N ) 

DO 65 NaN12,NJ 1 

DDFLb (S (n-M fl )-S(M ) /?. 

65 WP(mi # M) a r,Mu.TFi (v) *DDEL/BT C m 1 , M ) 

WPJMl , i )a-U(Ml, l ) 

DO 50 Nal ,N) 1 
DXIlsfUfMi ,N).ilO,N) ) /Ml 1 
DXVa( V(M| , M-VO ,N) ) /MJ1 
f;:Ve(P(Hl ,NI.F(|,N))/«n 
0 XTb(T (M , N).T O ,N) )/Ml | 

DXwPa(wPfMj,w).wPO,N))/Ml) 

U(1 ,N|M) 1 JsnflVlNTl, (HO ,N) ,DXLJ f l j ( 1 ,N;M| \ ) ) 

V(1 ,Nf MJ I )S0AV!NTl(VO , N) , oxv 1 V( 1,N|M1 1 ) ) 

PM »N|MJ i)sO«VINTL(Pfl #M) #DXP Jp£l ,N |Mt 1) ) 
TO,N|Mli)aOflVINTL(TO t *l),DXT|TO,NjMH)) 

50 WP( l»N|MO )s08VINTLfWP(l f N) ,DXwP)hPO , N,M 1 1 ) ) 

POO , l |NXM)sP(l # 1 jNX«)/TO , 1 |NXM)/GM| 

SH( 1 » 1 |NXM)sT(l,l )NXM)*GM5 

HO , 1 |NXM)aSM(l , 1 jNXM)* fUO , 1 ( NX M ) « *? + V ( 1 , 1 | N X M ) * *2 ) / 2 . 
TV2O»ll w XM) B vS0RT(TO,l|NXM),TV2O,l|NXM)) 

VISO * t |NXM)bT( ! , t |NXM)#TVgO, 1 |NXM)*r,M6/(T( 1 , 1 |NXM)+GH5) 

WRITE (6,39<>) ALPHA, re, FH 

399 F0RMAT(/,!X,*ALPHA',Ffl.5,5X,«ReVN NO , ' , F 1 5. 3 , 5X , *M A CH N0.',F*,J,/) 
WRITE C6, 035) (Y(N),THFCN),R(N),CUR(N) ,Hel ,N| 1 ) 

035 FORMAT (2X, * Ye' , F 1 0 .5. SX , 'THE*' , F 1 0 .5, 5X, 'fis' , F 1 0 .5, 5X , 'CURr ' , F 1 0 , 5 

DO 60 Nat, No, 2 
NRlTE(6#92t )N,Y(N) 


3 •; 


r 



WRITE(6,430) 

KRJTE(6,a«ni ( ZICN,N),U(M,N),V(M,N),HP(M,N),P(M,N),RO(M,N 

60 CONTINUE 

430 FORMAT!/, 6X,'Zl',<*X,'n', < ’X,'V',flX#'*P'. < >X # 'P',eX, , R0', 

15X,M*,9X» # M») 

400 FORMAT! J X» BF 1 0,5) 
pMn-l,/u'./Ml 1 

TFMP2( 1 |M! nsOAVlNTLCi .♦OM.OM|TFMPP( 1 im n > 

AtJ 1 C 1 » t ,1 |NXM10 )s1 . 

T r! '* , .|N51) s l^ 

TE9(ljN5t)=o. 

tf * on |N«5t )sr.«« n 
1^1 

c marching in time: starts here. 

1 CONTINUE 

MMal 

C LOCAL TT M E STEP CALCULATION, 

VMR(1 , 1 |NXm ? 1sG 4MA«PM , 1 jNXMJM/RfHI , 1 JN.YM2) 

VAIRM , 1 jNXM?)sVSO»T ( y/AIR( 1 , 1 jNX'ip) , VATRf 1 , ' |NXMp ( ) 

OTM # 1 |MXM2>bFOT* f ?N(P, 1 ,NYMP)«ZN( < , t fNXvpt ) /( VAR5(V( 1 , 1 |NXM2) J 
1THF2( 1 , 1 | NXM?) )*VAIR( 1 , 1 JNXM2 ) ) 

PO 70 Nal ,NH 

70 OTC!,N»mi n*nT<l,N,Mll)*TFMP2(t ,Mn ) 

TIMP( t »N1 t IsO«VGATHM(OT(Mlt , 1 »NXM1 > , TNT(1 t NI n fTlMPH ,N1 n 1 
DTfMl,1|NXMn*QflVSCATP(TIMP(l|Nl1),INT(1|Nn)|0T(M1 f l|NXMl)) 

C PRFOICTPP STEP, 

SS(1)b(S(?)-S(N12))/XX 
on 74 NS?,NJ 

74 SS(N)3fS(N*1)-8(N«n)2XX 
ssfNi2)ofS(Nih)-sn n/xx 
DO 75 NsMjt.N1 

75 SS(N)a-(u,*S(N-n-3,*S(N)-S(N»2))/XX 
DO 76 NSN16.NJ3 

76 SS(N) = C3CN*n-SCN-1))/XX 

DO 77 NaNJ2.N1 1 

77 SSfN)s-fa.*S(N-n-3.*S(N)-SfN-?))XXX 
CALL SHOCK 

call DFRV 
CALL VEC1 
DO 110 1*1,5 

AUl (2.2.1 |NXM6)3AUf2.?,I tNXM6).DT (?,?}NXM*0*( (AM(2.3»T jNXMfj).AM(2, 

12, I»NXM6) )/nVA(AN(3,2, I»NXM6)-AN(2,2, I JNXM6) ) /OZ + Aij C2.2. I iNXMft) ) 
AIM (2,N1 ,1 I M 50)*AUt2,Nl , l »M50)«0Tf2,Mt ,M50)*( (AM(?,NI , I »MS0)-AM(2, 

lN14*I»M50))/OY*f ANf3,Nl,I»M50)-AN(2,Nl, I j M50 )) /0Z + 40 ( 2, N 1 , I J M 5 0 1 ) 
AUl (2,N1 1,1 •m50)3AU(?,N1 1 ,! »M50)-0T(?,N1 IJMSOI a( (AM(2,N1 1 ,1 »M50)* 

1 AM( 2 ,N 15 ,I» M 50 ))/DY+(AN(S f N| 1 , T jM 5 <n.AM( 2 ,Nl 1 , I JM 50 ) )/OZ*AQ(?,Nll, 

21 ; H50 ) ) 

AUl (2.NJ1 , I }M50)aAU(2,NJl , T | M50 ) +0T ( ? , M Jt | *50 ) * ( ( 4 . * AM ( 2, N J 1 1 , I » 
1MS0)-3.*AM(2,NJ1 , I|M«jO)-AM(?,NJ12,I JM50n/XX-<A^l(3,^»Jt,I|M50)- 
^ANr?,NJl l IlM5<n)/nz..AfH^,HJJ,I|M501) . .. .... 


lit 



AIM (2# N jp , I jM50)sAllf?,NJ?, T | M 50 ) ♦OT (?, N J? j M50 ) * ( ( a .i AH ( 2, K jp j , j , 
1N50)»3 ,*a*(2,nJ2,t j mso ) - am (2, M.J22, Ij; 4 50))/XX«.(AN(3,NJ2,I|M50)«AN( 
22,NJ2, J jm< 50) l /OZ-AO(2,NJ2,I t M50) ) 

TAU( 1 » I |NH 5s08VRATHR( AIMMJ , J , J t MX“n , INT ft JN1 1 ) »TAU( 1 , I ji^M 1 ) ) 

TAMf 1 , I |N1 n*QBVr,ATHB<AN(FM , 1 , T »NVM1 ) , JNT ( 1 »N1 l ) |TAM( 1 , I - Ml I ) 

TAN1 ( 1 , ! |N1 MsQftVRATHPfAMfMl , t , I }NXM 1 ) . INT ? 1 S ' i 1 ' f TAN' r 1 , J j^n ) ) 
TAM?(1 , I »M 1 taOSVGATHPfAN'M i , \ , i ,n X mu , TNTfl |M n |TAN2( ; , I ,M 1) ) 
TAG ft , I |N1 1 ) rO.qvr.ATH'! ( AM 11 1 , 1 , t [Nkm^ , ’kiTf 1 | AM 1 ’ |TA7( t , I f M t) ) 
TA(M'?,!'N' = ;eTAUf?,!|MtSt"Tr'Pf2 S Nt«i)*tfTAH(3,I,MS)..TAH{p l I|N15t 
I ' / "! V + ( r A M (2,1 lMS)» T AN2(2.J I N 1 5 > )/07*TA0f2# T j M 5 1 1 
TAtJ< (fJt « I')=' r All£All, I ) — TIMPfM! >*C (TAM(Nt # I3»TAMC^ltl».T)')/L>Y+(TAA 1 t (Nl # 

tn»TA^’(Mt,nvDZ+ T A:i:'M.i)' 

TAlM(Nll,I)sTAU(N1t # It«TIMpfNll}*((TAf 4 (Nll»I)-TAM(M5 # J))/nV*{TANt 

i ( M 1 # ’ ) • T A N ? ( N 1 1 » I ) ) /l>Z + T A 0 ( M J 1 f I ) } 

TAtJt CNJ1 , T taTA'JfN M ,n*T!^P(MJn*f (<l'.*TAH(NJn , n» 3 ,y.'AHfNJi , I J. 

ITAN'NJip.n t/XX-fTA-: (NJt , t) -T A N2 ( N J i , n 1 /OZ-T A Q ( N J I , n ) 

TAJt (NJ2, nsTA»'(NJ2, T 1 ♦TINP(M,I2)* ( («1.*TAM(NJ21 , It-3.*TAHfNJ2,I)« 

? TAM (NJ22, I) ) /XX.fTAM fMJ2* n-TAN2(NJ?. I) ) /OZ-TAfJ (Nj2, j ) ) 

15 0 AIM ( M 1 T ;NyMM)-08VSCATR(TAUH?,I jNlSt.IN’f nM S) |AU 1 (“1 , 2 , I I 

CAU VFC? 

DO ISO Jat.a 

AIM (2f t , 1 1 u Sft)BA'm(?, I|KS0)»nT(2» 1 »h«> 0 ) * ( ( A^O f 2, ?, 1 1 MSO ) - amo ( 2, t , I 
1 »M50) )/0 v + camo( 3»I i M sro.ANn(2,i;M«50)t/nz4AQn(?,i|"50)) 

1 SO AIM (*1 , 1 , I )*AtlO(Ml ,1 ).nr f M 1 , 1 5 *( ( AMOf Ml , p, I J-AMtlfMj , 1 , I ) )/0Y+ ( ANOf 
1 m t » I )-ANn (MJ « , I ) )/nz*4Q0(Mi ,1 ) ) 

AIMn,3,1|NXv7)=GAVSCATPfTKPniNSn,TNTn,M5l)|Alll(l,3,l»NX v 7n 
AU1 (l*3,2|NXM7)S(3«VSCAT«f TFOf 1 ,M5t ), INTfl |NSn I AM n ,3,2,Nxy7) ) 
AU1(l > 3#3tNXM7)5GAN/SCATP(Tr<>(|{N5lt,TNTf1|NSlt|A(.llft#3#3|NX w 7)) 

AIM n .3,U|NXM7)*0flvsCATR(TFP( 1 |A.Sn , (NTf 1 |N5n I AIM ( 1 ,3,<4|NXN7n 
AU! (l,3.S|NXM7)s0ftVSCATR(TF10n |NSn,INT(| |N51) |At.M Cl #3,5|NXM7) ) 
IMUGT,I.TtJRfOGO TO 15U 
CALL VEC3 
fin TO 1^5 
ISO CALL VFCS 
1SS CONTINUE 
call VFCO 
CALL BODNO 

C CORRECTOR STEP. 

00 160 I«1,S 

160 AU2f2,2, 1 f Nxvfl) 0 AU(2,2, I iNXMBJ+AlM (2.2,! |NXMfl) 

00 170 Tcl ,0 

170 AI)2(2,l,I|M|nsAU0(2,I»Mtl)+AUt C2 # l,T|Mii) 
mm»2 

SSfl )e(8(2)-8CM2))/XX 
00 170 Na2.NJ 

170 ss(N)a(S(N*n.s(N*in/xx 
8S(N12)«f3(N|fc)-S( 1))/XX 
OP 175 NbNJi.N1 

175 SS(N)a-(fl,*S(N-l)-3.*S(N}-S(N-2n/XX 

00 176 NBN1B.NJ3 



176 

177 


160 


190 


200 


210 


ssfN)«(3(N+n-s(M-n )/xx 

DO 177 NeNJ2.M11 

$$(N)e»f<l,*S{M»l)*»3.*S(N)»S(N*»2))/XX 
CALL SHOCt' 

CALL DFRV 
CALI VEC 1 

AU1 c2°?!nNXMfl)*0,5*(AU?(2,2,I|MXMfl).im2,2|NXM6)*((AM(?,2,llNXH61 

j. A Mf 2 , 1 ,1 ,NXM«n/0Y + (ANC2,2,I jtgXMH).AN( 1 , ? , J |NXMfl))/DZ*AQ( 2 * 2 #I | 

^uuzInji .1 }Min*n.5*f AU2(?,NJ1 ,i|Min*nrf2.Nji|HHj*((a,*AM(2 # nji 

PA^n.^Ji . i» m u ) ) /nz-An(?,Njt , i jit n ) ) 

At'1 (2.WJ2, 1 j«1 nso.S*(Alt?(? # NJ2, 1 }M 1 n + DT(?,NJ?|MU)*f (<t.*AM(2,NJ2 

?an( 1 ,nj 2 , i »mii n /r>z-*o(P»MJ2, 1 1 * 1 * 1 1 1 1 
*mM . I , I iMHefiM (I 1 *A M ( 1 , t , I l M n 

Auif?.M?#i;fins,s*fAU?r?,Ni2,n M in-0Tf?,M2» M n)*((A M (?. , ' ! i2»it 

iMiii-AM(? f i f i, M iin/ov*(ANca f Ni2 f i|M,n.AN(i f Ni2 f i,Hii))/nz + A0(2, 

?Mi?,iifiun 

1 1 ( t f N 1 2 t M 1 ■' S r !J n , N 1 2 J n 

CAU VPC2 

fourth oder damping, cc is thf. damping ooefficif^t, 

Af2(3%! I fNXMMgjBAUtS, 2» 1 »NXMM?) + AI|( 1 , 2, t »nxmk?)-(i.M AU((I,2,I |NXMM 
121*AtJ{2,2. 7 jNXNm?}) 3. ?. T 

TFU i I NS n SAP VP at HP f AH (1 ,?. T |NXM7) , JM T (1 |M51) I TFl f ll^l) 1 
TF2C1 1 hSl) «0 AVG » Tmh ( AU(2» 2* 1 1NXM7),! NT r 1 jNSI) jTF2(1 l N M)) 

TE3( 1 J N5 1 )af)flVGATHfi(Af (3,2,1 |MXM7) , INTf 1 :.NSn |TE3( 1 I^ID 
TEU (1 »NS1 )Sf>BVr,ATM»( AU(*1#2, I |MXH 71 , I NT ( 1 I MSI ) t TEa ( 1 » MSI ) ) 

TFSniNm) = TFa(ljN e >n-TF1(1|NSi)-3.*fTF3ni.NSU» i TE2(!)N5n} 

Att?(2,2.1|MXM7) S n6VSCATH(TFS(lr^Sn,riTd |MSn|Ai»2(2,2, I»NXM7)) 

TF1 (1 jN 5 Os 0 «vr,ATHH(Al.l(MSl .2,1 » , ' ; XM7) ,TNT (t | *5 1) I TF 1 ( 1 | MS 1 )) 

TF?(1 tM5l)rOPVr,ATHfi( AH (f SO, 2, It k iXH7) ,1MTM |MSn?TE2(l |M5U ) 

TE3CllN51) = nfiVr,ATH(?(An(M 11 , ? ,il f| i< M 7),IMT(l,^Sn»TF3(l|MSl )) 
TF0niM51)sn8VGATHH(AUC^l,2,I»MXM7),INT(1»MSniTEa(lfNSm 
TF«:( 1 |fSn*TF1 ( 1 jMSn-TEHf 1 »M5n-3.*fTE2( 1 » N 51 )-TE3(l |MS1) i 
Iu2(Mn,?,I|NXM7)sOBVSCATR(TE5(ljMSn,IMTn |N51)|AU2(M11,2,I(NXM7) 
1) 

A02C3»t!l!'*S?)«AU0CS # I»M52J6AU0{l,!>MS2)-a f *{*U0C«,I|»152) + AU0(2#I| 

2M*!2))f6,*AlJO(3,I JM521 . . 

a u 2 ( 2 , i , n*Auo(a,n-Afon ,n-3.*(Ano(3, n«*uo(?, n i 

au2(mu, i#n»*«JO(Msi ,n-Auo(M\,n-3,*(AuofM5o,n-Auo(Mii,ij) 

t> AN(2?3!l|NXMMS>«AUt(2,3,I|NX«M3)-CC*(AU(2,5,I|MX^M3)4AU(2,l,I,MXM 

jH3)-0 t *CAU(2 f «,l |NXMM3)*AU(2,2,1 >NXMH3)U6,*AU(2, 3,IJNXMM3)+AU2{2, 
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23, ! jNXmmj) ) ----- — - . 

Ttl n 1^53) S(5B VGA Twt?f Ain (M 1 , 2, 1 1 NXMR) , INT ( t I N*3 ) f Tf 1 ( 1 1 M5S> ) 

AN(M1 ,2, I f NXMQ)b«J8VSCATR(Tf 1 fl } nS 3) , INTO fN^S) » AN(Ml , 2, I »NX*R) ) 
ANf?,N1,IjM5'Pl*AIJt(2,Nl,TfMs;o) 

AM (2, v ' a , I j M50 ) a AIM (P.NJ a , I , MSO )-CC* ! All (?, N52, 1 t M50 ) . A(J (2. N 1 , I ;M 
J50)-3,A(AUf?.Me;o,r|MSO)-AUf?,Mt<l, I | MSO ) ) ♦ AU2 ( 2, * 1 0 , 1 i *50) ) 

AN (2, N1<5 , : JMS0)BAIJ1 (2.N15, 1 » 0 1 -r C * ( AUf2,N53,I |MS0)-AU(2,M1 1 , 1 f M 

180)«3,*(AU(2«NS1,Ii M S f ')* , Aiif2 > ^l5»Ij”80))+A|i2(2#N15,It M 50)) 

AN(? f 2, ! »MFn^3AUH?,2,I jMSf))-rC*(Alif2,a f i»MS0)-r,M(I ) *AU(2,N12, 

I! jm50).«.*(aii( 2,3,I>MS0) + AI1(2, 1 ,1»M5Pn-»6,*AU(2,2, If M50)+ALI2C2,2. 1 

2fM50)) 

AN(2,wi?,I , m59)bAH 1 C2 f Nl? > l,MHO)-CC*(AU(2,Nl3,I|MSO)«GH(n *AU( 
12,2rI| M 5P)»a.*CAU(2,Mt6,l>M50),GM(i) *AU{2,l,IfM50))+6.*AU(2,Nl? 

2. 1 j M S0)*AU2(?,Nt2,i,HS0n 

AN(?,Nt*, I ,M50 ) a Alii (2,N1 k, I |M50)-CC*(Atl{p,Nl7, I |MSO)«GMU ) *AU( 
*(AU{?,N13,I|MS0)+AU(2,N12 # T|MSO))+6.*AU(2,Ntfe,IfM50) 

2+AU2(2,Ni6, 

22P AU1 (2,2, ! |NXM,>)eANf2,?, I JNXM61 

AiJ(2,N1?,2|MS0).*-Aii(?,Nl 2 # p|M5o) 

AU(2,Nt6,?f w «;n)B-Ai-((2,Mi6,2fM5n) 

DC! 23n IeJ ,0 
!IsI 

IF n »EO,a ) I Ta5 

?30 AU1 (2,1 , 1 |M50)«Allt (2,l,IfM*!0).cC*(Alif2 r l,IIf M50)*»UC2,M16,II|M50). 
10.#(A.U(2,2,!I|^50)*AU(2,N12,n t Mgo))+6.*AUO(2,I»M50)*AU2(2, i # I f 
2M50)) 

AlMU,3,1|N*M7)BQT»VSCAT«(TFfl(1,H5n,INT(i,N5l)|AUlCl,3,l»NXM7)) 

AU1(J,3,?fNXM7)*r)flvaCAT»(TF < »(1jN , ’U,TMTM|NSl)fAiiin,3,2fNXM7)) 

A^l ( 1 # 3, 3f N * M 7)«R.flVSCATB(TFQ( i jM5i j , imt (| |N 5i ) | aim ( 1 ,3, 3j VXM7) ) 

aim n , 3,a r NXM7)aHBvscAT«(Tr < » ( i iNSn , imt ( i f Nbn i aui n , 3, Nxi^rn 

AUl (1 , 3,S|NXM7)*ra8V8CATR(TE10(l jNS1) # INTM j NS 1) » AMI (1 , 3 , 5 | NXM7 ) ) 
TF(L.GT,LTURR)GO TO 231 
CALL VEC3 
GO TO 23? 

231 CALL VFC5 

232 CONTINUE 
CALL VECa 

C At L BOIJNO 

C EDDY VISCOSITY CALCULATION’. 

IFCL.lt. LTU9B1G0 TO 235 

LL»L/2S 

LL*LL*25 

IFCLL.NE'.LIGO TO 450 
CALL OERV 
CALL EDDY 
235 CONTINUE 
LL«L/500 
LL*LL*500 

IFCLL.NE’.LTGO TO «50 

C SURFACE HEATING RATE AND SKIN FRICTION COEFFICIENT CALCULATIONS. 

KRITEC6,a00) L 
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250 

C 

C 


300 


350 

060 

C 

C 


26! 


260 

000 

005 

020 

021 

031 


WR!TE(6,005) 

TF(l»NU)cCH(1 |»1U 
00 250 N*i ,Mf 1 

OSHt»(0,*SH(?,N).5.*SH(l # N).SH(3*N))/ZZ/3(N) 

CH(N)»VIS(1 ,M)*DSH1*XL 

DUl»(<l.*U!2.M)"3,*U!l.N)«U!3,N))/ZZ/S!N)*XLl 

CF(N)*2,*V!SM #N)*fPUl-UM /F M 0*(l,Nn/«C 

WRITE (6. 020) Y(NT,S!N),US(N),CH(N),CF!N),PW!N) 

TH1 CODE C ALC ULATFS ONL* PRESSURE 0 f? A G AT ANGLE OF ATTACK AS A 
FUNCTION OF LOCAL BODY RADIUS normal TO THE 80QY AXIS. 


cPi n ) sPw n ) 

CP?(l)®o'. 

tdpi ( l ) 3 o , 

C0P2(1 )«0, 

DO 300 N3l.N1 a 

CPI (NM)o(Pw(NM) +PW(N*M1 ) ) /2, 

CP2( N *nB(PW(N+l)«PW(N4.M1))/2, 

COPl(N*n=COPl(N)4'.5*DY*{CPl(N)AR{N)*SIN(THE(N)jACPlCN+l)*R(N-U)* 

lSIN(THF(N+n ) ) 

CDP2 (N* n *COP2 (S ) ♦. 5AOY* (CP2 (N ) *« (M) #COBC THE (NJJ6CP2CN* »?*«(»•♦ I )* 
icos(Thfcn+d)) 

COP( 1)02'. *PW ( 1 ) 

DO 350 N e 2 » N 1 

CDP(N)o(a.*CnS(ALPHA)*CDPl fN)-0,*SIN(ALPHA) 3 COP2(NJ/PI)/R(N)**2 

WRlTE(h,Oh0)f«(N) # CDPfN)#NBl,M) 

FORMAT!?*, 'ROnY RADIUS Rs ' ,F 1 2'.S, 1 OX , 'PRESSURE DRAG»',F12.5) 
CONVERGENCE CHECK, CONVERGENCE CRITERION IS THAT OVER 500 TIME STEPS. 
THF HEATTNT, PaTFS SHOULD NOT CHANGE BY MORE THAN l j PERCENT , 

TCI (1 » N 1 1 )SCH( t »Nt U-TEC 1 »N 1 n 

TF1 (l|NinaVAHS(TEt (i | Nil) fTEtMlM.il)) 

TF2( 1 »N 1 1 )srRTT*CH(l»Min 
1 I oQBSGE ( TE 1 (1 >Nin.TF2(l»Nll)) 

IF(II,FQ.NinGO TO 261 

LL«L/5000 

UbII*5000 

TFfLL.NE'.L )G0 TO 050 


CONTINUE 
DO 260 N*1 ,NJ 1 
WRJTE(6,021 )N,Y(N) 

WPlTEC6!«a5)(ZN(M,N),UCH,N),V(H»N),MPCM,N)»PCH, N|,RO(M f N),T(M,Nj f 


1H(M,N), VIST(M, N) , Mai, Ml ) 

CONTINUE 

IFtn.FO’.NinGO TO 551 

FORMAT!//, 10Y, 'NO. OF I TE« AT 1 0N8» ' , 1 5 ) , „„ 

FORMAT!/. IX, # 0IST, ALONG BODY * ,5X , 'SHOCK STANOOFF 01 ST , * , 5X , 'SHOCK 
ISPEED'.SX, 'STANTON NO. ' , 5X , 'SK IN FRICT ION' ,5X , 'WALL ‘ P*. # ) 

FORMAT !5X ,F 1 n, 5, 1 OX. FI 0,5. 1 OX, FI 0,5, 6X, FI 2,8, 6X, FI 1 ,7«6X ,F 1C ),5) 
FORMAT! /,5X, 'BOOY STATION NO . ■ ' , I 3, 5X, 'OT ST , ALONG BOOY«' , F 1 0 ,5. / ) 

F0RMAT!/,6X,'ZN',9X,*U',9X,'V',8X,'WP',<»X,'P',8X,'R0',9X, 'T',9X, 


1'H',RX,'VIST'} 
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NXM8*nx m -2*mi*j 

NXM9eNXM«a*Mj + l 

NXMtOafl*NXM 

NXHM8NXM«2*m 4| 

MXMM2aNXM-2*Mt«a 
NXMM3aNXH.a*Ht-2 
C STEP STZES 

DZe-J ,/mi 1 

xxs2,*nv 

7Zs2.*0Z 

YO »M )=f38VTNTL(n.,DV|Yf 1 »N 1 )) 
Z(ljMn-f)flV|WTLM.,DZtZ(1|Mn) 

*1oPI/2,-Thc 
YMoO Y*w 1 a 

RBa( (VM- A n*SI^(THC^COS(TMC))*RN 

C NJi IS THE “ESH POIN” UPSTREAM OF JUNTURE POINT. 

NJaAt/DY 
NJlaNj4l 
NJ25NJ.J 

NJJsMJt4J 

C INITIAL SHOCK shape ANO SHOCK SLOPE. 

S(1»NJna,17 
00 17 nsniJ3.ni 

17 S(H)aS(M.i UOYaSTNfTHESH) 

C BODY r.EOMETBY, 

00 10 Nat.Nl 
JF(Y(N),GT,AUGO TO II 
THE(N)aPI/2..V(N) 

CUR ( N ) a | , 

R(N)«COS(THE(NJ) 

CO TO 10 
11 THE (N ) aTHC 
CUR(N)sO. 

R(N)aCOS(THC)A(Y(M)-Al )#SlNfTHC) 

10 CONTINUE / 

DO 20 NaJ.Nl 
VS(N)aCOS(THE(NJ) 

20 VN(N)a.SIN(THE(N)) 

DO 35 Nal.Ni 
S1(1.N|M1)bS(N) 

CI»Rlfl,N,Mn*CURCN) 

R1 CI#N|M1 )3B(N) 

J5 THEl(l,N»MljaTHE(N) 

THEC(1#1 jNXMleVCOSITHEl (1,1 »NXM)|THEC(l,t |NXM)) 
THFSU.1 }NXM)sVSTN(THEl(l.l|NXM)}TH£$(l.t|NXM)) 

INT( 1 |N1 )aOSVINTL C 1*M| f INTM fNl ) ) 

C CALCULATION OF THE STARTING SOLUTION. 

00 00 Nal.Ni 

40 TFHP(N)aS!NrTHE(N))**2*(,94,!*fN«! )/N| ) 

P(l,llNXM)*OSVSCATRrTEHP(l|Nl ),INT(1|N1 ) |P ( 1 , 1 t N*N) ) 
Tf hPC ! jNi )aO. 



ro«M*Tnx,9Pio,5) 

A50 l«l*l 

IF (L.LE ,L M AX )f»0 TO 1 
551 CONTINUE 
RETURN 
FNO 
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SUBROUTINE BH VP 

DIMENSION CPI (15),CP2(15),rnPl(15),CDP2(l5),CDP(15) 

DIMENSION Yf20),THEr2<n,F»(j»q),cHC20).TeMP!(2«M. T EMP2non,Znon. 
iTznon,Tzu mn # r,Mfs) 

DIMENSION TAlif?q,5),TAM(2P,5),TANH20,5),TAN2(29,S),TA0(29,5), 
1TAU1(29,5)»CE(29) * 

DIMENSION Tf*(29) ,TF 9129) ,TE1 0(29) 

COMMON/E /ThEl f 101 ,29), THE? (t 01,20), VAIRfl 01,29) # DT( 101,29) 

COMMON /F 1/N1,Nlt,»J12,N)6,N15,M1,M11 # M50;Nia 

C0MM0N/F3/GMi,r,M?,GM3,GMii,r,M5,r;M6,GM7,GMfl,GM9,B2,BE,TW,GAMA,SIG,MM 

C0MM0N/F2/NXM,NXM1,MX M 3»NXM/4,NXMB,NXM5,NXM61 

COMMON/Fi)/XX,ZZ,FM,HETA,SIGT 

C0MM0N/F5/SC?O),sS(?9),R(29),CllR(29)»USf29),DS(29),VS(29),VN(29) 
C0MM0N/Fh/INTC?9) ,TIMP(?0 ),tEMP(20),tE(29),TE1 C29),TE2(29) , TE3(29) 

l,TE«(29),TES(29),TEfc(?9),Tf7(29) J ' (<f 1 

COMM0N/E7/U(101,29),V(lO),?9),T(101,29),P(tfll,29),RO(101,29), 

1 WP( 1 01 ,20) , VTStl 01 , ?9) ,Hfi 01 ,29) ,SH( 1 01 ,20) f PW(29) , VISTfl 01 ,29) 
C(1MM0N/Fn/HT(J01,?O),EM0A(trt1,29),7N(i01,29) 

1 2 O) M0N/|r< ’ /S1 n01,?9) * SS1(101 » 29, * r) 3inoi,29),CURinol,29),Rinoi, 

COMMON /F 1 O/THEC ( 101 ,29) ,THES( 1 01 ,?9) ,TV1 (1 0 1 , 29 ) , T V2 (J 0 1 , 29 ) 

COM MON/P J 1 /Dl.lf 101, 20) ,OSH( 1 0 J ,29) ,Dwp(l01 ,29), A6{ 1 01 ,29) , 
1AT(10I,29),ABU01,29) ' 

COMMON/E 1?/T7?(10 1,?0),ZJ (1 0 1,29), TZ3(101,29),TZ<|( 101,29) 

COMMON/E 1 3/AUH101, 29, 5), AU2( 101, 29, G) 

COMMON/E 1 «/All ( 1 01, 29,5), AM 1 1 01,29, 5), AH (1 01,29,5), AQ( 101,29, 5) 

C0MM0N/E15/AU0U01 ,il),AMO(J 01,2,0), ANO(101,a),AGf)M01,«) 
CnMM0N/F16/HTNF,REl ,RE2 
Vise 1 » 1 12929)00, 

VIST ( 1 , 1 j2920)ao. 

AU( 1 , 1 , 1 1 1 0605)80, 
am ( 1 # 1 , 1 ,10NO«5)30, 

ANJl , 1 , 1 | 10609)*0, 

*0(1,1, 111 Uhfl5)s0, 

AU1 (1,1 ,1 |10MS)«0, 

AU2(l,l,1 t 14(fcfl5)«0. 

CH( 1 |29)s0, 

CONSTANTS AND EREESTREAM CONDITIONS 

P|Bfl,*ATANCr.) 

GAMAS) ,0 

RG4S«2B7‘. 

SIGo.72 

CPbG*MA*RGAS/(GAMA«1.) 

SIGTb,9 
CCb,001 
EMbIO, J 
TEb06,26 
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nono 


'pFbIOO.77 * 

R0F=PF/TF/R0AS 
AUPH*s5,*0T /i 00, 

THCs«5 J .«P!/lB0, 

TWb 330,6 
TWsTW/TF 
RNs,03l 7S 
DYb. 18«5 
BETAS!. t 
F0T=1 , 

THE3Hs?,tPT n BO. 

VFbFM* ( GAMA*RGAS*T*)**.5 

VISFal.u«iR/in,**h*TF**t,5/<TF*110.) 

REbVFaOOFaRM/VISF 

C L*AX ! 3 THE MAXIMUM NUMBER OF TIME STEPS AFTER WHICH THE PROGRAM 

C WILL STOP EVEN THOUGH THg SOLUTION M 4 Y NOT BE CONVERGED, 

LMAXB20000 
LTURR*2Sono 

FOR LAMINAR FLOW, CALI. VEC3 and FOR TURBULENT FLOW, CALL VECS 
LTIJRB IS THE NUMBER OF ITERATIONS UP TO WHICH THE FLOW RE m A INS 
LAMINAR, 

CONVERGENCE CRITERION. 

CRITb.oi 
DO 5 lei ,5 
5 GM(I)Bf.i ,)**T 

Gmib1,/GAMA/f m **2 
r.M2*2,*GAMA/ (1 ,*GAMA) 

GM JbCPaTF /VF** 2 
GMflsSlN(AlPHA) 

G M B*1 10./TF 
GM6*1 ,+GMS 

GM7sCGama»! ,)/(GAma*1 , ) 

GM«b U .+GAMA)/(?.*GAMa*FM**2) 

GMRbCGAMa-1 , ) /GAMA 

GMJ osTw*r,MVGAHA 

HINFb.PRS*C.S>GM3) 

RF1bSORT(RE1 

REPa0.0168*RF 

C MFSH SIZES AND VECTOR LENGTHS 

NIb15 

Wit b2*N 1 ■ 1 

N123NU1 

NIUsNl-l 

NISaNI 1-1 

N50BN1-2 

N513N11-? 

NlSsNi+J 

N17BN1+U 

N52*Nl-3 

NS3SM1-3 

N16»nj+2 
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Ml * J 01 " * ' 

MU*Ml«t 

M508M1.? 

MS^aMl.u 
N X M a N 1 1 *M1 
NXMlBNXM.Ml \ 

NXM2bKxm. i 

NXMJbNXm.2 

NXMilCNXM.Ml 

NXM5=^1 «*mi 

NXM(,aNSl*Ml.g 

NXMbUNSl**!,.! 

►JXH7SNXM. J*M 1 + 1 
N X M 8 B N X M m M 1 . j 
NXMOaNXMadAMf +\ 

N X M 1 0bS*n;XM 
NXMMbNK«.P«mi ♦! 

NX MM lsSdMXMg 
NXMM2sNXM-2*Ml.a 
NXMM3jNXM«a*Mt-2 
C STF.P SI7FS 

DZb-I./mh 

XXa2,*0Y 

ZZn2,*0Z 

Y(t fNllaOftVINTL (0,,0Y|Y(1 |N 1)1 
7(11*1 )BQHVtNTLtl . ,DZ|Z(1 J-M5) 

C INITIAL SHOCK SHAPE AMD SHOCK SLOPE. 

8(1|N11)b,17 
DO 15 ns3,ni 

15 S(N)*S(N»1 )+r>Y*3IN(THESH) 

SS(1 )«(Sf2)-SfNl2))/xx 

DO 16 N a 2 » N 1 u 

16 ssfN)BfS(N*n»s(N«i ) )/xx 
ss(Mi2)ars(Ni6)-snn/xx 
DO 17 NaN l 6 # N 1 5 

17 SS(N)a(S(Ntl)*S(N«t))/xx 

00 18 N B M 1 ,N1 1 ,N10 

18 SS(N)a-fu.*SfW.i )-3.*3fM)-S(N-2) )/X« 

C 900 Y GEOMETRY, 

Ct'R(t)st'. 

THE(l)aPl/2, 

R(1)=0. 

Aal ,/TAN(THC)**2 

BaSORT (A ) 

00 10 Nb2.N1 

F1»0Y*S0RT((B**2*B(N-1 )**2)/fB**2+r r.*rA/B)**2)*R(N-!)**2) ) 
F2b 0Y*S0RT ( (B**2*(R(N.n+Fl/2,)*B2)/(B**2*(l, + (A/B)**2)* 
l{R(N-lUFi/2. )**2)) 

F3b 0Y*SORT ( (B**2+(H(N»1)*F2/2,)»*21/fBA*2+(l,+(A/R)**2) 
UfR(N-U*F2/2,)*a2)i 



' >‘««OY*SOBT( (p**2*(R(N-lH>31**2)/<B**2*n,f(A/H)**2>'*CR(N-l >*F3) 
1**?)) 

R(N)sP{N.n*rFt + 2,*P2*2,*P3*F<n/'6, 

XsA*SURT(R**?+R(N)**2I/H 

CUR(H)sA**a*fl**a/(A**a*R(N)**2*B**o*X«*2)*M,3 
XeB*X/( ASSORT (X#*2-A«*2)l 
10 THE(N)=AT4N(X) 

Y(Nl2jMitnsvr2jNia) 

R ( n 1 2 j n 1 0 ) = R ( 2 t *J t a ) 

THE(Nl?jNiaisTHFf2jN1U) 

CtJR(M!2|Nia)sCURf2|Ml«) 

on ?o ns i i mi 

VS(N)=cnsCTHE(N)-ALPRA) 

20 VW(M)aSlNf ALPHA-THEfN)) 

DO 25 NcN 1 2 * N 1 1 
VS(N)sCOS(At.PHA+THE(Nl) 

25 vN(N)e-STNtALPH4*THE(N)) 

00 35 Ns 1 # N 1 1 
81(lfN|MnsS(N) 

cum Cl .njmdsCurcn) 

ri n #n,mi)sr(n) 

35 THE! ( 1 ,NfMn*THE(N) 

THfC(l , 1 |NXM lt VCOS(THEl (1# 1 iTHECd * I iNXM) ) 

THE3( 1 « 1 |NXM}sVS1N(TWE 1 ( 1 r 1 1 N X H ) |THp«5f 1 , 1 |NXM) ) 

B1s(BET44l t )/CBETA»l,) 

B2*2.*RETA/AL0G(B1 ) 

XL1 3 -B2/(«FTA**2-1 , ) 

XLb-2.*R?/ (HPT A **2-1 ,1/SIR/RE 

TZ(J »Mi>aRl**fZ{ 1 t*n.l , 

TZ1 (1 |MnsBETA*(TZC1 l M t )-l.)/(TZ(l J M 1 > ♦ 1 .) 

00 30 Nsl ,N1 1 

so zi h #N|Mn* T zt (i i M n 

ZNfi#i|Nx«)«sici#t »nxh)a(i .-ztn.ifNxnn 

END A ( 1 , 1 ,NXM)«1 , + ZN(t ,1 |NXM1#C«IR1 (1 , 1 JNXM1 

btm,iimxmi.»ici,i,a.xm)*zncwi»nxm,*thfcciuinxm,- 

T?2(i * t |NX M )eBETA**2.»Zl ( 1 * 1 jNX M )**2 

7J3C1,1 *KXM3a2.*Zl l 1 • 1 

TZ«U # l jnxm)sR2/tz?( l , \ jnxmj 

JNT(1 |N 1 nsOflVIMTL ( 1 ,M1 j!NT( 1 1^1 n ). 

C CALCULATION OF THE STARTING SOLUTION. 

00 00 Nct.Nt 

40 TEMP (N)sSlN (THE (N) -ALPHA )**?*( .<>♦. l*(N«n/Nl 1 

00 as nsn 1 2 « n 1 1 , ; 

45 TFMP(N)«SIN(THE(N)+ALPHA>**?*(,9->.l*rN«NU/Nn 

P(l.l»NXM)«QBVSCATB(TEHP(l»N|n,INT(1»N1|)|Pll,l»MXM)) 

temp ( i i N i i)*o, . 

U(l,l»NXM)sOBVSCATR(TEMPClfNin,lNT(UNin jU(l,l|NXMJ) 

V(1,l lNXM)sOBVSCATR(TEMP(l|Min,INT(l »N1 1) 1 V ( 1 , 1 |NXM) 1 

HP ( 1 * l »NXM JsOBVSC ATR (TEMP ( 1 |N1 1 J f INT C 1 1 NI 1 ) I HP ( 1 » 1 1 A*XM) J 

TEMP ( l | M 1 1 ) STW 

T(1 # 1 |NXw )nQftVSCATRf tFHP M|MU*INTftjNi1)|Tn,l I^XM) ) 


2h 



tfsoini n*ssM }N!n/r i ,+sr i iNin*rup(i »ni in 

TF6d im t n = r. + tfs(i imi n * te 5( i j mi n 

TF 6 ( 1 f N 1 1 ) = VSORT f TFfc M t n t r 1 nil 1 n 

"5 71* * ! |V !)% f VN! I :*•! m»v S ' « 1* 1 1 )*Tfc5 ( 1 |N1 1 } ) /TEMl I Nil ) 

?" ( 1 | ' i i ) 3 ' T E 7 f 1 | N 1 1 ) * T F- 7 ; \ .m J 1 "> / 1 

Pint » i »Nyf.n s'j«v sc*TR(Te><pm rn i 1 ) # INT( i ini i ) f p (»i , t »NK*;t j ) 

TEMP' t , ' 1 ! ) =T(;mp( i |>t|i)/r,M| 

T M » Ml ! ) n ( tc hp n , n » n > r iH7 l / ( { *TfMpf»iNi i )*G« ) 
Trin;Mnsri,«t./mii'Mij)*Tr(iiNin/TEf>(iiNin 

te? n nji-sd imi ) - 1 e i r 1 i mi n *t? 5 f i > n t n 
TE3d I Nil ) rVf ' t Mt D + TEM 1 | hi ! ) 

THMiVIDsTFMOf.^iiD/TFfi^nt) 

l»(n i # 1 t M* M n = Q« vst AT W ( Tf-',? ( 1 I h 11 ) » T'lT n t N1 t ) fll - M 1 , 1 ,NX*M n 
i/(Hl,lf^Mi> 5 r)8vsCAT«fTE3M«ilin»TNTM»Ml)TV(Ml # i|NXMin 

7 (Ml * 1 I NXM1) oTPVSr ATRf Ti 0 ( ) |N1 n , I NT ( 1 jN) 1 ) J T (Ml , 1 ,NXM1) ) 

00 61 N=J,Mi 

OOFLs- (S (M-S (N*mi a n /?„ 

61 WPf*l ,N)s-Gm«.TF. 1 (N)*OI)F.L/nT(Ml ,Ni 

00 65 NSM12,M1 

DOELP ( S(N.\it<l)«S(Nn/«*. 

65 MP(Mi,Nisr,M«.TEt f Ni*onEL/nT ( mi # u) 
wp i » n =-ur » i , i ) 
on 50 n=i ,ni i 

0VII3 (U( M 1 ,N).lJ (1 ,N) )/mj 1 

DXVe(V(Mi *N))/mh 

0XP= (P (M I 

OXTs ( T (Ml ,N)-T ( 1 # Nn/**l 1 

0XWPo(u(P(M1 ( M),WP(1 ,1111/Ml 1 

U( 1 »Nf Nl 1 1 S 0 AVTNT 1 . HIM »N) ,OXl!ftl( 1 ,N|Mt t n 
V(l#NjMinP0AVTNTLCWn,H>«OXVlV(l,f!;Mll)) 

Pf 1 ,N |H1 1 ISQBVTNTL (P f 1 ,N ) , pup JP ( 1 ,1 ; M 1 1) 1 
T ( 1 * N j M J i noPVTNTL (T ( 1 , f ) ,0X7 |T{ 1,NjMll)) 

50 Wpfl ,N,M 1 11 sO«VTNTL (wPfl , m) , pxwP }'-P( 1 # KJ| Ml 1 ) j 
Ron , i |Nxm) 3 p ( i , i tN.xMj/T n , 1 1 nxmj /rmi 
8H( 1 , 1 |N*M) S T ( 1 , 1 |NXM1*GM3 

H( 1 , 1 »Mv )sSwM , 1 tNXMna.'U * 1 »MXH)**?*V{ l , I , fy XM 3 **21/2, 

Tv2n*i»NXN)svsoRT(T(i,irNXP)|Tv2n f i»NXMn 

V!S(1 .1 »NXMlsTf 1 , 1|NXN)*TV2(1 , t JNXM)*GH6/(T( 1 , 1 |NXM)*GM5i 
WRITE f 6,309) ALPHA, RE, EM 

399 FORMAT(/ # lX,»ALPHA , t Eft, 5»5X,* REYN NO . * , F l 5 .3 , 5X , »M ACH N0,*,FB,3,/) 
WRITE (6, 035) f Y(N) , THE(N) ,R(N) ,CUR(N) , Nel ,N1 1 ) 

935 FORMAT (2x,'Ys',Flfl.5,SX, 'THEs * , F 1 0 , 5 , 5X, 'R = ' , FI 0 .5 , SX , »C URc » f F 1 0 , 5 
1) 

DO 60 Nal ,NJ 1 ,2 
WRITE(6,021 )N, Y(N) 

WRITE(6,U30) 

WRITE (6, 000) ( 21 (M,N),U(M,N),V(M,N),wP(M,N),P(M,N) ,R0 (M,n 

n,T( M »N),H(M,N),M«1,Ml) 

60 CONTINUE 

030 FORMAT!/, 6X,»Zl»,9X,»U»,9X,»V',ftX,»WP*,9X,'P',«X,»R0', 

19X,*T , ,9X# , H*) 
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440 


C 

l 

C 


70 


C 


74 

75 

76 


"format ( i X,~8M n,5) 

D M *»1 ./u’./ H l 1 

TFMP2( 1 |M1 noQRVTNTLC 1 ,+RM.OMf TEMP2( 1 |M1 1) ) 

AU 1 Cl . 1 , 1 >NXMt OTal , 

TF8C1 JN51 ) = 1 ' f 
TER(1 1*51 )a0. 

TEion iNsn=r. u io 

Lsl 

MARCHING JM time starts HERE. 

CONTINUE 

MMal 

LOCAL time S’FP calculation. 

VAIR(l#lj M *' <: M o GTMA-''' 5 n,l)NXM2)/R0n,!|NX' < 2) 

V A I R ( 1 t 1 jHXM->)T!V'SORT IVATRCl*l|NXM?),vAIR(l,1|NXH2n 

0T(1 , 1 :NXM?) = F0T*C7.Nf?. 1 jNyMgJ.ZNn , 1 | NXM2 ) ) / < V ABS ( V ( 1 , 1 » NXM2) J 

1THF2C1 , 1 f NX«?)UVA TRft * 1 »NXM2n 
no 70 Mai ,M 1 

OTC 1 ,M *-1 UsnTCl 1>*TFMP?C1 | M 11) 

TIMPM , NU )=n B VGATHPfOTfM11 ,1 »NXMn , INT(1 |N1 n ,TTMPn,Nl 1)) 

OT(Ml f l|MXMM««nvse ATR(T?MP( iNTtljNin JOT (**1,1 »NXM1J) 

PREOICTOR STEP. 

SSC1)=CS(2)-S(N12))/XX 

DO 7 a Na2.Mii 
ss(N)«(S(N*n«scN.m/xx x 
SS(N12)a<S(Ml6)-8(t ))/XX 
no 75 NcNl , Ml l .N1 a 

$S(N)a-f u,*S(N»n-3.*S(M«S(N«21)/XX 
00 76 NaN 1 6 . N 1 5 
S3(M)a(SCN+l)«SCN«l 11/XX 
CALL SHOC< 

CALL DE R V 
CALL VEC 1 

00 110 Tal.S . .. 

AUlC2.2,TjNXMfO3All(2,2,I|NX M 61-0T(2 f ?iNXMfc)*(CAM(2,3,Il fJ XM6)-AM(2, 

l ? .J,NXM 6 M/DV*(AN{3,?,I»NXM6T-AN(2,?,T|NXM6n/OZ + AO{2,2, I jNXMfel) 

AUl(2.Nl.II M 50)aAU(?,Ml,ItMSO)-OT(2 r Nl|MSO)*CCAM(2,Nl,I|M5t.)-AMC2, 

IN 14, I | M 50) ) / O T ♦ (AN(3, M 1 . I l M 50)»AN ( ? , N 1 , 1 f M 50) ) /07+A0C2.N 1 ,11*50)) 
AU1 (2,N1 1 , 1 jmSO 1s 4U(2,N| 1 , 1 |M<50)"OT { ? , Ml 1 f M50 ) * ( C AM(2*Nl 1 , I f M50 ) » 

lAM(2#Nl5,IjM50))/DY*(AN(3,Nll,|jM50)«AN(2,Ml,!|M50))/D2+A0(2»Nll, 


TAUC 1 , I |N1 naOflVGATHRC AU(M1 , 1 ,1 jNXMl) ,INT(1 JNI 1) »TAUC1 , I |M 1 ) 1 
TAM( 1 ,1 |N1 naQBVGATMRCAMCMl , 1 , TfNXMl) , IMTtl |N1 1 ) |TAM( 1 , T |M l ) ) 

TAN 1 (1 .1 |N1 l>aOflVGATHR(AN(Ml , 1 , I J NX*M ) , I NT C 1 » N 1 1) | T AN 1 ( 1 , I | N 1 1 )) 
TAN2C1 , I|N 1 na04VGATHB(AN(MU , t ,ItMXM 1) , INTC 1 »N1 1> |TAN2(l, I |N11)) 
TA0C1 , 1 iM 1 )sOflVGATHB(A(5(Mi, 1 , I j Nxm n , I NT ( 1 I N 1 1 ) | T AO ( 1 , I » N 1 1 ) ) 

TAU 1 (2,I»N15)aTAiJf 2,I|N15)-TIMDC2jMt;T*( CTAM(3, 1 ,M|5)»TAM(2,I|N15) 
i)/OY*(TANl(2,I|Nl5)-TAN?(2,T|*M5))/OZ*TAO(?,T|Nl5)) 

TAUl(Nl,I)sTAU(Nl,I)-T!MP(Nl)*((TAM(Ml,n»TAM(Nla,n)/DY*(TANl(Nl, 


1 1 4»T AN2(Nl , I) 1 /02*TA0(N1 , II ) 

TAUl(Nll,I)3TAU(Nll.n-TJMPCNin*((TAM(Nll,n-TAM(Nl5,n)/DT*(TANl 

i (ni i , u-tan2(ni i , i n/o2+TAo(Ni i , m 
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110 AUl (M1,2,T|NXMM).QRVSCAT«CTAU1 (?, I > *15) , TNT ( 1 »Nl 5) » AUl ( HI , 2, 1 | 

INXMM) ) 

CALL VEC2 
DO ISO 1*1*0 

AUl ( 2 , 1 * T l M 50)»*UO(?, 7 | M 501»DT ( 2 , t » «S0 ) * ( ( AMQ ( 2, ?, T I *"50 ) - AHO (2, 1 , 1 
1 |M50) ) /Dy+( anOCJ, I |M«50)-AMn(?# t > ) ) /DZ + AQO C2 , I 1*50) ) 

ISO AUl CMl,l.n.AU0Ml,n-DT(M,,U*((AMOCM, # 2,n.AMOCM 1 # l,in/0Y + CAN 0 ( 

iMi.n-ANoo'i i . n )/nz+AMP(*i * in 

AIM(1,3,1 |NXM7)sOBVSC AWTFAC1, NS1), TNT (1,N51)| AUl (1,3,1 |NXM7)) 

AUin,3.?|NYH7)=OrtVSrATR(TFQ(1 ,^5n*T^T(1 jNSmAUtn,3,2»MXM7n 

AUl ( 1 • 3,3 jNXM 7)*DAVSCATR(TFR( 1 t N5 1 1 * I NT ( 1 j *5 1 1 I AU 1 ( 1 * S*3|NX M 7) ) 

AUl (1,3,0 |NXM7)*G8VSC ATP(TE9(1 |NSl),INT(1 |N51)»Ai. 1(1,3, U»NX“7)) 

AUl (1 *3,Sl N x>*7)a<3BVSCATR(TElO(l »N51),1NT(1 ,MS1) | AUl ( 1 ,3,5|NXM7) ) 

IF a.GT.LTURfllGn TO ISO 
CALL VFC3 
GO TO 15S 
ISO CALL VEC5 
ISS CONTINUE 
CALL VFCO 
CALL SOUND 

'C corrector step. 

DO 160 1 a I * S 

160 AU2f2.?*T|NXMfl)BAU(?,2 # !|NXM«)*AUt(2-,2,T|NXM8) 

DO 170 To 1,0 

170 AU2(2, 1 , T I m 1 nsAU0(2*I|Ml n + AUl (2, 1 ,!|M1 1) 

MMo? 

ssr i>«(Sf2)*scNi2n/*x 

DO 170 N a 2 . N 1 0 . 

170 ss(N)»(8TNAn»s(Af«n)/xx 

SS(M2)3(S(N16)-S(1)1/XX 

DO 175 NaM.Nl 1 ,nio 

175 sS(N)a-(0,*S(N-n-3,*S{fJT«>S(N-2))/XX 
DO 176 N*N 18, N1 5 

176 SS(N)«(S(N*n»S(N»l))/XX 

call shock 
CALL OERV 
CALL VEC1 
DO 160 Ial.5 

AUl (2*2* T (NXMfl)art # 5* ( AU2(?*2» 1 |NXHR)»0T (?,2 jNXM 8)*( C AH (2,2, 1 j N X *6 ) 
1«AM(2, 1 ,!|NXM8n/0Y + (AN(2,2,I|NXM8).AN(l ,2, I |NXM8) )/DZ*A(3(?,2, I» 
2NXM8))) 

AM( 1 , 1 , I |M!)ar,M(I) *AH(1,1 # T|H1) 

160 AUl (2,N1?,I|Hin«,5*(AU2C2,N12,I »Min«DT(2,N12»Min*((AM(2,N12,l| 
lMil)«AH(2,l,T|Hil))/DY+(AHT2,N12,Ii M 11)»AN(l,Nl2,l|Mll))/DZ+AU(2, 
2Ni2,itMinn 

WP(1#N12|M1>««WPC1 ,N12|H1) 

U(1»N12|M1)b.U(1,N12|M1) 

CALL VfC2 
DO 190 1*1,0 

190 AUl (2, 1 , T » H t ni,S*( AU2(2, 1 , ! fHH )-0T f 2, l jMi 1>*( ( amo( 2, 1 , I f *1 1 )»AMO 
if 2 ,?.i»M}in /ov*(ANOf 2, 1 1 mi n.Auoj i , i |Mi n )/dz*aqo( 2 , i jhi m> 
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c '■ 'fourth' OOF'R DAMPING, CC IS THE DAMPING COEFFICIENT t 
DO 200 Isl.5 

AU 2 CJ# 2 «I»NyMM 2 JsAUCS, 2 .IfNKHM 2 )*Aiin, 2 ,i f NXMM 2 ).fl,*(AU(«, 2 ,I»N*HM 

t2)+AU(2,2,IjNVMM?n*h.*An(3,2,T>Ny MM ?I 

TEtCl»NS1)»nAVGATMR(AU(1#2.I|NXM7I,INT(t»N51)|TElC1tN5n) 

TERM »N51 )sn»VGATHR(AU(?,2,I|NXM7),tNTC1 »N5n |TE2( 1 |NS1J) 

TE3 C 1 >NS1> = r 'flVGMM»( Ml { 3 • 2 « T s nxm7 } . TNT ( 1 jNSl) t TE 3 ( 1 I NS 1 II 
TE« M in^I Asf?.pvr,ATMR{ H'(U,2, 1 f MXM71 , I NT M ,NS1) |TE«n 1*511 ) 

yf 5 ■ 1 ( j 5 jAt?pa( 1 jNS 1 I»T^'!( 1 |N 5 i)« 3 ,*f T F 3 (l»N 51 )»TF 2 il|N 5 i)I 
•■.•'i? '2.?il |N*M7)sQRVSr ATRCTfSt'i |NS11 , T NT C 1 j m5 I ' | AH? (2,2, I |NXM71) 

TF1 U?N5«)s'lAvGATMHfAIJ(MSl,2»IlNXM71»I'‘T(1|N5nfTfl(1|NSl)) 

?£?( 1 |N , in=n«VGATHPf AM c HSO # , I . IMTf 1 |MSt) |TF?( 1 | MSI n 

TE 3 Cl|NSl) 3 nPVGATHHfAll(M 11 # 2 *Ir ,J XH 7 ),INT(t»NSn»': 3 (l»NSni 
TF«U ' *51 1 * 0 A VGA THR ( All f Ml I |NX *71 , I*T( 1 jMSn iTF'JM 1*51 n 

TESd )STF 1 M |NSn*TFa( 1 |M«;n» 3 .* CTF 2 ( 1 |NSn«TE 3 (l |N 5 l) ) 

200 Al'2( M l 1 1 iNXM7)sQPVSCATHf Tf 5(1 > *5 11 , I*T ( 1 | N51 ) I AU2'*1 1 # 2 . I f N X M 7 ) 


n 

on 210 ini ,a 

Ali2(J.1,Tl M '5?laA t '0(5,T|M52l+Alin(1,I»M52)««.MAUO(<j,I|M5?n.AUO(2,I| 

?M5? 1 /♦6,*A | J0(3» I J U S?) 

ah? ( 2 , 1 , nsAiiou, n-Ai'n . < ,ivs.*c»"oc3,n-AL'nr?,iii 

210 AU2 ( M M # 1 # 1 1 sAUO ( M S1 # I ) "Alin (Ml , J ) «3 , * ( Alin f H50 , 1 ) •AilO ( ^ 1 1 » I ) ) 

DO 220 1*1 ,S 

AN(2,3,T »n*m m 31*MII <2,3, 1 |NX*H31-CC*( au<2.5, I »nxmm3i*aii(2, 1 . 1 1 N X M 

lM 3 )-tt.*(AII( 2 ,a,IjNXMM 3 )*AU(?, 2 ,I jNXMM 3 n*O t *AU( 2 , 3 .I|NXMMJUAI) 2 ( 2 , 
? 3 ,ljNXMM 3 i) 

TE1 { 1 |NS3)sOAVRaTH 9(AH1 ( M 1 , ? # I f * X m Q 1 , XMT (1 | S 3 3 | T E 1 (1 f M 5 3 3 1 

AN(N| f 2,I|N»MQ)sl3gVSCATRCTF1 (l|N5n,tNT(l|N«;3),AN(Ml,2 # I|NXM9)) 

AN(2,*i , T t Mt?ft)sAUl (?.*1 , J JMSO) 

AN( 2 »N 1 (I, I |«SO)sAUH?»Nltt,I |MS0)-CC*(AU(?,N«j2,I|Msn)-AU(P,Nl # IfM 

1S01-3, *(AU(?,N50, I I M S0)»AI|(2,N1«> I |MSO) )*AU2(2,Nm, 1 JMS01 ) 

AN(2 f Nis,I: M 50laAH1 (?,*15,I »*5n)-CC*( AU(2,N53 . I »mS 01-AU(2.N1 l # 1 1 M 
150).3.*(U<:?,*51. 1 |MS01-AU(2,N1S,I »MS0n+AU2(2,MlS,T|M50)) 

AN (2, 2,11*50 Is AIM (?,?.I»MS0i-cr*(AU(?,a.I |M501-GM(I) *Ail(?,N 1 2, 

1IjMS0)»u.* ( A 1 1 ( 2 # 3 . t f M 50) + *l | (2» 1 , I j '• S 0 ) )+6,jrAU(2#2, 1 |MSO)+AU 2(2*2* I 
2 :* 50 )) 

AN ( 2 . N 1 2 « I j'mSOIoAUI (2, N 1 2 » T|M50)»CC*(AI.I(2#N13»I | M 50)«G*(I) *AU( 

12# 2* 1 1 M 50)«'j'. * ( A11(2,*1 6. 1 |MSO)«r.M ( I) *A'JC2.1#I|N50))*6.*AU(2#N12 

2,It M 501*AU2(?#Nl2,I jMSO)) 

AN( 2 iNl 6 #Il Mt i 0 ls AI M( 2 ,N 1 h,IjMr, 0 )»rC*(At | ( 2 »N 17 ,I|M 50 )-GMCn *AU( 

12.1.1 l M SO)-u'. *( AiJ ( 2, « i 3, 1 1 M50 )♦ All ( 2, n 1 2, 1 1 M50 n+t..* AU C 2. N1 6, 1»M50) 

2*AU2{2,*lfc. 1 5 HSO ) 1 
220 AUt(2#2,I |NXMft)aAN(?.2.I|NXMb) 

AU(2.N1 2 , 2 (Mso)a-AtJ(?.N 12 ,?fM«;o) 

AtJ(2.N16,2i M S0)a«AII(2,Nlfe.2»M5Ol 

00 230 lei ,« 

1 Is 1 

230 AU1 (2. 1 * I f *50 ) a AIJ 1 (2. 1 . 1 jM^0)-CC*(AU(2»3. II | M 501 + AIK2#N16. II | M 50)» 
10.*(AUC2,2, II jM«J0)*Ai.l(2.N12,II |M50))+6.*AUO(2,I»M«ioUAU2(2, 1, I 1 
2M501) 
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*1)1 (1 ,3, 1 |NXM7)«QflVSCATR(TF«(l |N5l), Wfl » N5 1) « *U 1 f I 3 , 1 | N X»7 ) ) 

A ' J 1 ( 1 »3,2tNX»'7),«8VSC*TR(TF<?M j NSl ) , JMT (1 | N5 1 ) | *1)1 ( 1 « 3 , 2 1 NXM7 ) ) 

*U1 (1,3,3 ) N * M 71an<1VSC ATR(TF0(1 I N51 ), JNT(1 |N5l) f*U1 ( 1 , 3, 3 I NXM? I ) 

*UI (l#3,«»fcXM7)B08VSC*TH(Tf P(i ,N51 ),INT(1 |»<SI) jAIII (1 ,3,4»NXM7) ) 

*U1 (l,3,SlNXM7)B0«1VSC*TR(Ttl0(l)N5n # lNT(l»N5n»*ni(l,3,S|NXM7)) 
ira.GT.LTURRJCO TO 231 
r ALt VFC3 
CO TO 232 

231 C*IL VPCS 

232 CONTINiif 
C*LL VFCU 
CALI BOUND 

C EDDY VISCOSITY CALCULATION, 

IE(L,LT,LTUPR)G0 TO 235 

LLbL/25 

LI.bLL*25 

JF(LL,NF.L)G0 TO <150 
CALL OERv 
CALL EDDY 
235 CONTINUE 
LLbL/500 
LLbLL*500 

IF (LL.Ng’.l 1GO TO U50 

C SURFACE HEATING RATE AND SKIN FRICTION COEFFICIENT CALCULATIONS. 

WRITE (6»«00) L 
WRITE (6,<J05) 
tem »Nin«CH(i *Nin 
DO 250 Ns 1 , N 1 1 

DSHla(<|,*SM(2,NJ-3.*SMn,N).SM(5,N))/ZZ/8(N) 

CH(N)aVIS( 1 ,N)<rOSH1 *XL 

DU1b(«.*U(2,M)-3.*U(1 »M).U(3#N))/ZZ/S(N)*XL1 
CF(N)b2,*VISM #NlA(m.M«U(l,Nl*CUR(Nl/FMDA(I,wn/RE 
250 WRITE(6,«20) Y(N) r S(Nl,US(N),CH(N),CF(N),Pw(M) 

C THE CODE CALCULATES ONLY PRESSURE CRAG AT ANGLE OF ATTACK AS A FUNCTION OF 

c function of local body Radius normal to the body axis, 

CPimaPwU) 

CP2(1)bO'. 

C0P1 C 1 >«0. 

COP2(1)bO. 

DO 3«0 Nal,Ni« 

CPI (N*l)a(Pw(N*n+Pw(N*Ni ) )/2, 

CP2(N+1 )a(PwfN*l ).PW(N*Nl))/2. 

COP1 (N>1 )bCDP! (N)*,s*dt*(CP1 (N1*R(N)*SIN(THE(N))+CP1 (n+I )*r<n+1 )* 
isiN(THE(N*nn 

300 C0P2 ( N ♦ 1 laCDP2(N) + .5*DY*(CP2(N)AR(N)*COS(THF(Nn+CP2(Nfl)*R(N+l)* 

icos(The (N+nn 

C0P(ns2>Pw(|) 

00 350 Ma2,M1 

350 COP(N)a(o,*COS(ALPHAl*COPl(N)-R,aSIN(ALPHA)*COP2(NJ/PI)/R(NJ*»2 
WRlTE(6,«60WR(Nl,CDP(N),Nsl ,U1 > 

460 F0RMAT(2X f 'BODY RADIUS Rb',F12.5, 1 OX. 'PRFSSURF DRAG*',F 12,5) 
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r> oi 


“CONVERGENCE check, COWVFRr.EMCP criterion is that over soo time steps* 

THE HF4TINC, RATES SHOULD NOT CHANRF RV MORE THAN i PERCENT, 

TE1 ( 1 |NU ) **CH C 1 1 Nil )»TE ( 1 »N1 1 T 

tei tiiNi i i*vaps(tei n ini n itei n iNim 

TE2M»M1)aCBIT*rH(ltNin 

IT«Q8SCF(TEl(UNU)#TE?Cl|N1in 

TFOI.EO.Nl nGO TO 261 

UsL/^fOO 

LL»LL*50no 

ir(LL,NE.L)r.O TO 050 
261 CONTINUE 

DO ?60 Na 1 * N 1 1 
WRTTE(fc,0?nN,Y(N) 

llRlTEt6',o55!f7N(M,N5,l)(M,Nl # V(M,N),WPfM,Nl,P(«.Nl,RO(M,N),T(M,N}, 

1M(M t U ) , VIST n ,Hat , M 1 ) 


260 

<100 

405 

420 
4 21 
41! 

/■' as 

450 

551 


CONTINUE 

j c n * i f * K ‘ 1 1 ^ T ^ ^ ^ ^ 

rnait' iT f / / . 1 OV . # ^0 fir T T £ »f & . T * f 1 N ? ss 0 ♦ T 5 ^ 

FORMAT!/, IK, 'OIStIaLONG «nov*, r iX"'3HnrK STANDOFF 0 I ST , * . 5 x * * S HOC * 

« opccc» , 5K f 'STANTON NO , * , 5 X , * S* I N FR I C T I ON * , 5X , * * ALL PR , ■ 
r Tn\, ATfRV.Flft 5»1''''*F10,S*1 n> w !r 10« t ',6X-Fl?.8*NX,F1l«T,bXfF1C,51 
FORMA" (/, NX , *MODV STATION Mn,s'#I3#.Sy* 'OIST, ‘LONG HOD^VtO V) 

fSrmatS/Ux.'ZNSRX^U'.RX.'V^.SX.'HP'.RT.'PSRX.'ROSRX^T'.RK, 

1 »H**8X, 'VIST*) 

FORMAT ( IX, OF 10,5! 

L e L + 1 

IFfL.LF,L M *Xir,n TO 1 


CONTINUE 

RFTUSN 

FNO 
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3URR0UT I K’E EnOY 

DIMENSION PE f?9) ,DPF C29) ,BLDC?9) ,BLEC291 ,UE(29),CUT2( 1 M ) 

rOMMON/F 1/Nl,WU,N|?,N16,NlS,Ml,Mlt # M50,NHl 

COMMON /F t/r> k, i ,r,M2,RM^,r,Mu f RM^,r,M6,r,MY,CM(i #C M9,B2,RE,Tw ( GAMA f SIG.MM 
COMMON /p« 'XY < /Z#PM t RETA,SIRT 

C0MK0N/FS/Sr?9) , SS(291,G(20),CllRf29),MS(29) ,0S(29) ,VS<29) ,VW(?9) 
C0MMnN/F7/'J C 1 0 1 ,29) , V( | 01 ,29) ,TC10 t ,29) ,P( 1 01 ,29) ,RO(101 ,29) , 

1WP(5 01 , f 9) , VTSdOl ,?0),N( 1 oi ,291 ,SH(i 01 ,29) ,PK(29) , VJST( 101 ,29) 
CnMMON/Fn/RT(10t.29),F M OA(inl,?9),ZN(101,29) 

COMMON/F1 0/THFC (1 01 ,29),THFS(101 ,29) ,TV1 noi ,29) ,TV2( 101,29) 
COMMON/FI 1 1 01 ,291 ,nSM( 101,29) ,DwPfiOi ,29) , A6(101 ,29) , 

1 A 7 f 1 o 1 ,29) ,AR(lOt ,29) 

fOMMON/F1h/HTNF,9E1,RE2 

00 3075 Na 1 , N 1 1 
00 3070 Msl5.M11 
IF ( H ( m , n ) ,GF. H INF)G0 TO 3072 
3070 CONTINUE 

3072 RLE (n)sZn(M.n) 

UE(N)atl(M,N) 

PE(N)sP(m,N) 

BLO(N)oO, 

00 3073 J=2,M 
J 1 B J" 1 

3073 B|_ 0 C n ) a8i 0(N)*f2.«(U(Jl ,N)+ll(J,N) )/UE (N) )*(ZN(J,N)*>ZN(J1 ,1) )/2, 
3075 CONTINUE 

On 3085 n a 2 , N 1 5 

3085 OPE (N) a (PE (N4l )-PE (N.i ) ) /XX 
DO 3095 NrN i , N 1 1 , N 1 U 
3095 OPE(N)o(PE(N)-PE(N-l ) )/XX*2, 

OPE Cl )o(PF (2).PE(N1 2) ) /XX 
0PF(N12)sCPFCN16)-PEM))/XX 
00 31 1 0 N 3 i , N 1 1 

CUT2C1 |Mn = r. /Cl .45.5*CZNC1 ,N)M1 )/BLEfN))**6) 
TV2Cl,NjMi)-9E2*«0C1,NjMl)*UECN)*BL0(N)*cUT2n)Mi)/viSCl,Ni M l) 

TV! (1 ,N)aO, 

CUTTsSORTCVISCl ,N)*0UC1 ,N)) 

CUT2C1 |Mi)=VAHStrHICl,N,Mi),cUT 2 Cl |M 1 )) 

00 3115 Ms2 , m l 
CUTsCUTT /SORT (PO(M,N) ) 

UTAUaCIIT/PEl 

PPLIJS»"OPECN)*VISCM,N)/RE/R0(M,N)*#2/uTAU**3 
CUT1*1 .-1 1 ,8*PPLUS 
TECCUTI .1 E.O'.IRO TO 3137 
APLUSs2b’, /90BTCCMT1) 

CNPLUSsZNCM,N)*BOfM,N)*RFl/VISfM,N)*CUT 
ANPLUSaCNPUlS/APLUS 
JP(ANPL'J3.GT’.12.)G0 TO JUt> 
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VI8TL*.a«ZN(M,N)*n ..l./EXPUNPLUS)) 

60 TO 3115 

5116 V!STla,<i*ZN(M,N) 

3115 TV1 (M,maHE*BO(H,Nl*V18TL**2/VIS(M,N)*CtJT2(^) 
5137 mbsm-1 

DO 3130 Mai , mb 
DVIS aTV2(M,\)-TVl ( M ,N) 

IFtDVIS.LE.o'.JGtl TO 3133 
3130 V!$T(M,N}eTvi (M,Mi 

MeM+l 

3133 mbisMJ.m+1 

3110 VI8T(M f NtMBn*TV2(M,N,MBl) 

PF”H«N 



SUBROUTINE SHOCK 
DIMENSION A2(29),VNF(29) 

COMMON /n /N1,Nll,Nl?,Nlft,Nt*5,Ml # Mtt,MSO,Nia 
COMMON/F7/NXM,NXM1,NXM3,NXMU,NXM8,NX“5,NXM(,1 

C0MMON/F5/r.Mi ,GM?, GM J,GMU,r,H5,GM6, r ,M7 t GM8,CM9 f B2,RE#TW,GAMA # SIG,MM 
COMMON/F<j/XX,ZZ,FM,HETA,SIGT 

COMMON/FG/S(?o) ,8Sf?9),Gf?9) , CUR (29 ) , US (29) , OS (29) , VS ( 29) , VN(29) 
COMMON /Fa/ I NT (29) ,TTMBf?P) ,T(-mP(29),TE(29),TE1(29),TE2(29),TE3(29) 
1. TE9 (?9), TF5(29),TF*.(?9),TE7(29) 

COMMON/F7/IJM 01 ,?<n , V(ini .29) ,T<101 ,29),P(101,29) ,RO(101 ,?9), 
JHPM01,?«), Vise 101, 29),Hfl01.?9),SHC 101 ,?9),PW(29),VJSTd01,29) 
COMMON/FB/BT ( I0l,29),EMDA(101,29) ,ZN(1 01 ,29) 

COMMON/Fq/Sl flO l,?9),SSl(101,29),OSt( 101,29), CURl (101,29), 91 (101, 
129) 

C0MM0N/F1 0/THFC( 101, 29), THF8(101,?o),TVl (101 ,29),TV2f 101,29) 
COMKON/Fl 2/T7?( I'M ,291 , Z1 O 01 , ?9 ) , T73 (1 0 1 , 29 ) , T 2(1 (1 0 1 , 29 ) 

A2(l I Ml 1 )sSRM »N1 1 )/(t , + S(l|M11)9CUfi(l|Nt D) 

G(i»Mi)ai, + A?n»Nin*AHfijN»n 
GC1 »M 1 )«VSQHT (G( t |M 1) |G(1 f N t t ) ) 

Tt MP ( 1 1 N 1 l)8Q«Vr.ATHP(P(Mi,i,NXMl),INT(1|Nll)|TEMP(ltNU)) 

TFMPd | N 1 1 )sTFMP() t Mf 1 )/GM1 

TF(1 f N 1 1 ) b f Tf m P( 1 »M l)*r.M7)/(l>TFMP(l,Nil )*GM7) 

VUF(1 » N 1 1 )sr,MB*(TFMP( 1 | N 1 1 ) + GM7) 

VNF (1 |N1 1 )eVS(3HT(VNF f i t N1 1 ) f VNF( 1 jNU ) 1 

USM »N1 1)*VNF(1 »NH)*fVN(l |M1 1)«VS(1 »N1 1 )*A2(1 |N1 1) )/G(l |N11 ) 

TE1 ( 1 (Nil )0(1 ,.t ,/TFU t M 1 n )*VNF(1 |N1 n/G(t »N1 1 ) 

TE2C1 1 M 1 1 )*VS(1 jNj D.TF1 (1 tNH)*A2( 1 ,NH ) 

TES( 1 (Nil )SVN(1 |N1 1)*TE1 f 1 JNI 1) 

TF<I(1 »N1 1 )bT£MP(1 * m 1 1 )/TF( 1|NJ 1) 

TEMPI 1 t N ] 1 ) sVSOPT (TF0C1 > n 1 1 ) j TEMP d j N 1 1 ) ) 

TE5<1 |N1 1)«TF/i(l | Ml 1 )*TEfPd |Ni 1)*GMb/(TEaCl |M11 )AGM5) 
TEA(l|N11)*TF0d»Mll)#GMJ 

TF7d|Nll)*TFhd!NH)*(TE2d»Nll)*TE?d|N!l>*TE$d»Nll)*TE3d I Nil) 
l)/2, 

RO (M 1 , 1 ) NXM j ) bQAVSCATR(TF( 1 |N1 1 ) ,INTd ,N1 1 ) |R0(M1 , 1 ,NXM1 ) ) 
U(Ml,1|NXMl)aORVSCATR(TE2(l tNl 1 ) , TNT ( l , N \ i ) ,u (M l , i , nxm ) ) j 
VI Ml , 1 |NXMl )sO0VSC*TR(TF3d |N11),TNT(1|N11)|V(Ml,t,NXMl)) . 

TCM1 ,1 |NXM1)bOBVSCATR(TE«( 1 |Nl 1 ) , TNT d fM 1 ) |T(Ml , i |NX mj ) ) 

VIS(M1,1 jNXMl )*QAVSCATR(TE5d jNl 1 ) , I NT d j N 1 1) f V I S < m 1 , 1 jNXMl )) 
SH(M1,1 |NXMna(3flVSCAT«(TE6d »N1 1 ) , INTd |N1 1 ) |SH(M1, 1 f NXMl ) ) 

H(Ml,l|NXMl)aQ8VSCATRfTE7d|Nll),INTd|NU)fH{Ml,t ,NXMJ)) 

00 60 Na?,Nl 

HOELP* (3 rO-S<N*Rl «>)/?. 

60 KP(Ml,N)8.GMfl.TF1 (N)*00EL/RT(M1,N) 

(in 65 NbN 12.M11 
DOFL»C3(N-Nia)«S(N))/2. 

65 HP (Ml ,N)»GM(i.TEI (N)#onEL/BT(Mi ,N) 
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wpcMi,nc.u(M, n 

IF(M*. FQ'.2)G0 TO 70 

OS(l I N 1 UbUS( 1 |M1 1 )*5M » >J 1 U 

GO TO T1 

70 08(ifNtll*(f'S(t|Ml ) 708(1 |N1 t)*G(1 | N 1 1 ))/2. 

7i sniNinaSd ; Nin708(i*Nin*TiMP(i 
00 75 7!a1 ,N11 
81 (!»N|Mt)sS(N) 

831 (l|N|Mt )sSS(M 
75 OS! (1 ,N|M1 )sOS(N) 

EM0*(1 , 1 |NXM^al ,7ZN( | , 1 |NXH)*CU»1 d . 1 |NXM) 

BTn*iiMXM).Rin.i f NXM)7ZNn,i,NXM)7THEcn.i.‘JXH) 

RETURN 

ENO 


3U 



70 


suBROiniNe oe»v 

C0HM0N/F1/N1 ,NH,M12,N16,N15 # M1 ,H11,H50,N!« . 

COMMON /F 2 /K Xm, NXMt »NXM3 ,NXmu,^xmR,N)(m^,M*M61 

COMMON/F7/U(iOi,pO),Vf 101 , 29) . T ( 1 0 1 , 29) , P (1 0 1 , 29) , RO ( 1 0 1 . 29 1 
COMMON/Ffl/BT(10t,29),EM0A(!01,29),2fc(i01 # 29) * 

^C0MM 0 N/F9/S1 not ,29),SSl(101,2«n,0S1 riftf,29) # CURlM 01,29) ,Rl ( 10 !, 
COMMON/F 10/THFC(!01,29),THFS(1 01,29). TV1(l0t,29).TV2 (101 291 

tni(2,l |NXM3)a(IJCS,l |NXH3).u(t , 1 1NXM3))/TV1 ( 2, 1 ,NXM3) 

0WP(2, 1 (^ , XN3^*(irP(3. j |NXm3).mP(] , j |NXM3))/Tvj (2, t » N X m 3 ) 

DSH ( 2 , l|NVM3)a(SH(S,t jNXM3)„SH( 1 , 1 jN*MS) )/TVl (2, 1 ,NXM 3 ) 

DO 70 N*t |M 1 1 

Dtl(l,N)»(a.*i|(2,N).3 # *u(l,N)-U(3,N))/TV! ( 1 ,N) 

0MP(l,N)B(«.*WP(2,\)»3 t *wPrl,N)-MP(3,N))/TVin ( M) 

0SH(l»N)s(a,*SH(2,N)»3,*SH(l,N)-SH(3,N))/TVt(l,N) 

1 #*0-3. *11 (Ml ,N)«M(MSO,N) )/TVl (Ml ,N) 
'm ,= " (U ** SH(M1 1 » n, - 3 ^*SH(M1 ,M).SH(M50,N)) /TV1 (Ml , M ) 

OhP(Ml,N)a.fa,*wP(Mt t ,N).3.*«Pf»'l,M)-wPfMS0,N ) ) /TV i f i- 1 # m , 

*8(1, 1 |NXM)b r 1 Zl(l,l| NXH))/S1(1,1|NXM) 

RFTUWN 

end 
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SUBROUTINE BOUND 

COMMON/ri/M ,NU ,NI?,N16,N15,M1 ,M1 1 ,M50,N1(I 
COMHON/F?/NyM,NXMl ,NXM3,N*H<l,NXH8,UXH*>,N*M6t 

COMMON/F5/G*1 ,0 M 2,GM3,GHa,G»‘S.r,M6,GM7,r;MB,G^O f R2,Pe,TW,CAMA,SIG,MH 

C0MM0N/F<s/tNT(29),TIPP(29),TFMP(29),TF.(29),TEt (29) , TE2 ( 29) , TE 3 ( 29 ) 
1 ,TE«(29) # TE5f 29) ,Tf fcf?9) ,TE7(29) 

COMMON/F7/UM01 . 29), V M D1 ,29),T(U)1 ,29),P(J 01 ,29),R0( 101,29), 

1 wp f 1 0 1 ,29) , VIS( 1 01 ,29),H(101 ,29) ,SM( 1 01 ,?9) ,p* (29) , VIST (101 ,29) 
TFHP ( 1 |N( 1 )*03VG»THH(P(?, 1 |N*M1) ,INTf1 jNIt ) »TENP( J |N1 1 ) ) 

P(I #1 |NXMI1bO«VSCATR(TEMP( 1 |N| 1 ),INT(1 t Nt 1 3 |P(I» 1 |N*M|) ) 

PW(1 |N) 1 )«TPHP(1 |NJ 1 ) 

TEHP(1 |M t)3TFMP( ) | Ml 1 )/TW/r,M) 

R0(1,l |NXH1 )sOSVSCATR(TFMP( j |N11),INT(1 > N t 1 ) » RO () , U N*M 1 )) 

RFTURN 

END 
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SUBROUTINE VECt 

DIMENSION OP ( 1 0 1 , 29) , A9 ( 1 0 1 , 29 ) , A 1 0 { 1 0 1 , 29) , A 1 3 ( 1 0 1 » 29) , A A 1 ( 1 0 1 , 
129),AA2M01,?9) 

COMMON /E l/Nl ,N1 t .NI? # N16,N1 5,Mt .Mil ,M50,Nt « 

CnMM(JN/F?/NXM,NXMl,NXM),NXMa,NXM8,NXM5,NXM<,l 

COMMON/F^/GMI ,0M?,r 3 MT( t t;M/j,t;M5,r.M8,r,M7,GMB,GM9,B2,RE ,TW,GAMA # S1G,MM 
CnMV0N/F/j/XX # ZZ.FM,flET4,SIOT 

C0MM0N/F7/UC 1 0t.29),Vfl01,?9),Tf 101,29),P( 101,29), RU< I 01.29), 

1 wP (1 0 1 ,29),VTS(101 , ?9 ) , m ( l n j ,29),SM(ini ,?9),Pw(29),VIST(10l 29) 

CnMM0N/FB/BTfl01,?Ql,EMDA()ni r ?'l),ZNn01 < 29) 

C0MM0N/F9/SI C101,?9),$S1 (lot ,?9) ,0S1 ( 101 ,29) ,CUR1 ( 101 ,29) ,»1 (101 , 
129) 

COMMON /Fin /TMPC M01»?9),THf$(10l,29),TVl(101,29),TV?(101,29) 

COMMON /F - 11/ Dl 1(101, ?O),D3M(10l,29),DNP(t01, 29), 46(101.29), 

lA7(101,29),AP(t01 .29) 

C0MM0N/ri2/TZ?n 01 ,?9) ,Z1 not ,29) ,TZS(1 01 ,29) ,TZO(101 ,29) 

COMMON/F1 «/AU( 101 ,?9,S) , AM( 101 ,29,5) , AN( 1 01 .29,5) , AQ(101 ,29,5) 

AA1 (1 ,1 |NXM)aVISn,1 tNXM) /Rf*VlST( 1 , 1 |HXM) 

Afct 1 # 1 »NXM)SAA1 ( t , 1 |NXM)*0U(1 , 1 |NXM)*4fc(1 , 1 |NXM) 

A7(t,2|NXM<J) sA A1 ( 1 ,?»NXM«)*0 kP( 1 ,2jNXmU) *A 7 ( 1 ,2|NXMfl) 

TV1 (1, 1 |NXM)sSl (1 ,1 |NXM)*tMDA(1 . 1 |NXM)*HO( 1 , 1 | NXM ) 

AU(1,l,t|NXM)sTV1 (l,t|NXM) 

AU(1 , 1 ,2|NXM)bTV 1 (1,1 I NXM) Aim , 1 1 NXM) 

AUM,1,3|NXM)*TV1(1,1|NXM)*V(1,1»NXM) 

AU(1 » 1 ,<*|NXM)bTV 1 ( 1 , t >NX M )*wP(1 , J jNXM) 

AU(1,1,5|NXM)bTV 1 (1, 1 |NXM)*(H(1, 1 ,NXM)-P(1 , 1 |NXM)/R0(1 , l f NXM) ) 
TV2(1, 1 |NXM)bS 1 (t , 1 |NXM)*R0(1 ,1 »NXM)*|J(1 , 1 |NXM) 

AMM » i # 1 |NXM)oTV?(1 , 1 JNXM) 

AM(1 , 1 ,2|NXM)«T V2C 1 , 1 | NXM) *u (1 , 1 JNXM) +s I (1 . 1 (NXM) *P(1 , \ )MXM) 

AM(1 ,1 ,3|NXM)bTV?( 1 , t |NXM)AV(t , 1 fMXMt 
AM(l,l.a»NXM)«TV2n,1»NXM)*MP(l,1 INXM) 
AM(1,1,S|NXM)aTV2(t»l|NXM)*H(l,lfNXM) 
DP(l#2»NXM5)e»(P(l,?|NXM5).Pn ,M2|NXMS))/2, 

DP(1 ,NJ2|NXMs)a(P(1 ,2»NXM5)«P(i ,Ni2jMXM5) ) /2'. 

A9M,2|NXMa)sEMDAn ,2|NXMa)*THFS(1 ,?|NXMfl)/RT(l ,2tNXM4l) 
A10(1#2|NXM«)*S1 f 1 ,2|NXM«)*THFC( 1 ,2|NXH0)/BT (1 ,2|NXM<I) 

AA1 (1 ,2|NXMO)sEM0A ( 1 ,2»NXMU)*Rn( 1 ,?JMXM0)AV(1 ,2»NXM0) 
A13(l»2|NXMa)*Afl(1 ,?|NXMU)*fOS1 ri,?|MXMa)«AU(l,2.1 |NXM9) + SS1 (t.2l 

1NXM0)*AM(1 ,2, 1 |NXM0))-AA1 (1 ,2|NXH«) 
AN(l»2,l|WXMfl)BTZa(1,?|NXM<J)*Al3(l,2jNXMU) 

AOf 1 »2, 1 |NXMfl) JA9C1 ,2»NXMa)*AM(l,2, 1 f NXMa)*At 0( 1 ,2|NXMft) *A A1 ( 1 , 2) 
lNXMfl) + TV1 (1 ,2tNXMu)*wP(1 ,?|NXMU)/BT(1 ,?»NXM«)-TZ3(1 ,2jNXna) « A 13 ( 1 , 


?2»NXMU) 

AA?(1 ,?|NXM0)»AA1 (1 ,2|NXMa)*tt(l ,2f NXMfl)«EMDA( 1 ,2|NXM4)*A6( 1,2* 
1NXMU) 

A13(1.2»NXMU)sA8(t ,2|NXM«)*(5Sin.2|NXM(l)*AU(l,2,?|NXM«)+SSl (1 ,2| 

lNXMa)AAMn,?,2)NXMU)).AA2() ,2)NXM«) _ . .. ... 
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*N(t»2.2|N* M in *T2U(1 ,?JNXM 0 )*A 1 3(1 , 2 j^XM<i) 

AO(i» 2 , 2 iNx , -<o)BAR(i,?tNXMa 5 *AMni?»?jNx*u)'fAio(t,?tNXMU)*AA 2 (i# 2 i 

lNXM«)+TVl(l,?>NXM«)*ll(t#2jNXMU)*^Pfl,2|NXMU)/BT(t # 2»NXMU)+TV2n»2| 

2NXH(l)*ru«t fl ,2>jMXMa5*v/(t ,2|NVMfO-St (I ,?»NXMU)*Ci'Rl fl ,2 jnxmu)*A 6{1 
3,?,NXM4n -Pn,?»NXHU)*St ( 1,21 NX^«)*A9( 1,2 rNXHU)-TZ3 

<Kl,2t N X M «l* 4 l3(l,2l N * Mt ‘J 

AA2U*2tNXrtU>sAA1(1,2»NX M *O*V(l,2»NXMa)*FM0A(l,2|NXMa)«Pn.2|NXM«) 

Al3n»2»NXM4n*AHM»2|NXNOl*(OStM,2|NXMU)*AUn#2,SjNXMU)+SSl(l»2l 

lNXMU)*A*n#?.3»NXW<n)-AA?M,2jN*MO) 

AN(l,2,3tNX M (nsTZU(1,?|MX^U)*Al3fl.?) w XHU) 

AQn*2#3!NX k '«lsA9(l,2|NXMUT*AM(l,2»3!NXNtt)AA10(l,2|NXMa)*AA2(l,2l 

1NXM£1) + TV1 ( 1 ,?»NXMU)*Vfl ,2|MXHU)*wP(l # ?|MXM0)/fiT(l,?|NXM«).TV2(l » 2 J 

2NXM0)*C'IR1 ( 1 , 2 »NX“ 01 aM(i ,2|MXMa)-Pn ,2 »mxmu)*( a 1 on , 2|MXMU)*EM0A c 

3l,2tNXM<j}*Si ( ' ,2»NXM/I)*C0R1 M ,?jNXmu))*TZ3(1 ,2i fJ Xna)*Al3(i , 2 > N x n u ) 

AA?{l,2|MXM«)BAAin,2tWXM<|)AWP(l,2|NX M <O»tW0A(1,2?NXM<t)*A7tl,2> 


A! S(l,2|NXMU)aAA(1 ,2|NXH(|)*(r>S1 <1 ,2|NXMfl)*AU(l ,2,a|NXK4)+SSl (1 *2l 
lNXM«)*A*M»?,<JfNXMa)l.AA 2 (l, 2 |NXMO) 

ANd ,P,a|MX“ 0 lsTZfln ,?,NXM 0 )*A 13(1 ,P)MXMfl) 

AO< 1 # 2 »U|NX*u)sA 9 ri,?»NX k 'U 5 *A*M,P,U:NXMa) + Aio(l, 2 |NXMU)*AA 2 n, 2 | 
iNXMOJ+EMDAd , 2 jNXMU)*Sl (l, 2 |NXMO)/HT (1 ,P|WXMin*{ 0 Pfl, 2 |NXMU)+Hnci, 
2?;wxMa)*wPn,?jNXMu)*(? I *wpn,?jNx^o)*un.2r KJ xMU)*THtS{i,2j^x^'a) + 

3vn,2^xMio*THFcn,?^^an-A7n,?i.Nx^ai*THEcn,2i 
U N X M 0 ) ) • T 7 3 ( \ , 2 »NX*tt)«A 13 n, 2 |NXM<n 

AA2<l,2rNyM<n«AAt<l,?tNXM<O*H(1,2,UXMU)-e»0Af 1.2,NXM«)*{VIS n ,2> 

lNXMU)/SIG/RF*n.+VlST(l,?,WXM«)*SlG/Sir,T) *DSH (1 , 2 j MX MU ) + U ( 1 , 2 j NXMU 

2)*A6( 1 * 2 I N X M ti ) ) 

At5(t f ?|MXHa>*(OS1 ( l >?jNXMil)*AU(l #2#S|^X M *0+S$t O 


JNXMU)*AMn*?,S|NXMAO)»AA2(l,2|NXM«) 

AN(l, 2 , 5 i^X k 'iJlsTZun', 2 t k iXMu)*Ai 3 ( 1 . 2 |NX M U) 
A0(l,2,5|NXM«)aA9(1 ,?|NXMU1*AMM ,2,S|MXMa)*A10tl . 
JNXMU)ATVin,2fNXMa)*HM,2jNXMU)*>'IP(l,2|AIXMU)/9T(1 


PjNX m U)*AA2(1,2> 

,?jNXM«)-TZ 3 ( 1 , 2 » 


2NXM«)*A13(1,2»NXM«) 

RETURN 

ENO 



SUBROUTINE VFC? 
dimension t?ouon 

CflMMON/Fl /N1 ,N1 | ,N12,N16,N15,M1 ,M1 1 ,M50,N10 

COMMON/FVr.Mj , 0^?,r,M < s,RM U>r ,HS.f;M6 # r,M7 > (;M fl ,GM'»,B?,RE»Tw,G4MA f SIG.mm 
COMMON/F ii/XX,ZZ.EM,RETA,SIGT 

C0MM0N/P7/UM ot l ?<n l vMfM,?ot,Tnoi,?‘)),p(ini l ?«),ao(!Oi l ?<n, 
lhP( 101 ,29) , VISMOl ,29),H(ioi ,?9),SH(iot ,29) ,P*(29) , V1ST(101 ,29) 
COMMON /Ffl/BT(1 01 ,29) ,E*r>4<tni ,?P5,ZNftoi ,29) 

COMMON /F 9 /SI r 1 01 ,?9),SSt (101 ,29) ,OSl ( 1 0 1 ,29 ) , CUP 1 ( 1 0 1 , 29 ) , R 1 ( 1 0 l , 

1 ?9 ) 

COMMON /F 1 1 /nil f T 01 ,29),DSH(1 01 ,20) ,r>WP(101 ,29) ,46(101 , 29) , 

1 47(1 01 ,29) , AA( 1 01 ,29) 

CnMM0N/Ft2/T7?( IM ,?9),Z1 f 1 01 ,?9) ,T73(10l ,?9) ,TZ<l(1 01 ,29) 

COMMON /F 1 5/4lin( 101 ,<i) ,4M(I( 1 01 ,2,U) , 4N0( 1 01 #«) , AGO( 1 0 1 , U) 

on ioo joi,? 

ITn! 

IF (I. EG.?, ANO.hm.fO,?) HsN12 
420(1 »mi)eS 1 ( 1 , 1! » M n*RO(l , II » M 1 ) 

4M0( 1 , T , 1 j M neA20(l )M1 )*(2.*l.l(l , II »M1 )*KP{1 , II | M 1 ) ) 
AM0(l,I,2»MnaA20M»Ml)*(Pn,IT>Mt)/«n(l,II|Ml)+2,*un,II|M1)4 

lud , n i mi ) 411(1 , it , iii M i n 

420( 1 I HI )s4po(t|Mi )*(2,*U(1 , IT |Ml) + wp(i ,II|Mi)) 

4M0(l,I,3|Ml)e4 20(1|M1)*VM,IT|Ml) 

IOO AM0(1,I,u|Mi)cA20Mt M U* w M,TIrMn 

A ?0 ( 1 I M ', ) *S 1 ( 1 ,’! t'M 1 ) *P0 (1,1| M 1)*E. M 0* (t,l|Ml) 

AUO( 1 , 1 (Ml )e4?0(1 |M1 ) 

4uo (i,2iMUsA?o(i| M n*u(i,i»Mi) 

A(m(l,Y|M1)aA20(1|Ml)*V(1,l|Mi) 

AIIOM ,<l 1*11*4 20(1 JMJ )*(M(1 , 1 |M1 )-Pf 1 , 1 |M1 )/R0(1 , 1 (Ml ) ) 

A 20(1 |Ml)aFMnA (1 , 1 |Mt )*R0(1 , 1 » M 1 ) 

4*0(1 , 1 |Ml )8A?ft(1 »M1)*V(1 , 1 t M 1 )-AH(i , 1 |M1 )*JDS1 { 1 , 1 |M1 )*4U0(1 , 1 )MJ 
D+SS1 (1 , 1 jMI , 1 , 1 , M t ) ) 

ANn(J,2|Ml)BA?0(1|Mn4(lin,ljMn*V(1,l»M1).A6(i,lTMl)/R0U,!|M1))» 

lAB(l,1»M1)*(OS1(l,1)Mn*4IJO(1,?)Mi) + sSlM,l|Ml)*AMn(l,l,2|Mj)) 

ANO( 1 ,3|Ml)B42fl(1 jMI )#(V( 1 , 1 |M1 )**24P(1 , 1 |M1 )/R0( 1 ,ljMt)).6flM,l» 
1M1 ) * (OS 1 (1 ,1 • Ml ) * AUO ( 1 ,3|MU*5Sl(I,ltMl)*AMO(l,1,3|Ml)) 

AN0(1 ,<I|M1)S420(1 JMI )*( vn , 1 | Ml )*H( 1 , 1 »M1)-VIS( 1 , 1 | Ml )*DSH(1 , 1 JMI ) 
l/PE*M',/SIG*VTST(l , 1 |M 1 )/Sir,T.)/RUn , l|Ml)-U(l,1|Mi)*Ab(l,ljMn 

2 /won , 1 1 hi ) )-ab(i, 1 1 mi > * t ns l (i , i jmi )*m,io( i jmi ) + $si ( i 
3, 1 | Ml )*AMO( 1 , 1 ,0j Ml 1 ) 

a 3 o = c u p i n,n*si n,n 

aoo( 1 , l »mi )eA3o*Rnn , t |Mi)*vn, i |Mi )4ANr)M , i |MD*tz 3( i , l |Mt) 

400(1 ,2)Ml)s?,*4 10*(»0(1, 1 |M1 )*u(l , l,Ml)*V(l , 1 |M1)-A6( 1 , 1 }M1 ) ) 
UAN0(1,2|M1)*TZ3(1,1|M1) 

A0O(l,3fMl)*430*(-P(1,l I Ml ) *B0( 1,1 |Ml)*(V(t,l »M1)**2-U( 1, 1 fM) )**2) 

i)4AN0(i,3i«n*TZ3ri,i,Mt) 

A00Cl,«|Ml)«A30»(PO(l , 1 |M))*V(1 , 1 »M1 )*M(1 , 1 |M1 )*VIS( 1 , 1 |M) )*0SHC1 , 
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11 « M t )/RE*(l./3IC*VlST(i » I »Hn/SIGT)-U(l , 1 |Ml)*A6d» 
UAN0(l»a|Mn*T75d#ll M n 
AWOfl « 1 |MUs«TZ«M , 1 |Ht )*ANO(l« 1 |M1 ) 

ANOd t?1»\ )s«TZ«( i , HMD* AMO ( 1 ,2»H1) 

ANOd .llMlls.TZa ( 1 . i t M d*AMUd ,i|H1 ) 

ANOd #«lMn*.TZU tl# 1 |*n*ANOd »«»M1) 

RETURN 

END 



subroutine VFC? 

COMMON /El /Ml ,Nlt,N12,N16,N15,M1 ,Mn,MSO,NII» 

C0MM0N/F2/NXM < NXMl,NXM3,N)(M0,NX k *8,NXM5,NyM6t 

COMM0N/F3/GM1 ,RM2 f RM3,GMa >f ;M5 > r,M6 < r,M7,GM«,6M9 f B2,RE,TW,GAM4 > SlG,MM 
C0MM0N/F7/UM 01 ,29), VM 01 ,29) ,T (101 ,?9) ,P( 1 01 ,29) ,RO< 1 01 ,29) , 
1WPM01 ,29),VTS(19t ,?9),H(101 ,29),SMfl01 ,29) ,Pw(29) , VIST (101 ,29) 
COMMON/FA/BT M 0 1 ,29) , F “0 A C 1 01 ,?9),ZNf 1 0 t ,29) 

C0MM0N/F9/S1 ( 1 01 ,29) ,ssi ( 101 ,29) ,0S1M 01 ,29) ,CURl (1 01 ,29) ,R1 ( 1 01 , 

129) 

COMMON/F 1 0 /TMEC f 1 01,29 ), TMES( 101, 29 ),TVlf 101, 29 ),TV2( 101,29) 
COMMON/FI 3/4111 (lot ,29,9), AII21 101 ,29,5) 

R0(2,2»nxmb )«AU1 (2,2, 1 iNXMfl ) /S 1 (2 , ? > nxm® ) /Fmoa ( 2, 2 j NXMR ) 

U (2» 2 1 NXMB JaAtM (2,?,?|NXMfl )/AII1 (2,2, 1 jMXMfl ) 

V(2,2|NX*8 ) bAU) f 2, 2, SfNXMfl ) / A 1 1 1 (2,2, 1 I NXMfi ) 
wP(2»2|NX*'fl )bAM1 (2,2,«iNXmb ) / AUK? , 2, 1 t NX“8 ) 

TV1 C2, 2 » NXM0 1 a {')(;>, 2 1 NXMB ) *LI(2, 2 1 NXM A )+v (2, 2 1 NXMR)#V( 2, 2 1 NXM8) )/2, 

SHf 2,2|NXM0)s(AM1 (2,2,S»NXMB)/AtJl (2,2, 1 jNXMR)-TVl (2,2)NXMfl))*GAMA 
Ht2,2»N*MH>=SH(2.2»NX*A)+TV1(2.2lNXM0) 

T(2,2»NXMfl) = SH(2,?|NXMB)/GM3 

P(2,2rNXMB) 3 G“9*R0f?,2|NXNBi*gH(2,2|NXMfn 

TV2 (2, 21 NXMA) aVSGR T (T (2, 2fMXMBi, TV? (2,2| NXMB)) 

VIS(2,2iNXM n)B T (2,?jNXMfl)*TV?(?,2jNXMfl)*G*6/{TC2 # 2|NXM8)*GM5) 

RETURN 

END 


Hi 



SURROUTIMF vfc« 

COMMON /FI /Nl,Nn,Nt?,N1»,Nl5* M l* M U»M5n*Nt<» 

COMMON /F 3 /G Ml ,OM?,r,M^ t GM<j,r,MS,OM6,GM7,G M ft,G M, >,02 f RF,Tw,GAHA,SIG,MM 
C0MM0N/F7/HM 01 ,P91,V(1«1 ,?9),T(101 ,29>,P(101 ,?9),BO(101,?9), 

lWP{101,29>,VtS(101,?9),M(lot,29),SM(101,29),Pi»(29),vlST(lG1,29) 

COMMON/FP/BT (1 01 ,2«>) ,EMOA ( 1 01 ,29) ,7Nfl 01 ,29) 

COM MO N/F 9/31 (!01,?O),S81(101,?9),0S1 ft01.29),CURin01,29).Rl (101, 


129) 

COMMON /F tO/TMfC. (101 ,2°),TM?S(101,?9),TV1 (101,29),TV2(1C1,29) 
COMMON/F13/4U1 (101, ?9, 5), *112(101, 29, *5) 

R0(2.1 |M1 1 la*U1 (2, 1 , 1 |Mi l)/si (?, 1 1 M n )/fMDA(2, 1 I Ml t ) 
U(?.l|Mina4Ml(2,l,?|Min/AlM (?,l,l!Min 
V(2,llMll)s4iJl(2,1,3|Mil)/AUl(2,l»lJ M ll) 

WP(2* t 1 M 1 1) (2# 1 1 ¥ 1 1) 

TV1 (2, 1 )M1 1)9(11(2, 1 I Ml n *!»(?, 1 I Ml U+V(2 , 1 |M11)*V(?, 1 | Mil) )/2. 

8M(2»ll M 11) 3 (*01(2#1,0i M 11)/4U1(2,l,l|Mll)-TVl(2,l|Mli))*G4MA 

H(2. 1 »M1 1 )aSH(2, 1 |M 1 1 )+TVl ( 2 , 1 |M 1 1 ) 

T(2»ll M l 1 ) e SM(2,l l M l 1 ) 7 0*3 

P(2»1) M 1 l)sG M 9*Rn(?,1 | M lt)*S H (2* 1 jMlj) 

TV2(2. 1 |M1 DsVSORT (T(2, 1 |M1 1 ) » TV2(2, 1 |Ml l )) 

VI S (2 # 1 ) M 1 1 )«T (2 , 1 l M 11 )*TV?(2,1 | Ml 1 )*GM«>/(T (2, 1 |M1 1 )+Gm5) 


RETURN 

END 



SURpout i ni- vfcs 

COMMON/FI /Ml ,N11 , MJ 1 ,HS 0 ,W 1 « 

CnMMOM/F?/NXM,MVM1,MXM1,H*Mi4,M)( >, H,MXM5 f K.yM61 

COMMON n VG* 1 ,G“?,r.Mj,r,M«,RMS,r,Mh t r.M7,r,Mfl,G M 9,R2,»F,Tw > GAMA,SlG,MM 
COMMON/f 7/U( 1 01 ,29) ,Vf 1 01 ,?9) ,T f lot ,?9) ,P( j 01 , 29 ) ,RC( 1 01 ,?<}) , 

1 wP( 1 01 # ? 0 ) , vtsn 01 . ?Q) ,M (1 0 1 ,29) ,s«(t 0 1 , ?<n , P*(?9) , VIST ( 1 01 ,29) 

CQMMOM/F R/BT fl 01 , 2 9 ) , F M 0 A ( 101 , ? 9 ) , Z M f 1 0 1 , 2 9 ) 

C0MM0N/FQ/S1 fl 01 ,29) ,SSI f t 0 1 ,29), 081 ( 1 0 1 , 29 ) , CUR 1 ( 1 0 1 , 29 ) , R 1 ( 1 0 1 , 

i?<n 

COMMON/ Ki O/TWEC ( 101 ,?9),THES(i ot ,?9) ,TV1 (1 01 ,29) ,TV2(101,29) 
COmmqm/f 13 / aim (toi,?Q,S), a i!?(ioi, ?Q,S) 

R 0 ( 2 , 2 j NyM^nsAUl ( 2 ,?, 1 jMXMfcn /SI ( 2 , 2 |NXMhi ) /emoa ( 2 , 2 jMXM 6 i ) 
Uf?,2|UXMM)sAU1 (2,?,?»NXMM)/AII) (2, 2# 1 ) NXMM ) 

V(2, ?t NX MM) sAUl (2,2, IfMXMM ) / A 1 1 1 ( 2 , ?, 1 I MXMfcl) 

WPf?,?|NXMfei leAlll (?,?,4>NXMh1) /AU 1 ( 2 , 2, 1 }MXMf,n 

TVl(2,2|NXM(,ns(U(?,?)HXMF,l)*U(2,2jNXMM)*V(2,2lNXM61)*V(2,2fNXM61 

n )/?. 

SH(2,2|NXM6ne(AUl(2,2, t ;jNXMfc1)/AlJlf2,P,l|NXM6l)-Tvn2,2»^XM6l))* 

1 RAM A 

H(2,2|NXM6nsSH(2.?,NXMM)+TV1 ( 2 ,?jNXMhl ) 

T(2,2|N*MM )sSH(?,?,NXM<>n/r,M3 

P(?#2|R* M 61 )«G M 9*Hn(?,PiNXM(,t )*SH(? l p,NXMhl) 

TV?(2,2»^y u M)=VSnRTfTf2,2|NXM»,i)|Tv?(2,?|NXM61)) 

VIS(2,?»uxmm )oT (?, ?t NXMbl)*TV2(2,2|NXM(,i 1 *GM#>/(T(2,2 jNXm^i) + gm? 5 

DO 100 NcNl ,N1 1 ,N| a 

MN1 «N-1 

NN 2 SN -2 

MN1=n»3 

u(p, f «)Mii)3?^s*i»f?,NNi|Mn)-2,*ii(?,NN?|Min+ f s*ur?,NN5fMin 
Vf?.N|MJ1 )*2j,***V(?,*JWt I Ml 1 )-2.*V(2,NN?|Mt 1 )*.5*V(?,NNJ; ><■ n 
P(2,N|M1 ns?'. 5*P(?,NM1 I Mil )-?,*B(2,NN?|M1 1 ) + .S*P(2, * N ( J MJ t ; 

HP (?»N|M 1 t)S2,S*V»P (?,MN1 |M1 1 )«?,*WP(?,MN?jM 1 1 ) + ,S*«P (.?, NM? ;t 1 1) 
SHC2,N|M|l)B? t 5*SM(2,MMl|Mll)«?,*SH(2,NM2|Mll)*,«:*‘5H(?,NNJ.-i«'l 
H(?,N|Ml 1 )sSH( 2,N|M1 1 )*(IIC?,N|Mll )**24V(2,N|Mll )**2)/2, 

T(?,N'|M 1 1 )*SM(2,N|M11 ) /0M5 

B0(2»M|M1 1 )»P(2,NjMl 1 )/r,MO/SH(?,N|Ml 1) 

100 VIS( 2 »N|MH) 8 T( 2 ,N|M 11 )»*) ,5*GM6/(T(2,NfMl 1 )+r,M5) 

RFTUHN 

2 ND 



non on nn 


APPENDIX B 


LISTING FOR ZERO ANGLE-OF-ATTACK CODE 


C 


PROGRAM RPODV c I NPIIT, OUTPUT, TAPFStf INPUT, T APFfeaOUTPUT ) 
TIME-OEPFNOEnT FINITp-OIFFERENCE M E THOO FOR CALCULATING 
AND TURBULENT HYPERSONIC FLON9 a Hot ( T HLL'NT AXISYMmeTRIC 
ZERO ANGLE OF attack. LOCAL time step IS USED TO MARCH 
SOLUTION IN TIME, 

REFFRENCF ATAA PAPES no. 77-17 ? (KUMaR and GRAVES) 
NHODYel FOR S p HERE-CONE , CALL RCON, 

NPODYei 


NPOOYs? 

iriNBooy'.EO.ncALL rcon 
IF(NBODY.EQ.?)Call PHYP 
STOP 
END 



kb 


The Laminar 
BODIES at 
THE 



SUBROUTINE HCON 

DIMENSION Cr)PnS),COF(lS),CDU 5 )»Ow(i 5 ) 

.?^Mo^ I ?Tzrr^ t o^;:^E^ l ?^o?^ f 4?:^Itr^ <l5, ' TEMp^ctol5 • 2€loi ^ 

ss 

COMMOfi/F 1 / N t , N1 U,N5n,NS2 # Nl ,H1I ,MS0 
COMMON/Fp/NXM^yMt ,NXM5,NX M a,NXM5, NXMfl#N)(M 2 

COMMO /F3/r,Hi # r,H? ( (;Mj ( i;M5 (GMj) RH j 5M(( GH(? GAM* Sir MM 

C0MMClN/Fu/XX,7Z,FM,BETA,r)V,0Z * C#T»,GAMA f 3IG»MM 

«C0MHaN / F s / S f, 5).S s nS),G < lS) f c UR(15).USn5).0Sn5).VS(15j,VNnS) 

Es:p:aj}:!::!i!:K?;{i?:;i;. , sir,«i : ; ; r> 5 , 

CH{1)15)»0. 

V1STCI. 1*1515)80. 

constants ANn freestream conditions 

PTs<I,*ATAN(|,) 

GAMAel ,u 
RGASs?87. 

CPeGAMA*RGAS/(GAMA«l ) 

SIG*,72 

SIGTs.R 

ALPHAsO. 

CCe.OOJ 
FMb 10,3 
TFeR6,?6 
PFnlOO. 77 

ROFsPF/RGAS/TF 

THC=«B.aPI/ibo, 

TWa330.3 

TWsTW/TF 

RMs.03175 

OVs.185 

BETAS 1 ^ i 


U5 


I 



noon 


c 

c 

c 

c 


FDT»1. 

tHFSHs2,*Pi/,fl 0 . 

VFsPM*ffiAMA*BGAS*TF)**,5 

VISF*1 . BSR/ ] 0,**6*TF**1 ,5/CTF+l 10.) 

RFbVF*R0F*RN/VI5F 

LMAX IS the maximum number OF TIME STEPS AFTER which THF poftro*- 
LHAXs2oonO EVFN ?H0UGH THE SOLUTIOn m * v NOT BE CONVERGED, 

JI J'iTS 6 *' 0F ME8H P0INTS ‘LONG BOOT SURFACE UPT 0 WHICH THF Flow 
IS IAM]naR, THE FLOW BECOMES FULLY TURBULENT FOR NT*J. L 

N T Be * 

lTURBa2S000 

for laminar flow, call vfc3 *no for turbulent flow r ai i veru 
t™" ’J ,M ' NU -" E » »' ”EB*TIONS UPIO -HKH THE fuS RehmJJ- 

CONVERGENCE CRITERION. 

CRJTa.m 
GM(J)3J , 

GM(?)3-1 . 

GH(3)a|, 

GM(P)aj, 

GM 1 e 1 ,/GAMA /fm**2 
GM?b2,*GAMa/M .♦gama) 

GM3oCP*TF/VF**2 
GM5* 1 ! 0 , / TF 

PMfjBl , + GMS 

GM7s(GAMA-t,)/(GAMAfl.) 

GMflad .♦G4MA5/(? 1 *gama*fm**2) 

OMpa(GAMA-l , 1/GAMi 
GM1 OsTw*GM3/GAMA 
HINFcO,9<35* ('.S + GM3) 

RFlsSQRT(RE) 

RE ?a, 0t6A*R£ 

HFSH SIZFS ANO VECTOR LENGTHS 
Mst5 

N](|SN]>1 
NSOcNl-2 
N52*Nl«i 
MialOt 
■M aMl.j 

"< 4 «• 2 

M51a«l-3 

«52 = mi *<4 
NXMjNS 

N X M 1 3 N '/ M rn M I J 
IJXM..V Ni-1.1 

N»3CNXM«4 
r.i'H , '■> X X .. M * 

»; : ;..aNXM»M j * \ 

'4 < rt f, S K x M „ ^ | „ 2 

NXM*/SNXM»3wM 1 * 1 


1 + 6 ' 



NXM8sNX*«2*m i • l 
NXMPeNXM.u** i + i 
NXM1 03il*NXM 
NXMM3MXM„P*M \ 4 1 
NXMM2aNXM-2*Mt-<| 

NXMM3aNiX IJ -«*M 1 «2 
C STEP SIZES 

DZa-1 ,/hh 

xxc 2 , *nv 

?z=?.*nz 

V(1 jNl)snsVlNTLtO,,DrjVt 1 j Ml ) ) 

Z(ljMS5aP?»VT»iTLri, f nZfZM|«n) 

MsPI/Z.-iC 

X (•’ : • !?Y'-V t ti. 

RHs ( ( v *<!• '■ « ) *S I Ni ( T>-* 5 *c^S ( T MC } I 
C '■ J ' lo 'HE “k* v EOT ' T !i 0 STP£AH £ie JUsTUPF PTINT. 

MJs»l/0T 

NJ3 = "'J 1 >1 

C I X -’. 7 I*'. SHOCK shape AMO SHOCK SLOPE. 

5*1 t"Ji )s t ? r 
PO *7 M«MJ3 ,m| 

17 .i( V' 0 ' k r- 1 1 ♦ 0 V * ? T M ( T h E S h ) 

5 POOY ij E 0 M L 7 P Y , 

or ’o nbi,mi 
I E i V ( *J } , r. T . A 1 )G0 ’0 11 
THFfM= r T/2,»Y(M) 

O’. 0 , N ; =1 , 

R(N)3C0S(THp(NJ) 

GO TO 10 
11 Tup (M ) =THC 
CUR'NJsO. 

R(N)=CnS(THC)4(Y(M)-Al)*SIMfTHC) 

10 CONTINUE / 

DO 20 Nat , N l 
VSCN)3COS(THF(M) ) 

20 VN(H)a.StM(THE CM) 

DO 35 Mat, Ml 
SI f 1 ,M|Hl)oSfM) 

COR ! ( 1 , M j M J )*CUR(M) 

R1(!»N|M1)br(mj 

35 THE 1 ( 1 , MjMl ) s The (M) 

THEC(l,ljNXM)eVCOS(THFtfl,ijMXM);THEC(l#lfNXM)) 

THFSd # 1 t w X M )«VSTN(THEl (J , 1 ,NXM) jTHESC l, 1 |NXM)) 

JEJT (1 jNnsORVINTL ( 1 ,H) j IMT fl »M1 ) ) 

C CALCULATION OE TMF STARTING SOLUTION. 

DO <10 Na 1 , N i 

40 TFMP(N)aSIN(THE CN))**2* (,«>♦. l * (N-i ) /mi ) 

R(1 # 1 >NXM)*Q8VSCATRfTEMP(l | N 1 ) # INT(J,N1 ) |P( 1 , 1 , NXM) 1 
TfHP(l,Ni)ao. 


bl 



V(l,l jNXM)*RHVSC*THtTEHPn»Nt ),lNTd»Nl ) | V ( t f 1 f NXK) ) 

U(1»i jNXM)sGavSCATRfTEMP(ltNl ),INT(1}M 1 | U ( 1 , 1 t NX* ) ) 

TEHPCljMI )arw 

T(l* J JNXM) s 0RVSC*TR(TEMP(1|W1 )#WdfNl )»Td.l|NXM)) 

SS ( t ) so , 

DO 15 NS2.NJ 

15 ssfw)afS(N+n-S(N-n )/xx 

DO 16 N=WJ1,W1 

16 SS(N)s-(y,*Sf^-l)-?.*SfN)-SfM-2))/XX 

TE5C1 (Nijsssm rn > / c i ,fSd |Nn*cuR(i »wt n 
TF6C1 jMllst .♦TESM »fn*TPSn r n 

7F6( 1 tW1)aVSfJPT(TE6n fNI) |TH6(1 jNm 

7F7d J NDa-f vud jNI )-VSd |M )*Tfc5d }Nl ) )/TE6d |M ) 

TE MP( 1 j Nl )a(TF7(l JN1)#TF7( 1 (Nil /G ms.g m 7 ) *G * 1 

P(M1 , l'jUXMn-ORVSCUPfTfcMPd |N1 1,INT(1 |Nn |P(M1 , 1 | NXM 1 ) ) 

TfMPd ;M1)aTFMP(t jMd/GMI 

TE(! m )»(TEMP(l »Ni)fGM7)/d . + TEHP d n n *GM7 ) 

7Md»M)ad’.«l,/TFd|Nni*TF7d>M)/TF6d|Nl) 

TF2d jN1)aV8d jMlWTEt d 1*1 )*TE5d |N1) 

TE 3 C J 1M1 jrVMd IM )+TEl d |N|) 

7FO(l » N t )sTF“Pd J k| 1 ) /TFO INI ) 

U(M1 * 1 J N!X*narj0VSCATWCTT2( 1 ) , IhT(l jNl) » U < M t , 1 |NXM1 ) ) 

Vf''l # lJNX' 4 naf3flVSC*TRtTE3(1 fWUrlUTtl |Nn|VfMl,l|K!XHl)) 
7C*M # 1 J FiXHl)afJ»VSC*TP(TEft(l |N1 )#INTC1 |M1) »T CMJ , 1 fNXMl) ) 

DO 50 Na 1 , N 1 

DXUa(U(Mi ,N).Ud ,N) )/«ll 
DXV=(V(H1 ,Nl«Vd,N))/Mlt 
OXPr(P(Mi,N).P(l,M))/M11 
OXTaCTfMl,N>«T(l,M))/Mt1 

I'd 1 )BOpVlNT(.('UC1 ,N) #f>*U»IJ(l ,N jMl 1 )) 

vd .NjMi nsuRvmLtvd ,M)#nxvivd , NjMn n 

P(1 ,*1*1 1 JSOPVINTLfPf I »N)»riXP|P( 1 ,N JM11 ) ) 

50 t : # n •,« n )30AviNTL(Tci#M>#nxT|Td ,M| mi i ) ) 
r?5d » 1 t MXMJsPf 1 , 1 {NXM) / Td , 1 jNXMl /G M 1 
5> , l iNX^JaTCl , 1 }NXM)*GM3 

Hd,n NJ x M ) B 3 H n,iJ N| xMwnid,i|Nx^)**?+vd,itNxM)**2)/2, 
TV2d 1 1;NXM)-VSI5RT(TC1»1SNXM) »TV2(1,1 jNXM)) 

V18(l.tlMXM)sTd,1>NXM)*TV2d,); N XM)*GM6/(Td, 1|NXM)+GMS) 
pl=(BETA+l,)/fHFTA»l .) 

0? = ?,*RET A/AL r> G (hi ) 

TZd :*m )sri**cz< i »*n ) 

TZ! d dM)3«ETA*(TZd ?Md-.l ,)/{TZf : |M) )*r. ) 

05 30 Nsl.Nl 

3P Z1 fl »N.M1 )3 " l 1 djMI) 

?N(! , 1 jN<V)=< 51 d , i jNXM>*d .-Zt d • \ JNXM) ) 

CMOA d , 1 |NXM)sl .♦7Md , 1 ? NX* )*?: IIP* d , 1 »MXM) 

BT( 1 # 1 »NMK)b«? 1 ( l , 1 |NX«*?+ZNn , 1 |NXM)*TMFC( 1 » 1 »NXP) 

TZ2( 1 , 1 )NXM)bRETA** 2-71 d , 1 ,NXM)**2 

T?3d»lJNXM)e? t *Z1d#lfNXM)/TZ2d<l»^XM) 

7ZU(l»l;NXM)a82/TZ?d.ljNXh) 

WRITE (6,3<»9) ALPHA, RE, f» 


U8 



399 FORMAT (/, IX, iALPHA',Ffl.5,5X, 'RE YN ' NO , F | 5 , 3 ,'5X i 'MACH N0,',F8, J,/) 
WHITE <6,u35)fY (N) , THE fN),P(M),CUR(N) ,Nal ,N1) 

00 60 Nel • Nj 
WRITE( 6 ,« 2 nN # V(H) 

WRITE ( 6 , 30 ) 

WRITE ( 6 , 0 ftO)fZiCM«N),UCM»N) # V(H,N),Pr“,N),RP(H f N),T(M,N),H(M # NJ, 

1VTST ,Hnl , Hn 

60 CONTINUE 

1130 F0RMAT(/,6X,'Zl',9X,'U',9X,'V',«X,'P',8X,'R0',9X f *T',9X,'H',9X, 

I 'VIST') 

<135 FORMAT ( 2X , # Vb'«E 10 ,5# 5X , 'THE*',E 1 0 , 5 , 5X , 'Rs' , F 10 # 5»5X* 'CURs'.F 1 0,5 

n 

0 Hb *1 1 

TFMP2( 1 1“1 noOeviNTL t J .+OM,OMjTEMP?( t I Ml) ) ) 

00 2 1 * 1,11 

2 AUlCl,l f I|NXM)»l, 

TERM »N1)*1 , 

TERM IM )s0, 

TE10M JNI )*GH10 
L° 1 

C MARCHING in TfMf STARTS HERE, 

I CONTINUE 

MM =1 

C LOCAL TIME STFP CALCULATION, 

VAIRM , 1 »NXM?>bGAma*P(1 , 1 JNXM2) /H0( 1 , 1 »nxm?) 

VAIRM • 1jNXM?)bVS0PT(VAI«( 1 , 1 »UXM2) |VATH( 1 , 1 |NXM2)> 

r>T< 1 . 1 |NX u 2)sE0T*(7N(2, 1 ,NXM2)-ZNn , 1 »NXM2) )/(VABStV( 1 , UNXM2) , 

ITHE2C1, 1 »NXM 2 )) + VAIR(t, UNX«2)) 

DO 70 Naj , N i 

70 DT M » N j mu )sOT( 1 , N»mi i)*TFmp2M »H1 1) 

TIMPMjNl )BORVGATHRrOT(MU,llNXHi),TNTM|Nl > > Tl HP ( 1 r N t )) 

0T{Mt , IjNXMi )BORVSCATR(TIMPn ,N1 ),INTM|N| ) | DT ( h i , l ; N XM i ) ) 

C PRFDICTOP STEP. 

S5M ) = 0 , 

00 7 u Nb2,NJ 

7<i SSfN)s(S(N*n-S(N-l >)/XX 

DO 75 NcnJI.M 

75 SS(N)a-(<i,*S(N-n-3,*S(N)«S{N-2))/XX 

call shock 
CALL DERV 
CALL VECl 

oo no i b t , <j 

AU1(2,2,T >NXM6)aAU(2,2,IiNXM6)-DT(2,? f NXM6)*( (AH(? f J, I ,NXM6)»AMC2, 
12, 1 }NXM6))/0 ya( AN(3,2, 1 jMXM6)-aN(?, 2, I |NXM(,))/0Z + AQC2, 2, 1 jNXMb) ) 
AUt C2*Nl, I )M5ft)BAMC2,Nl,I|MB0)-DT(2,Nl » M 50 ) * ( ( ah ( 2 , N l , I}M50)-AM(2, 
1NU1, I |M50) )/0V+(AN(3,N1 , I JM50).AN(2,N1 ,T » «50 )) /O Z ♦ Af) ( 2 , N 1 , I JH50 ) ) 
All 1 C2,NJ1 , 1 >MS0)*AUf ?,MJt , ! |M50)+OT(2,NJl »M50)*( (R,*AMf2,NJ, I |M50) 
1.3.*A M (2,NJl,IjH50).AMf2,NJ2,IlM50))/XX-(AN(3,Njl,ljM5a)-AN(2,Njl, 
2! »H50) ) /0Z-A0(2,NJ1 , I |M50)) 

TAUU.IjNl )bQRVGAThR(AUCM),i,IiNXM 1),INT(1|N1 )»T*U(1,I|N1 )) 
TAM(!,I»Ni )BOflVCATHR(AM(Ml,I,IjNXMl),INTM)Nl )|TAM(1 ,I,ni )) 


lt9 



150 


TANKlflfNl )oQflVr,4TH0(*N(M»,l,lTNXHn,IMTnfNl ) f T AN 1 ( 1 # ! I M 1 ) 
tan2(1,IjN 1 leQBVGATHRC AN(H1 1 , 1 , 1 jMXMl) , t^'Tt 1 |Nl ) J T A*2 ( J , I 1 N I )) 
TAQ(1»1*N1 lafl»V6ATHPfAO(Mt < l f l*NXMn#I*JT(l JM )»TAQ(t,IjM )) 

TAU1 (2. 1 J N50)c TA'J(2, 1 |N50)-T I u p (? j ' * 5 n 1 *( f T A f 3 . 1 jNSO)»TAM(2, I »n50) 

1 Wf)Y*(TANl (2, 1 jn50>-TAN?(2,I|*<S0) ) /DZ + T AO ( 2, I t n50 )) 

TAUl(Nl,n*TAUm,n-TI^<NM)*C(TAM(Nl,n.TAW(NlO,in/DY*{TANUM f 

1 I ) ml AN? ( N\ f I ) > /DZ + T AO (N 1 , 1 ) ) 

TAUl (MJ1 , ncTAU(NJl ,I)4.TIHP(NJ1 )*f (U.*TAM(Nj,n-S.*TAM(NJl,l )-TAM( 
1NJ2 # n ) /XX- (TAM (N.M , l)-TAW?(MJ1 , T 11 /DZ-TAQf NJ1 , I) ) 

110 At! 1 (Ml ,2, I jNXMM) = Q0VSCAT«n A til (2,I|M1U) # I\tn|NlU) f AUl(Ml,2,I| 
1MXMM)) 

CALL VFC? 

AM ( 2*1 |l;Mso)sAiiOf2* I j M 5fi)«r'T f?# 1 | m*>0 1 * ( (AMO(2»2* 1 }^50)»amO( 2# l , 1 

l»M'>o])/nYA(A^n(3,ijM5ni-ANO(?.T: M 5t')l/DZ4»QC(2,l|M50)) 

a mi cm , t , 1 ) = auo(mi ,n-oT(»M ,n*( t amo (Mi,?, i ) -amo (Mi , i,d)/oy+(ano( 

lMi,n-A- 0(^H.ni/n7*A'jn(Mi ,m lv _,, 

Aii1(1,?,1j‘'X UM 3='*'fl v SCATR(TFLn: > 'ltJ5,I , 'JTniMU)|Aum,2,l| f YMM'l) 

AIJl C 1 , 2, ?{^X MM )= f )8VSr ATWf TE<5( 1 |M1 U) , T* T( 1 jM 0) I AIM ( 1 (NXMM! ) 

At’1 Ci,?,5)^XMM)~nMVSCa T -'f T,rc, t t fMti),iNTflf^l^)|&Jl ? 

A’ll ! 1 f 2, A 3 ^.<vv)a ;BVSr.ATK(Tf:iO( 1 J N) a) , IAiT ( 1 )N 1U) | AU1 ( 1 ,2, tt| ^ 

SV(L.GT a !.7'J^H)G0 TO ISO 
CAL 1 . V2C5 
0 ' TO 155 
r ' ■ i. «'ec 5 

7 * N I - T 

A ! » COUMf) 

Cnpr>LCT.-B STFP, 

• > ! 6 o I s 1 » u 

A"?(?,2,Tj N YM5)aAU(?,?,Il^XM5)AAUl(2,2,I|NXM5) 
A | l2(2*l»I) M H) a AU0(?,Ij N, ll) + Aiii(2,l,ljMli) 

MM»2 

SS(D=0, 

DP 174 N=2,NJ 
5S (N ) 3 ( S ( n* 1) “S ( 1 1 /X X 
DO 175 NcNJl.Nl 

S5(t)o ,, (4.*S(AJ«l)"3,*S(W)»Sf^“?)WYX 

CALL SHOCK 

call ofrv 

CALL VEC ) 

AUI C2??! T JNXM5)B0,5*(AU2(2,2,I|NXMS)-PT(2,2)NXM5)#((AM(2,2,I ,nxm5) 

l.AM(2!r,I»NXM5)>/DYXAN(2,2,I,NXM5)-AM(1.2,l|NXM5))/0Z + An(2,2,n 
180 ? AU1 (sIn.M , I fMH )«.5 a(AU2(2,NJ1 , 1 »Min+PT (?,NJ1 »hi 1)*( C4.*AM (2,NJ, 

U,i}l).3>AM(2,NJ1,I,Min-AM(2,NJ2,I,Min)/XX-(AN(2,NJl,l,MlU.AN( 

2 l,Njl,IjMlU)/DZ-AR( 2 ,NJl f l,Minn 
un.2i*n«-un,2i M i> 
call YfC2 

190 Jui(2!l!l|Mn)«.5*(AU7(2,l,I»Min-DT(2,1fMll)*(UMn(2,l,X)Hll).AM0 


150 

155 


160 

170 

1?«; 
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If?,?# 1 1*1 l))/DV*ir*N0(?,I|Mtl)-AN0(1 # I }MilJ)/0Z + *Q0(2‘,IfHin)) 
c FOURTH OOER OAMPING. CC IS the damping COEFFICIENT. 

DO ?00 I si, a 

AU2(3,2,T »NXMM2)sAU(5,2, I jNXMM2)*AU(1 ,2,T |NXMM2).U,*(AUU,2, I JNXMM 

12)*AU(2,2. I JNXMM2) )*6,*AU( 3, 2, n NX MM?) 

TEt ( 1 »N*>n)si)«VGATM«( AUC 1 #2, 1 JNXM7) , iNTtl jN50) fTEi ( 1 »N50) ) 

TE?( t »fc50)»f>AVGATHKC Allf?,?, I |NXf»7l ,INTM |N50) fTE?(l IN50) ) 

TE3( 1 tN50)*OAVGATH»(AU(3.2,I |NX*7), INTO ,NSO) |TE3( 1 1*50)) 

TEO ( 1 »N50>s0*VGATHk(AU(a,2,I |NX“7) , INT (1 jN50) >TEU ( 1 »N50)) 

TE5 1 1 I^SOjaTpan jwSOl-TE 1 ( 1 fTF3n |N501-TF?f 1 f N50) ) 

AU2(?# 2, 1 »nxm7)«08VSCaTR(TfS( 1 jNSO) , I*T ( 1 1*50) | AIJ2(P,2, 1 |NXM7> ) 

TE 1 ( 1 JN50}=QPVGATHW( AUfMSl ,2, I |NXM7 ) , I NT (1 »N«iO) f 7 F 1 ( IjNSO )} 

TF?f 1 fNSMsfjsVGA THRf AiK M SO, ?, I jMXM 7) , INTf 1 1*50) |TF?( 1 ;N5C ) ) 
TE3(l»NSO)8nAVGATMp(AU( M 11,?#I|NXM7),lNTMfN50)tTE' l Ifl|N50)) 

TF A ( 1 f N5P)»0AVGATHRf Atl(**l I |MXM7) , IMT ( 1 | N50) jTtaM |N*»0) ) 
TE5(llNSo)aTElCliM«?ni«TEU(lrN50)-3,*fTI2(l|N50)»TE3niNb01) 

A U 2 ( M 1 1 ,2, I jNXM7)sQ(WSCATR(TE5(l JN50), I*T( 1;N50) | AU2(Mlt I I NX* 7) 

n 

AH? (3, 1 , ! j )* AllOf <5, I f ««>?) ♦AUO f 1, I jH5?)»A,*f AUOfA, I | M52 ) + Al t) { 2 , I f 
?M5?))76.*AUn(5,I|M52) 

AU?(?, 1 , I)sAUO(/l, n-ADOf 1 # T)-3.*<AUO(3, I)-AUOf ?, 1 1 ) 

Alt? f Ml 1 , 1 , I) a Allot MSI , I)»AIIOfMl , I )-3,*f AUOfNSO, I)*AUO( M ll#M) 
an(?,3,Tjnxmm3^bAU1 f?,3,I|NXMM3)-CC*f AUf2#5,I |NXMM3)*AUt2, 1 , I J*XM 
1M3)-a,*(AU(?,«,I>NXMH3)+AU(2,2,I>NX M M3)) + 6,*AUf?,3,I>NXNM3)+AU?(?, 
23, T;NXHM3) ) 

TE 1 ( 1 |N5? JsOPVGATHHf AIM {HI I f NX*R) , I NT ( 1 |N52) »Tf 1 ( 1 }N5?)1 
ANf Ml ,2, I |NXM<J)BGAVSeAT«(TE 1 f 1 | MSP) ,TWT( 1 |N5?) | ANf Ml I |NXMR) ) 

ANf 2, Nlfl, I JMS0)SAU1 f 2,N1U,I |M«iO)-rC*f All (2, *52, I jMSCJ-ALIf ? , N 1 , I |M 

!5<n-3,*( AUC?,.m 50,I |HS0)-AUf?,N1 0, 1 |NCion*AU2f2,Mlo, I > M bO)l 

A N ( 2 , 2, I j M 50 )aAUl t?,2,I »NSn)-CC*f AU(? f 0 , I | *50 ) *GM ( I) • AU ( 2, 2 , 
lljMSO).u,*f All(2,5,I»M50)*AUf2,l,I;M50))+6,*AUt2,2,nM50) + AU2f2,2,I 
2fM5on 

AIM (2,2.1 1NXM6) a AN (2,2,1 ;NXM6) 

IF (I ,EO,?K,n TO 200 

AU1 (2,1,1 |MS0)sAU1 (2,1 ,I|M50)-CC*(2.*AU(2,3, I»M50J-»,*AU(2,2,I |M«>0 

n+6,*AU0(2, 1 jMS 0 WAH 2 ( 2 , 1 , !»M50)> 

200 CONTINUE 

AIM ( 1 , 2, 1 |NXMM)sOBVSCATP(TER(l fNlfl) , INT ( l |N1«) f AU1 (1 ,2# 1 |NXPM) ) 

AU1 ( 1 ,2,2>Nxmm)b 08VSCAT«(TEO( 1 |Mla) , TNT ( 1 »AMU) | AiJi (1 ,2,?f NXMMn 
AU1 (!,2,3|AiXMM)aORVSf.ATR(TEQ(t |Mfl1 , INT(1 »Nia) | aim (1 ,2#3|NXMM)j 

AIM ( l,2,a»NXMM)sQflVSCATR(TF10(l >N1«),INT(1|N1«) JAUI (1 ,2,UjNXMN) ) 

IF(L.GT.LT'JRR)GO TO 231 
CALL VFC3 
GO TO 232 

231 CALL VEC5 

232 CONTINUE 
CALL BOUND 

C FOOT VISCOSITY CALCULATION’. 

IF(L,LT.LTU«8)G0 TO 235 

LLeL/2S 

LL°LL*25 
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235 


C 


250 


340 


350 


355 


460 

465 

C 

C 


261 


IF (LL.NE.LJGO TO 450 
CALI DfRV 
CALI EDDY (NT ) 

CONTINUE 

ll«l /SO 0 

L* cll*500 

lFrLL.NE.LTGO TO 450 

SURF'CI HF STING RATE AMO SKIM FRICTION COEFFICIENT C ALCUL A T I ON.-, , 

WRITE (6, a TO) L 

WRITE! 6* 4051 

XL 1 * •F*? / , ) 

Xi =-2 ,*R?/(RF t a**2-i , ) /SIC/RE 


TE ( 1 I N 1 ) 3 C H f 1|N|) 
on 250 Nsi.Nl 

0SHla(4,*SH(?,N)«3.*SH( 1 ,N)-SH(3#N) VZ2/S(N) 

CM(N)3VIS(1 »n)*OSHI*XL 

DUia(4.*li(2,N)-3.*un ,N)*l»C3.N)) /Z7/S(N)*XLl 
CF(N)r2,*VlSM*N)*(0lM-'l(l . N) *c‘' lW (N) /E*0;. ( 1 , N) ) /RE 
QW(M)aCH(N)*'. 5*ROF*VF**3 

WRITE' (6 , 4 20) v(N),S(N),tlSCN).CH(N),CFCN),Pw(N),OW(N) 

COPUJaO. 

00 340 Na 1 # N t 4 

CDP(MYnsC0PfN)4.5*ny*CPw(N)*R(N)ASTN(THE(N)) + PwtNO.)*P:N*l)aSJN( 

1THF (N+n ) ) 

CDF ( llao. 

ECa-B2/(BETA**2-t ,)/ZZ 
00 350 N a 1 . n 1 a 

C0F(N*1)oC 0F(N)4£C*, 5.*OT*(VIS(1#N)*R(N)*COS(THE(N))*(4 # *U(?.N)«3,* 
lun,N)-U(3,N))/S(N)4V!S(l,M4U*«(N+1)wC0S(THECN4n)*( .*UC2,Ntn- 

i3,*uu ,N+n-u(3,N+n )/st N+n ) 

00 355 Ns2,N1 
CDP(N)b4,*CDP(N)/R(N)*a 2 
C0F(N)r4 # *C0F fN)/R(N)*»?/PE 
CD(N)aCOPfN)+COF(N) 

WR1TE(6,465) 

WRITE (6# 460 )(R CN) ,CDP(N) , CDF C N ) , CD ( N ) ,Nel,Nl ) 


FORMAT nX#4Fll .6) 

FORMAT (/, flX, 'COP', AX.'COF'.flX, 'COM 

CONVERGENCE CHECK. CONVERGENCE CRITERION IS THAT OVER 500 TIME STEPS, 
THF HEATING RATES SHOULD NOT CHANGE RY MORE THAN j PERCENT, 

TEI(l|Nl)aCH(1tNt)«TFn»Ni) 

TE1C1»N|)*VARS(TE1C1 |N1)jTE1(1 ,N1)) 
tep(i»ni)b e»iT*CH(iiNi) 

IIb08SGE(TF1 M |N1)»TE2C1 ?N) )) 


1F(II.E0.N1)G0 TO 261 

LLsL/5000 

LLaLL*5000 

IF(LL,NE‘.L)G0 TO 450 
CONTINUE 
DO 260 No 1 , N 1 
WRITE(6,421)N,Y(N) 
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260 

(100 

005 


020 

02t 

0O0 

031 

050 

551 


WRITE (6, 031) " . ' - - - 

WRITE(6*OO0HZM(M,N),U(WfN),V(M,N),P{H,»J) i R0CM,N) # TfM.N) J H(W,Ml < 

lVIST(H,N),M»t,M) 

continue 

IPfll.il.MlGn T 0 55 1 
FORMATflOX, 'NO. OP ITERATIONS', 15) 

FORMAT ( / , 1 X, 'OIST .ALONG 50n Y • , 5 x , « SHOCK STANDOFF DTST, ',5X, 'SHOCK 
1SPEFI)',5X, 'STANTON No, ',5X, 'SKIN FRICTTON'.SX 'HALL OR.'. 
25X,<’*ALL ^FATING, w/Mgr) 


FOP M ATf5X,Fln,S,10x,cio.5,10X,F10,5,6X,Fl?,3, 

1 5 X , r 1 5 , ' ) 

FORMAT f / ,5X, * w 00 Y STATION N 0 . = ' , I 3 , 5 X , ' n T 0 T . 
En'JNA T ( 1 X, «F1 0.5) ” 


ax,Fi;.-.',ox,Fio.5, 

along eoo v =' f f io,5. 


/: 


FnRMAT(/,6X,»/N» f VX,'U', < ?X,'V',flX,"’*,fly 
1 'VTS t * ) 


RO ' , <?X , » T ' I , f H ’ , e> X , 


L*L ♦ t 

IF (L ,LE ha x T GO TO 1 

CONTINUE 

return 


end 



SUBROUTINE bwyp 

DIMENSION COP(15),CDFMS),C0(15)»OW(lS) 

dimension vnBj.THFMSi.cHf is ), tempi t i5i .TEMpanon »znon , 
iTZMon.izi (iod.tsei noi,i5),r.M(a) 

DIMENSION THF2(101,15),VAIH<1 01 , 15) ,0T(1 01 , 15), AiJ2<101, 1S,U) 
DIMENSION TAU(15,a),TAM(l5,4),TANin5,4),TAN2n5,a),TAQ(15,a), 

IT AU1 (15. a) ,CF ( 151,TE»n5),TE9(15),TElO( 15) 

COMMON /FI /N'l ,N1 U,N50,N5?,M1 ,M1 1 ,M5n 
C0MM0N/F?/N/m,NX« 1 ,rXM3,NXMU,N<NS,NXM8,NXM2 

COMMON /F 3 /G*M ,GM?,f,M3,GMB,RM6 t ".M7.GMfl,G M <),R2,RE,TM,GAMA,Sir, f MM 

COMMON/Fa/XX, 77,FM,RETA,f)Y,OZ 

C0MM0N/F5/S(15),SSn5)#G(15),CUP(15),US(l5),0SC15) # VS(l5),VN(15) 

1 ,Rf 15) 

CnMM0N/F(,/lNT(151,TIMP(15),TF.MPnsi,TEn5),TEl <15),TE2(13),TF?d5 
1,TF4(15) ,TF5f 15),TEh(15),TET(15) 

COMMON /FT/U( i 01 , 15), V( 101 , 15)»Tnoi. 15) ,P( 101 , 151 ,RC’:i 01 , *'i> , 
lVIS(101,15),H(101,15),SMMni.l5),P*«(15),VIST(IP';,>5) 
C0MMON/Ffi/rtT(101 .!51,FMDA(1(11,l5),ZNnf»l,l5),RlC1PI,15> 

COMMON /FP/S 1 (1 01 , 15),SS* Cim , 15) ,031 ( 101 , 15) ,FUPt (t O' , IS) 

COMMON /FI O/TWFC ( 1 0* , 15).THFS(101 , 15), TV! f 1 01 , 15) , ~V2(101 , ) 5) 
COMMON /FI 1 /OIK 10 1 , 15),0SH( 1 01 . 15) , Afc( 1 01 , 15) , ABC 1 01 , 15) 

COMMON/FI ? /T 72 ( 1 01 ,151,7* (101 »15)»TZ3(101 , 15),T/0( 101, 15) 

COmmOn/f 1 J/aiii (1 01 , 1 5, a) 

CQmmon/F l a /ah ( i o i , l5,a),AM(ioi # iS,a),AN(ioi,15,a) f AG!(lol,i5,i? 
COmmon/f i 5/a no ( 1 01 ia) , amo( 1 oi ,2,a) , anoooi ,a) , A'lni 1 01 ,a; 

COMMON/F 1 6 /HIMF.RF1 , RF2 » S I GT 
CH(1 |15)a0, 

VISTtl.t t1515)»0. 

C CONSTANTS and FREESTRFAM CONDITIONS 

PI»a,*ATAN( 1*. ) 

G A M A u 1 .a 
RGAS3287'. 

CP»GAMA*RGAS/(GAMA-1.) 

SIGb.72 

SIGTb.P 

ALPHAbO. 

CCa.001 

FMatO.3 

TFb06,?6 

PF«100.77 

ROFbPF/RGAS/TF 

THC*a5. API/180. 

TWb330.3 
TW*TV»/TF 
RNb, 03175 
DYs.185 
BFTAbI.1 


5 1 * 



o o o o 


fotsi, 

THE8Hb?.*PI/i * o, 

VFaPH*(r,AMA*RGA8*TF)**.5 

VlSFat ,UB«/10,**f>*TF**1 .S/CTF + 110,) 

REaVF*P0F*RN/Vl3F 

C LMAX IS the MAXIMUM number nF TIME STEPS AFTfP KHJCM THE PROGRAM 

C WILL 8T0P EVEN THOUGH THE SOLUTION may NOT BE CONVERGEO, 

LMl YagflOQO 

C NT IS NUMBER OF mesh POINTS ALONG ROoY SURFACE UPTO WHICH THE FLOW 

C IS LAMINAR. THE FLOW BECOMES FULLY TURBULENT FOR NTM, 

NTs? 

LTURBo2500P 

FOR Laminar FLOW, CALL VEC3 AMO FOR TURBULENT FLOW, CALL VEC«, 

LTI.IRB IS THE NUMBER OF ITERATIONS IJPTO WHICH THE FLOW REGAINS 
LAHINaR. 

convergence criterion, 

CRlTs.ot 
GH(1 ) a t , 

CM f 

GM(3)at . 

GM(A)al , 

GHIal , /GAMA/ em#*? 

G M ?*2.*GAMA/n .aGAMa) 

GM3sCP*TF/VF** 2 
G M 5a 1 1 0 , /TF 
GMA.nl ,*GM5 

GM7a(GAMA-t , )/(GAM a*1,) 

GMfl a ( l , +G A m a )/(2.*GAMA*FM*a2) 

GMQa (G.ama-1 , ) /GAMA 
GM10aTw*r,M3/r,AMA 
H|NFaO,9Q5*(.SAGM3) 

REIbSORT(RE) 

RE2*.016«*Rf 
C MESH SIZES AND VECTOR LENGTHS 

N1 a 1 5 
NlflaNl.l 
N50SN1-2 

NB2SN1-3 
MtrlOt 
Ml laMlpi 
MSOaM 1 *2 

M5laMl-J 
MB2aM 1 mlt 
NXMaN 1 *M 1 
NXM18NXM.H1 1 
NXM2BNXM.1 
NXM3BNXM-2 
NXM<ieNXM«M 1 
NXMSBNXM.HJ.i 
NXMfceNXM.?*Mt.2 
NXM7«NXM.3*mi*1 

* 


i 
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NXMBsNXm.2*mi -1 

NXMQbNXM.Ubmi 
NXMlOoiloNXM 
NXMMCNXM»2* M 1 ♦! 

NXmm2oK'Xm«2*m i •(( 

NXMM3snxm.4i*mi«2 
C STFP stzfs 

OZa-i , /mi i 

xxe2.*cr 

ZZa?.*OZ 

Y(1 jNnsnPVIMKO, ,DY|YM )Ntn 
Z(t ;Mi )sfj«viNTL(i..r>z»Z(trMn) 

C INITIAL SHOCK Shape 4M0 SHOCK SLOPE. 

S ( 1 ) s , 1 7 
DO 17 Ni2,M 

17 S(N)sS(Ki-lunY*SI M f*H5SH) 

C BODY OOPmptc'y , 

cuRttisi; 

7HE(n5PI/2, 

Ao) ,/TAfJ{THf> *•>? 
naS*lH T f A 1 
00 1C Ve?, Ml 

01sDY*SnRT(rR**2 + R(M«n**2)/(M**?+(V i + fA/B)»*2)*n(M«n**2)) 
02cDY*SORT ( (P**?+(H(N«1 ) +0 1 /?.,)**?)/( R *#2+ (1 .♦( A /R) ** 2 ) * 

i(p(M-n+oi/p'.i**?n 

OSaOY*SQBT ( (n**2+(Wf N-l) +0 2/2, )**?.)/ (B **2+ ( l'« + ( A /B ) *.*2 ) 

J*fR(U-1 )*0?/?, 1**21) 

0«sDY*SQPT( (H**2*(R(N.l) *03 )**2) /<«**?+ Cl , + ( A/B ) **2) * (R (N» t ) +03 1 

1**?)) 

P (n)bR(m* 1) 4.(0 1+2, *02*2, *03 + 0'n/6. 

XhA*S0RT(R**?+R(n)**2)/R 

CURfN)sA»*u*P,**« / ( a + •/(*» (M)**? + B**U*x**2)**r,5 

XbB*X/(A*S0RT(X**2-a**2)) 

10 THE(N)=ATAM(x} 

DO 20 No ) , mi 
V8(N)eC0S(THFfNn 
20 VN(N)e-STN(THE(Nn 
DO 35 Mai, Mi 
8 1 <1 # N I M 1 ) a S f M ) 

CUR! ( 1,N;M1)bCUR(N) 

R1 Cl» N jH| )sR(N) 

35 THE l (l,MfMl)sTHEfM) 

THEC (1,1 1 NX M) evens (THE1 fl , 1 »NXM) JTHEC (1 , 1 |NXM) ) 

THF$ ( 1 , 1 |NXM)aVSINtTHEl (1#1 »NXM) jTHESd* 1 fNXM) J 
INTO * N 1 )o(3RVINTL Cl, Ml jIMTnjMl) J 
C CALCULATION OF THE STARTING solution’. 

DO ao Mb 1 , M 1 

00 TFHPCN)aSIN(THE(N))**2*(.<»+,l*(N-n/Nn 

PC1,1)NXM)»0RVSCATR(TEMPC1,N1 ),INT(1|N1 ) t P ( 1 , 1 | NXH ) ) 

TFMPU »Nn»o'. 
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V(1 , 1 »NXM) cQAVSC*T«< Tf Mf>(l ,N1 ) * I NT ( 1| N t ) I V ( 1 , 1 J N XM ) ) 

U(1 » 1 »NXM)aQPVSr AT»(TEMP( 1 |W1 ),INT(1»N1 ) | U ( 1 » 1 t NXH) ) 
TpMPCl'»«l)*TM 

TM,ljK’XM>rnpV8CATR(TEMP(1tNl ) f lNT(l|Nl )|T(1,»|NXM)) 

ssn )*f>. 

DO 15 M = ?,NU 

15 S 3 { W ) a f [ f M + 11 • S f '» • 111 / X X 

S3(f ; * )=-(ii,*3(N1 fl1-1,*SiN1)*SD r 'in ' /XX 
7£S£ * (Ml *=SS t 1 »M11 /k 1 ,+S (1 jMl *C‘JB ( 1 f N1 ) ) 

T F t> < ; * N ‘ ) - \ . * T F 5 ( 1 I M 1 * T y <5 ( i ; m 1 ) 

?E6( i ;»i n svsntJT ( re sntMiiTFhdrun 

TE7(ll M 1lS"fVM(ljM)«VS(l|Hl)*K5(ij v Ml)/7fct>(l»Ml) 

TfMP (1 ,ni)s(TFT( 1 : M ) *TE 7 ( 1 j M \ ) 1 *G* 1 

PCMl,lt^X u nf!OHWS'ATB(TF.MP(l{Nn,INTM>MnjP( M t|1|MX M 1)) 

TfMr;i i m n = tp MP (1 j '■ ' ) /G* t 

7E( 1 » Mils ( TpMPft |NU*RH71 / ( 1 . + Tf;MPn »N1)*GM7) 

TF 1 n |M):( 1 ' -1 ,/TF f 1 jMl ))*Tfe .7 (1 ;»'1) /TFDfl |M ) 

TF?M l N 1 ) = VS f 1 T Mil - TF 1 f t : M 1 ) *TE5 Cl | N 1 ) 

T F 5 n | N n a v M ( 1 J M n ♦ T K 1 (1 • M n 

T F 0 ( 1 lM 1 sT£mP(1jVJ ) / T p ( 1 J N 1 1 

U(M1 1 1 j WXMlIeHBVSr AT»fTF2(1 f Ml) » INT(1 {Ml) |U( wi , 1 |NXMJ ) ) 

V(V1 , I I My-^n sOftVSCAT«(TF!3( 1 |M1) , IM n JN11 | V(M 1 , ! I vx«n ) 
i(Mi,iiMy u ncosvsrATPrTEy(i»MUfi s JT(itNn»T(Mi,i|Nx^n) 

DO 50 Nsl.Nl 

nxiisCHMt ,N).un,Nn/Mn 
DXVsCVfHl ,N1.V{ 1 ,N))/H1 1 
DXPe ( P ( u 1 ,N).P(1 1 

OXT=(T fMI ,N1.TM .N11/M1 1 

UM ,N jht 1 )=0AVINTL (')( 1 , N) ,nxiljii( ) .Nt'Mi 1 1 1 

VU f N|M1 1 )sOpVlf:TL (V(1 , Ml ,DXV j V( 1 ,N jMl 1 ) 1 
Pf 1 . N jMt nanAVINTI. (Pf 1 ,N1 , OXPf P( 1 . NjmH ) ) 

50 T f 1 ,N jMJ 11sf)SV!NTL(Tn ,M) ,DXT,T( 1 ,H»Mt 1) 1 
ROM , 1 »MXM)bP( 1 , 1 f NXM)/T Cl , 1 J NXH) /GM1 
SHC1 * 1 »NX*)sT( 1 , t |NXM)*GMJ 

Hf1,l»NXN1aSMn,1rNX v )*('lf1,ljNXMW*? + V(t,l|NX k ')**2)/2, 

TV?(1 . 1 >wxm1sVSQRT(T( $ , 1 jNXM) |TVi»n , 1 jfJXM) ) 

V I S < 1 #1 >NX“lsT(l, 1 jNXM)*TV 2(1 , 1 ;NXM)*r,H6/(T (1 , 1 »NXM)+GM5) 
fi 1 a ( BE T A > I .1/(HETA-1.1 
B? = ?.*DFTA/&LOr,(R1 ) 

T2CI »"j ) S «1**fZ(l i«l>1 

TZi(i)MnenETA*(TZ(ij M n-r.)/(T?!(irMi) + i.) 

DO 30 Ne 1 , N 1 

30 Z1 Cl »N)M1 1=TZ1 Cl »W1) 

znci * i i N x M )=s i ( i . i » Nx^)* n ,-zi n , i |Nx M i ) 

EMDA(1,1jNXM)c1,*ZN(1,1 jNXM)*CIIR1 (1 # 1 jNXM) 

BTC1 » 1 |NXM)SB1 Cl , 1 |NXM)tZNf 1, 1 |NXM)*THECC1» 1 |NXM) 

TZ2( 1 # 1 |NXM)aRETA**?«Zl (1 , 1 ,MXm)**? 

773(1*1 jNXM)s 2,*Z1 (l.i |MXMi/TZ?(l»l jNXM) 

TZ«(1 * l)NXM)BP2/TZ2n , 1 jNXM) 

WRITE (6.3RR) alpha, RF.Fm , 

3RO PORMATC/, IX, 'ALPHA', F«.S,5X,'REVN MO . ' , FI 5, 3 , 5X , 'M AC H N0,',F8,3./) 
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W(mE{6#43Sm<N),TMFC*0»R(N)«CUH(N),Ns1 ,N1) 

00 60 Nst.Nl 
WRTTF (6, U211 N, Y (N) 

N&jTF(6*uaOI(Zl( M ,N),U(M,N),V(M,NJ,P(M,Ni,RO(M,Ni,T(M,N},H(M,N), 

1VIST (*•* N 1 .Met ,Mt I 

60 continue 

430 FORMAT f /,6X , 'Zl # »PX,*U # ,PX,'V*,aX,*P',RX,*R0 # ,9X, # T'»9X,*M*,oX, 

1 'VIST' 1 

4,35 FORMAT (2X, »Y= »,F 1 0.5, 5X,'ThFb'.E 1 0.5. 5X,'R=',FI0.S,5X,'CURs',Fl 0,5 

n 

OMa-1 ,/4,/MI 1 

TFMP2 ( t , m 1 t)a(3RVINTL(t. + OM,OM,TE M P2(l 
DO 2 Isl ,4 

2 AU1 ( 1 » 1 » T f NXm ) a 1 , 

TERM |Nn«l, 

TFRUlNDaO. 

TF 1 0 (1 , N 1 T =G M 1 O’ 

Lei 

C MAPCHIWG IN TIME STARTS HERE. 

J continue 

MMel 

C LOCAL T T M E STEP C Al.riJl ATTON, 

VAIR ( 1 , t l NXMJ»)*G AmA*P M , 1 | NXM2) /TTO ( 1 , 1 1 n*m?) 

VATfid , 1 jNX-?)avS'-}»T(VAI«(1 , 1 , N X M? ) , v A 1 R (1 , 1 , NX M?) ) 

OT ( 1 » 1 jNyH2)aF0T#(7N(?, | jNxh2T«ZN(1 , J ,NXM?) )/(VARS (V(! , 1 ,nxm 2) j 

1THF2(1,1|NXm2))4VAI»(1,1jNxm2)) 

DO TO Nal,Nl 

70 DT ( 1 , N , m, U*r>T(1 ,N,Mi t )*TFMP2(, ,*11) 

T I HP ( 1 , N ] IsORVGATHR ( nT {M,1 , 1 |NXH, ) , INT n I M t )|TIMP(1,M n 
DTfMl,t,NXMl)sQBVSCATR(TIMPn,N! ),INT(1,N1 ) »0T(M1 , 1 ,NXM1 ) ) 

c PREDICTOR STFP, 

SSM )S0. 
on 7« n s ? , n j a 

7a ss(N)«s(SfNtn-8(N-i n/xx 

ss(Nns-{a.*s(Nia)-5.*3(Ni).s(N50))/xx 
CALL SHOCK 
CALL OERV 

call veci 

AUl(2!2!l»NXH6)eAlj(?,2,IjNXM6T-OT(2,?,NXM6)*((AH(2,3,I,NX«h).AM(2 > 

!2,IlNXM6n/OV*(AN(3,?,I ,NXM6)- AN (2,2,1 ,MXMfen/OZf A (,c 2,2,1 1NXMM ) 

AU1 C2,N1 ,I,M50)aAU(2,Nl ,I|M50)«OTfP,Ml *M50)*(CAH(2,NI,I|M50)-AM(2 # 

|N1fl,Ij M 50))/0T'»(AN(3,Nl,T,M50)-AN(2,Nl,TjMS0))/0Z+AO(2,Nl,I,H50)) 

TAUn.IfNl TaRRVGATHR( A'J(M, , 1 , T jNXMIT , INT ( I ,Nl ),TAU(1,I,N1 1) 

TAM(1 , I ,N1 )b 08VGATHRCAM(H, , 1 , I jNXMl) , JNTM »N1 ),TAH(t,I|Nl )) 
TANl(l,IjNt I sQAVOATHR ( AN (Ml , 1 , I ; NXM,1 . INT ( 1 ,N1 I, TAN 1(1, I, Nt 1) 
TAN2(1 , 7 J N l TsQRVGATHR(AN(M 1 1, 1 , I • NXM 1 ) ,INT(t,N! ) | TAN2 ( 1 , I , N l )) 
TAO( 1 , I l Wt )*08vr,ATHR( , 1 , T SNXMII , INT (1 |N1 ) , T AQ ( 1 , I f'N 1 )) 

7AU! (2, I jN50)eTAM(2, I |NSn)«TIMP(2:N5ni*( ( TAM( 5, I |NS0)*TAM(2, I ?N50) 

n/DYT(TANt(?,I»N50)-TAN2(2,I,N50) J /OZ+T AO ( 2, I » N50 ) ) . 
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TAUn^l.I)sT4lJ(Hl,n-TI«PfMn*(CT4Hfpgl,I)-TAM(NlU,li)/'OY+ffANl(M* 
i n-TAf:?f ■ i , n )/n? + T4DfMi , t ) ) 

11 0 All! I rMXMM)B0«VSCATR(7AUl (3 # I|N1«), INTO |N1<I) « AU1 CMl. 2,1 I 

1 n x m y } i 

CALL Ml C? 

00 ISO Til ,U 

AUl(2*1»t> f '50)sAi|0(2»I|H50)»DT(?,1jMi;p)*((A M Of? # 2»I|H501«AMO(2#l#I 

l? M Son/OY*(AMnn,lT*-sn)«ANn(?,T;^5on/02*AOO(?,l»Mso)) 

150 A |i l(MJ > !,T) = 4l!0(vi,T)«OT(Mi,i>*f(4Mn(i<j f p f IT»AMO( M l,l,I))/0Y+l'4N0f 

1 M 1 , 1) - a K 0 ( » H , I) ) / 0 Z ♦ Aon ( “ 1 , n ) 

AIM ( 1 , ? , 1 J M(«K ) = 0MVSr ATP r TfP f 1 t N 1 U ) , TMT f t » N1 a) I AUI ( 1 , 1 I a. Y«M ) T 

Al'1 ( 1 >?,?» NY«M)sfJBVSCATPf TF«j(t » M1U1 . TNT f 1 j N 1 ii ) | Ain ( 1 , 2 , 2 > ) ) 

AUlM.?.TiM‘" J )=08VSCATR(TPQMi’Mil),T'jT(l|K'liJ)|A!ll(1»2.3rM''*)) 

Alll (1 , ?,/|^XM!4j-o«VSCAT»(Tp | Of l jNl «l) ,INTn JMM) | AU1 (1 ) 

fL.GT .LTUPp )G0 TO 1 Sil 
CALL V f C 3 
GO TO j^ss 
1 5tl CALL ”5 
155 C 0 m T I *• U E 

CALL ROKv.0 

C COR it; C TOR STf-O, 

00 160 I c 1 , u 

AH? (i t P , T j vx *s ' sa*. ( t , a, T j fi*; i- .*> v <• a jt r ? , 1 1 kxh" ) 

160 AU?f2. 1 ,T» M 1 1 ) = !..' r? f T ; 1 ' • , T ? ‘ M 

MMe? 

SSM 1=0, 

DO 17 0 6 a?, MO 

17« 5S ( N ) = ( S < M 11 "3 ' i ) f ' 

SS(M )e - f a, *s (N1 <!)•"? , *5 : N n •5 'Wi'j' ) ! y X 
CALL Shock 
call DfPY 
CALL vrci 

DO 180 Isl ,a 

180 AtJl{2.?,I|NXH5)«n,«?*T4ij? ??,?,:•' ♦ ('••( ,'.t 

1-AM2, 1 , T »NXM«;n/nVM AK!(2f2. i *6! : • ' }" ,; U/ '■ 

2WX»5 ) ) ) ' " ' ' ' 

tlfl # «? * H 1 
CALL VEC? 

DO 190 Iol.ii 

190 AIM f?, 1 , 1 |M ] 1 )a ,S* ( All? f ?, 1 , 1 1 u i M..rr f?, 1 1*1 1 * « 1 • >0f2, .1,' ni-AMQ 

1 (2. 2 > l | MU 1 ) /OY+C AMtlf ?., I JM1 1 )„460M , » | *\ n “ >4 1 ' P, I f ■ , ) ) 

C FOURTH ODER HARPING, CC IS THE Oahptag COE ; F I ! ' I *' \ '* 

DO 20 0 lei , li 

AD2C3.2, I >NXHM2)a4ll{5» 2* I ?NXMM?)44ii( , T , ( A y j 

l?) + AU(?,2.IjkjyMM21)*6.*4iJ(3,2.TjUXMHp) 

TE 1 ( 1 » Ai50)ro«VGATHR f Aii( i ,2, j |Nxm7) , imt fl t N50i | TF 1 ( 1 ; M50) ) 

TE2C1 »M50 1af38VGATHR ( 4UL?,2 P I f MXM7) , INT n j^iSO) ? T E P(1 »N5C)) 

TE3C 1 >M50)Bt38VGATHRf AU(3.2. I fNXM7) , TMT( 1 |M50) f TF 3 7 1 »H50 n 
TE«(t| K t5O)Ef}aVGATMPtAij(ii,2»I|NXH7),lNTM»A.50)jTF(;f1iMS0)) 

TE5C1 JN50)«TF«(1 JN501-TE1 (1 »N501-3.*(TF3(1 »^50l-r r ?fl »N50n 
AU2(2#2#T |NXM7)®Q#VSCATR(TE i 50 jAI5o)»TNT(i f Ki S 0 1 |Au 2(2.2# I r x h* 7 5 ) 



TFt niM50)»n»»VGATHRr AU(H«5t ,2,I» N X*^?).l N Tn JNSO) iTf 1 (1 »N50)) 
TE2(l|N5(O*G fl VGATMR(AU( M 50,2#IjNXM7),INT(1|N50)jTF2n|N50)) 

TE3 ( 1 »N5o )sO«VGA THR ( All ( *11. 2. I »NXM7},1NTM |N50) fTETSfllNSOJJ 
TF0C1 |N50l=0PVGATHR(AUf-M *?« I J NXM7 > , TMM »NSO) |Tfc t>( 1 »N50n 

Tri^i*NS0l«TF1 (1|M50)-TF<I (li*50 J-3.*fTF2*1 l«50)-TE 5(1 »N501> 

Iu2(MU»?.I jN<*7JaQBVSCATR(TE5<l JN501.INTC1 »NS«) t AU2 (Mi l , 2, 1 »NXH7) 

^ AU? ( 3 # 1 » T l “5? I a AtJO (5 » I | *52) ♦ AUO ( t , I }*5?)»U«*(AU0(ft, I |*52)+AU0(2, 1 1 
2M5?) )+#>.*AUO( 3# I 1*5?) , 

AU2(2, t , I ) “AOO ( '<, I )«.AU0(1 , J)«3.* < *110(1, 1 )»AUO(2. T ) ) .. , , 

AU2 ( M 1 1 » 1 , I1nAUO(Mt )1 , n .4UnfHt,n-J,*fAliOfM« ) 0,I ) -A'lO(Mll,n) 

AM (2.1. 1 |NyMM3isAlll (?,1,I*NXM*3)«CC*( AH(2.S.T|NXNM3 )faU( 2, 1 #I* XlXM 
< I,j(p # a # I s mxmmhi + & U(2.2. I S^XMM3i 2,3. 1 *NXMM3)+A U 2(2. 

2 TM (1 |M 5 ?)innv 6 ATM«(AUl (Ht ,2,1 |NXM0),INT(1 1*52) JTF 1 (1 1*52)) 

AM(M1 ?# IiNXMP) S ORVSCATRfTFt(l»‘>52i,I^n,M5?),AM(«1.2.I|^XM9)) 

iM (2, M L , I s MS ft ) s 4U 1 (?# M u f 1 1 - cc* ( AU Ci? # KS 2 , 1 J KSO ) *AU C2, 1 # I I 

1 J ,MS0)-<J .* ( All (2.3. f l M 50)+ Alt (2, \ , T ! M 50) ) + 8 ,* All (2, 2, I | M50 ) ♦ AU2 < t # 2# I 

2 j M .50 ) ) 

AiM(2.2»Il*X k, 8) s AN(?»2»Il ***«>) 

ALII (2,1 ) n*MM = Aim?, 1 , I » W «G 0 1 -rc * ( 2.* AU (2. 3,! I M 50)«8, *AU ( 2,2,1 S M 50 
1 )+h t *AUO(2#T ? M 5^WA n 2f^# 1 1 T 

200 C-jMTIM t ■'jCATRATFAd |N1«) , JMT (1 |N1 R ) I At 1 1 ( 1 , 2. 1 |NXMM) ) 

; ,r *T'( ( TFO ( 1 I Hi a ) , 1 1- T ( 1 |Mt <1 ) I AMI ( 1 « ?., P» ) 

(’AT- ' T T Q f 1 | M t 'O > I N T ( 1 | M 1 <1 ) I A U ' (1»2#3( N WM )) 

f T'- i r n M j) , IMT(t JM1U) I A’ 1 ' J ,2,u f XHM) ) 


AIJ1 (1,2,1 |NYMMT3fJ*V( 


: r; A v r 

; a v C 

: '• v ; 

T 0 


AT- 


LA T I )*. 
2 35 


A M (1,2.?1 NVmh) 

1 ' 1 ( 1 , 2 , IjN* 1 -'*) 

A'JI (1 ,?,ti 

r ; :> ."t.l’i 
r.’.LL Vf.;l 
rn vo ?:? 

211 .. 'LL V 1 cs 

? 7 ? r n u 7 ! *■■ ' ' E 

r l| H?’ 'MO 

* F~nv v i sens: v v C AtC» 

IFll.LT.LTVRfnGO 

LL=L/25 
U = U*?5 

!F(LL.nE.L)GO TO <J50 
CALL OCR V 
CALL EOOV(NT) 

235 COMTINIIF 
l LsL/500 

LLnLL*S0n 

IF (LL.NF.l )RO TO a 5« 

c SURFACE HEATING RATE AMO SKIN FRICTION 

WRlTECb.aOO) L 
WRITE(6,U05) 

YLia •P2/(BETA**2-1,). .. 


COEFFICIENT CALCULATIONS, 


6o 



250 

300 


350 


a 6" 1 
0 65 


C 

c 
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XLe* 2 ,*H 2 /(BPTA** 2 ^l , l/SIG/RE “ “ 

TE ( 1 J N i )sC«( 1 >ni ) 

On 250 Nb 1 .ni 

DSHlaf<!,*SH(;>,N).^ i *gn(l,N),{!Mf 3 # N}’>/ZZ/S(N) 

CM f N ) r V ?f) { 1 , Ml *n$H • * tl. 

P^l = C a . * ' N >3 . *u M , N) -U(3, N) ) /ZZ/S (N)#XL 1 
" T (N)=^,*Vl3M»N)«( r MMr'UM. , N)*CHRfN)/EHDAn»N))/RE 
ow i *■ ) = ;N)^.s*RnP*vF**3 

WRITE f 6,4120 ) Y{N),$fNJ,liS(N),CH(N) ( CFtN),o*(N),C*fN) 

COP . Hr o . 


THAT 

Than 


OVER 500 T I M £ STEPS. 
1 PERCENT, 


260 

aoo 

005 


020 


PO JOS Nal.NIU 

ithe(N’ i !“j r)P(N) + * 5 * 0V * (PW(tl} * RfN) * 8 i f ^ TH E(N) ) +p.''fN+i :*R(N+n*siN( 
CDF ( 1)S0. 

Frc-B2/';RETA**2-1 , )/ZZ 
DO 350 Ns 1 . N 1 a 

Cnr(NA()sCDFfN)+FC*.S*OV*(VISn,N)*R(N)*rOSfTHECN))*(« *!lf?,Mn3 * 

lUf 1 3.NJ )/S( N)+VTS t 1 ,N+n»R(N*n*C p S(THE(M +J , )*(«[*,) f p N+ n ’ 

1 3, -Ml ( 1 , N+T ) - II f 3, M >1 ) ) / 5 r > J + 1 J ) .. .» 

DO 355 N s 2 » n i 

CDP f n ) = a ', *c DP ( n > /» ( u ) **2 

CDF (N ; 30 . *CPF (N ) /a ( N) It ;:/Rf 
CDf N)2CCPfNl4.Cf)F(NT 
WRTTE (b, 065) 

WOT TE (6,060 ) fR(N) ,CDP(N), cor fNJ.cn f'O.NxI ,N 1 ) 

pnrj.«AT(ix,0Fn ,6) 

FORMAT (/, A X, »R*,AX, *CDP*. AX, 'COE', AX.'CC*} 
convergence CHECK, CONVERGENCE criterion is 
THE HEATING RaTFS SHOULD NOT CHANGE hi PORE 
TE 1 (I |N 1 ) =CH( | |M ).tf f 1 |N1J 
TEl M |N1)SVARS(TE1 (1 ,N1) |TF1 (1 i'nijj 

T F ? ( 1 i N 1 ) a CRTT*CH(1 ,N1) 

IIcOfiSf.E (TFl ( ! ; N 1 ) , T F 2 ( 1 r'Nj )) 

TFni.EO.Nl ) GO TP 261 
lid / 5 poo 
U cLL*5000 

IFfll , Ne’.L )G 0 TO (ISO 
COnTI NliF 
00 260 Na 1 , N i 
WRI TE (6, a? 1 )n,V(NJ 
wRTTE(6,a31 ) 

tVTSTC^NJ^Ma[^M l ( ) 1 ' N) ' U(M#N5 ' V<M ’ N, ' P(M ' N) ' POfM ' N, ' TfM ' N) ' H(M, ' N >* 
CONTINUE 

TFf ri.E'O'.Nl JGO TO 551 
FORMAT ( 1 0*. *N0, OF ITERATIONS*. 15) 

.At.OKor; HOOV'.SX, 'SHOCK STANDOFF otst *, 5X 'SHOCK 
tSPFFO »5 X, STANTON Nn, »,5X, 'SKIN FRICTION*. 5X, * w ALL P » ' 

ZSX.'wALL HPATTNG.W/H2*) ' L ' ' 

icv^e-«c'^i» ,P1 ^ ^ * 10^»F10,5, 10X,F10,5,6X»F 12.6, ax, FI 1 , 7 » 0 X , F 1 0 .5, 

O k #f 15, 3 ] 
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<121 FORMATS ,5*, »R00Y STATION N0.«»,I3,5X,'DIST, ALONG bODY* * , F 1 0 , 5, / ) 

aao format M x.flFi o,5) 

<131 FORMA TC/.fcX, 'ZN',9X,»U',RX#'V',flX,'P',BX,»R0',<>X,»T'-« , X,'H' , ,-'X, 
l*VIST») 

<J50 l.sl.+ l 

IF(l .LF.LMAX1G0 TO l 
551 CONTINUE 
OF TORN 
P-ID 



siiqrni'T i m? 

I nu P^f t 5) #t>PP(15> *''l rivIS) #RLEM5) ,U€( IS) »CtlT2Clon 
rnt-< M 0 'j/F i /N i , M u ,»'5r ,*52. * i . m 11 , ^sn 

COMMON /F 5/r,M, ,G m ?.GM3.C m 5,G u 6,G m 7,GM« , G U P , B2 , RF . TW , ' ‘“A, S!G, M* 

C : r'CM /r 11 / YX , ZZ , F» , HET \ , D v , 'l 

<;nH"(v«/r7 /i in ot , is) . vn ">i , t si , T c m , i *5) ,pn m . i5) ,pnu oi , is) * 
tVTS<!1C',?5>,n(i<M»15),SM(101,tS)#PMl5),vi5T(lM,15) 
CrwiO'i/FF/HT f 1oi,(S),EMDM1Pl,l5),ZNMni,t5),Rinoi,15) 

*OM>iON/F 1 O/TPf'C f ' •'I » 15) , Tt«rS(l 01 , IS) ,TVl ( 1 0 1 « 15) » Tv2( 1 01 , IS) 
CfiM u O>:/r 1 1 /|)I ' 1 •11,15), OS H( J 01 , 1 S) , A6M ni , 15), ab( 1 oi , ! 5) 

CDHMUM/r - * S/N 7 NF,Wtl ,»E2,S',.- 
NT \ =MT -1 
1 

'V' 3075 NS^JT.MI 
,) n i o "’o * * i s . M ' ; 

TFfHCM.Ki; , 6 E.HlNF)f.O TO 3072 
3070 PJ T T V'Jt 
^ r , 72 h_F(N)sZN(M,N) 

UFtN)5U(M,N) 

PF f N ) s P ( m , N > 

lf(N,ff],l!GO f. 3075 
DO 307 » 

J1=J~1 

307 j nin(N)=Pi 0(K W2.»ti)(Jl ,N)4il(J,N))/UEm )*(ZN( J#N)-Z 7 J1 ?/2. 
3075 ro^TIMUf 

Dn 30"5 mbNT?,niu 
300,5 DPE (N)s (PF (n* 1 ) -PE f N-l ) ) /** 

D n F(Nl jafPEf Ml )-PE fM1<«) )/XX*2. 

DO 3110 N=NT?,N1 

CUT 2 ( t Y m * .♦5,S*CZN(1 f u- Ml) /RLF (M) )**6) 

7V?(1 ,N jMl )e»F2*H0( 1 ,Nf MJ )*|)g (N)*BLDfNi ♦CllT2(l »^1 ) /VISC1 ^)M1 ) 

TVlUiNJcO, 

CUTTsSQPT { VT? ' 1#M)*D'K1 ,M)) 

CUT2( jMI )sVAHS(OU( 1 ,N|M1 ) jCUT 2( 1 ; Ml ) ) 

DO 3115 ms 2 ,mi 

CUT=CUTT/SOOT(ROCM,M)) 

UTAUcCUT/Rfc'1 
VPlUSeVM ,N) /UTAU 
VPLUSs-VPLUS*5,P 
APLUS326'.*EKP(VPLUS) 

CNPLUSsZN(M # n)*«0(m,n)*RE1/VIS<m,n)*cUT 
AMPIUSsCMPLUS/APLUS 
ir(ANPLUS.GT'. 12 .ir ,0 TO 3116 
VI8TLB.a*ZN(M,N)*{1 .-1 ,/EXP( ANPLUS) ) 

GO TO jllS 

31 16 VISTl*.<l*ZN(M,N) 
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3115 TV 1 CM,N)«RE*RO(M,Kj) # viSTL*VISTi./Vi8(M,N)*CUT0(M) 
313T mr b m*1 

DO 3130 Ms 1 , MR 
DV!9aTV2(M,M}„TV1 (m,n) 

IFCOVIS.LE.o'.IGO TO 3135 
31 30 VIST(m,kj)bTV1 (m,n) 
msm+1 

3135 MR1=M!-M + 1 

3110 VIST(",*J»MBnaTV2(M,N|Mflt) 

RETURN 

END 



SlIRPOUT I NE SHOCK 

3TMENSIOM *?(15) » VMF (15) 

COMMON /F t /►.'! ,N1U,toSO,to52,Mi ,MJ t ,h«5 0 
COMMON NX M2 

C0hm0N/F3/.';m 1 ,r,H? # GM^,r,M5,GHft,r,M7,GMB,r,MQ,H2,RF,Tw,G* u A s:g,hm 
COM 14 OM/FB/XX,ZZ»F m »PFT4#OY,02 

C0MM0N/F5/r>'l5),SSM5),r.(1S),CUR(15),US(1S),0sn5),VSnS),VKi(15) 

l.RMS) 

COMMCM/F#>/IMTtl5)»TI> , P(15)*TF.* 10 ft5)*T^tl^>#TFl(15)# 1 ’E2M5)*TEJC15) 
1 .TEU(IS) ,TE5( 1 5) , T^iSdSI ,TE7(l‘d 
COMMON/F 7/"Jf 1 0! , 15) , V ( 1 0 J , IS) ,?MM , 1S1 ,P(I 01 , 15) ,90(101 , IS' , 
1VJS(101. IS), H(101, 15), S M (tftl,lS)»Pto(l5),VtS '001,151 
COMMON/Ctt/RT n ot , l*i * ,FMOA (1 01 , ' 5) , * l p 1 , IS) , »\ ( 1 0 t « 15) 

COMMON /F P/Si ( t 0? , IS) ,S51 ( 101 , 1 5) fOS 1 ( 1 M , 15) -CU»1 ( I Cl , ' .3) 

C 0 M '* 0 n / F lO/THFCnoi, 1 S),THFSC 0 1,lS),TV1(10t,'5),Tv2nOi,15: 
common /f i?/Tz?notns),zinoi,i5),Tzi{ioi,i5),TZunoi i i , j) 

A?(n^i ) = ssm;ni )*cupn»M n 

c n ; n n e i , ♦ a ? m j n n * a 2 ( t > n o 

G ( 1 I N 1 )eVSQR.’(GO |N 1 > I c; C 1 J M 1 )) 

TFMP 1 1 f N i )afJ»VRATM»(P(Mi,t|HXMn t !NT(l|N! ) r TF.MP ( 1 | to 1 )) 

TFMPdjNl )sTFMP(1,ni ) / C- m 1 

TfdiM )s(TfMPn,>i| )*r,M7)/n.+TEMPM jni )*gm7) 

VNFtljMl )SGMB*(TE"P< 1 ,M )*r;M7) 

VNF(1|VJ )SVSORT ( VNF ( 1 1 Ml JfVNFfljNI )) 

US(1»N1 )SVNF(1|N1 >»fVN(1jN1 )«VS(I|'M )*A2(1|Nt ) )/G(l |N1 ) 

TF 1 ( 1 J N 1 )s(l ,-1./TF(1 »M1 ) ) * VNF f 1 i n 1 )/GOlNl ) 

TE?(1»N1 ) =VS ( 1 » n 1 >-TFiCt|N! )*A2(1 jN| ) 

TF3(l|Ni )*VNf t jNt ) + TF1(J|Nt ) 

TF<1 O | N 1 ) s T F M P (1 jNI )/TEU|Nt ) 

Tf MP ( 1 f M 1 )=VSGRT (TEttd Jtol )»TfMPdrNl )) 

TE5(1|M JaTPttnjMt )*TEMPfljNl )*GMfc/(TF<ltl >N1 )+GM5) 

TF6(1»M1 ) a Tpa ( 1 » M t )*Gmj 

TE7(1 |Vn»TFfc( 1 rMn*(TE2f 1 |Nn*Tfc?(1 jNn+TFK 1 |M1 )*TE5(1 »M) )/2, 
R 0 (M 1 , l|NXMi)aQPVSC ATP f TF f t ,N1 ) , INTO |M ) J RO ( M 1 » 1 l NX w ' . ) ) 

U(Mi, J |NXMl )sOBvsrAlW(TF2(1 »N1 ) f I bT ( 1 i to 1 ) 1 1) (M 1 , 1 | NXM i ) ) 

VJMt, 1 |NXMn3QflVSCATR(TE3(1 »M > ,INT(1 jf-'1 ) I V ( m i , \ , NXM 1) ) 

T(Ml#l|NXMl)aQ»V$CATRfTF4(J |N i ),ItoT(1ftot ) I T ( m 1 , j , NX* 1) ) 

V I S ( M 1 t 1 | M X M 1 )af3BVSCATR(TFR(1 | n ) ) . JtoT O |N1 )|VlS(Mt,1tNXMm 
5M(M1# \|NXMnaORVSCATR(TF6M f M t >,INTd»Nl ) » SH f M J , 1 | to t M \ ) ) 

H(M1 , l |MXMnaQRVSCATR(TE7( 1 |N1 ),lNT(J|tol ) | H ( m I , ! . NXM !) ) 

IF(MM,E0’.2)G0 TO 70 

0S(1|M1 JsUSMjM >*GO|Nl ) 

GO TO 71 

70 OSdJNl )a (ns d j tot )«ilSfl)N| )*Gd»Nj l)/2. 

71 Sd'lM )8S(t?Nl )*OS(l|Nl >*TIMPd|N1 ) 

00 75 Nat,Ni 
81(i»N|Mns3(N) 
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S8l(l#N|M1 )sSS(NO 

nsi ) = o s t ^ 3 

ZN(i # l |NXM)CS1 (1 , 1 »Z) (1 , 1 jH*H) ) 

P'Wd, 1 jNX M )=1 , + ZfJ( 1 , 1 ( 1 , 1 f NXM) 

RT(1 , t |NXH)SR1 (1 , J |MKM) + ZN(1 # 1 ,»iXH)*THgC ( 1 . 1 JNXM) 

RETURN 

END 



SUBROUTINE oerv 

COMMON/M /Nil ,N1(l,N‘jO,NS?, H l,Mt1,M50 
CnMM0N/F2/N¥u,NX'«il,NXM’S,NXM«,NXM5,*gXMfl,NXM2 

COMMON /F 3 /G“i ,GM2,GM3,GMR,GM6,GM7,GMB,GM9,B2,RE,Tw f GAMA,SIG,MM 
COMMON/Fd/XX, ZZ,F‘ < ,HFTA,nV',nZ 

COMMON/F7/U(l01,lS),VflM ,l5),TCln!,iS),PM0l,l t >).PO(10t,15I, 
ivis(if>i»is>,Mnoi,i5),srt(ioi,i5)*p^(t , ii»visTnoi,is) 
C0MM0Ni/FA/HTnnt«1ST f EHnt(1At,1S).ZMrtni,15).Rl(101«lS) 
fOMMON/t<?/St f 1 01 , 1S7.SS! (lot , 1S)»0S1 M 01 , 15), curj not , 15) 

COMMON/FI 0/THFf ( 101,t c i),TMFnno5,1S) < TVl(101,15),TV2(101,lS) 
COMMON/F ' t/r>U(10t.l*5),OSH<101.t5>#A6n<ll,15),AR(10t»iS) 

COMMON /FI ?./M?A 1 0t.lS5»Z1(1ft1,l5),TZ3fl0 1,ri),TZU(t01,i5) 

TV1 { 1 1 : 1NXM-I3-ZZ*S1 (1, t }NXM 5*722(1 , 1 jNXMi/H? 

out?* 1 inxm’jI =(U(3» t »NXM3J«U( l , t iNXM3n/TVl (2. 1 »nx*3) 

0SHC2? t |NXM3)e(SHfS,l|NXM3)»SM{t , J JNXMJ7J/TV1 (2 # 1 JNXM3) 

00 70 No l , N l 

0U.l!.N)3(fl,*W(2,N)»3,*L'(l#N)ot!(5,N))/TVl f 1 , N 5 
03 MU t N)c(U,*SH(2,N).3,*SWf 1 ,N)«SH(3, k n J/TV1 U '*) 

O'JfMJ ,N1-3.*IJ(M1 ,N)«SI(MS0,M) )/TVl ( ^ 1 ,N) 

DSHU*. ,N)*-(«,*5H(M|t,M)»3,*8M(Ml,N)-SH(M50,Nn/TV1(Ml f N) 

CONTINUE 

AfcU - 1 tNv»)-vTS( \ , 1 jNXM)/RE*(OIIU , 1 >NXM).W (1 , 1 t NXM ) *CUR 1 f 1 , '.'NX' 1 )/ 
1FM0A(1 , \ }MXM) ) 

ARU » 1. fMXMlef 1 ,«ZJ U , 1 »NXM) )/S1 (1, 1 ;N*m} 

RETURN 

fno 



subroutiwf bound 

COMMON/F 1 ,M<l»NSD,N5?» M t f mi 1 » M 50 

c o*-i MfjM / f / r.M \ , <■>*? 3 . o* m b » r^*>, G k1 7 , o>* n , g m r . b ? * » T w « 

COMMON/FM'J«X,.ZZ»F m #AE-'*'' V '‘'4 

i,TE«(lS^,"FSM5>,TE<'(ISJ,TF;7M'5: 

COMMON ^F7/ , - i rii')i#lS)»V , jf''l.,lS5»T(10?*lS^rF > C101 # l55» t *OC10J»lS^» 

JVTSf "31 • t 5 ) ,M( 1 01 , 1 ,SH(l ni , 15} ( t'O . VIST ( t 01 , l SI 

T E M P ( 1 j N ‘ 1=(3nvGATHP( 3 ( ?-. 1 t kJ X’M) , TNT ( t iNi )jTE HP (ll K 'l )) 
rw . i.uVui>.PIl^P*rt)fTf“9|1l , l! 1 . I N T f 1 i F' 1 )»P(l#t|NXMD) 


GAMA,S IS, mm 
T t'2(lFj f 7E3(15) 


Pt'. 1 { i • N1 ) = T F^P t 1 1 N t ) 

TFMPU|N1 ) =TFNP ( 1 5 N1 1/Tw/nMj 

Ron # 1 |N*m; ) sr.PVSGAT^CTFMPt 1 |N1 ),TNT(1{N1 ) j RO 1 1 , 1 f NXM 1 ) ) 
PET-'PN 


F NO 
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SUBROUTINE VECl 

DIMENSION 40(101. 15), *10(101, 15). A13(l 01. 15), AAl(lOlitS), 

1 AA2( 1 01 , 15) 

COMMON /FI /Ml ,Nl«,N l S0,M52,M1 ,H|1 ,M«jO 
C0MM0N/F?/NXW # NXM1 # MXM3,«lXw/i,Nyf'5,NyM8,MyM? 

COMMON /F 3 /GMi ,r,M2,RM3,r,“r, <( ;M (>f rif ,7 tR Mf»,r,MQ,n2,RE,Tw,CAMA,SIG,MM 

C0MM0N/F7/UM 01 , 15), VfUM . 15) ,T(tr i , mi ,p(( 01 , 1 5) ,l?0( 10 J , ; 5) , 

1 V I s ( 1 0 1 , 1 5 ) , H ( 1 0 1 , l 5 ) , S H ( '■ 0 1 , 1 5 ) , r ,■/ M »j 1 , V I S T (1 0 1 , t 5 ) 

COMMON /Ffl/UT f 1 0) , 15) ,F'*0 A ( 1 01 , 1 5) , 7N ( 1 01 , IS' , 01 ( 1 0 t , 1 5) 
CnMM0N/F9/$1 ( t 01 , 15) ,SS1 ( 1 01 , 15) ,nsi (1 "1 , 15) .Cl’Rl ( « •.' 1 , 15) 
COMMON/PI O/TwA'C J i 01 , 1 5) , TM? o (1 01 , IS) ,TV1 ( 101 , 1 5) ,TV3( l ?1 f 15) 
C0HM0N/F11 /minoi , 15),03M(1 01, «5) , Ah ( l 0 1 , 15) , AB( 10! , 15) 
COMMON/Fl ?/TZ?r 1 o; , 1S>,Z1 ( 1 C‘. , 15) ,T73noi , 15) . r 7U> 1 01 , 15) 
COMMON /F 1 «/A|l( 1 ot , 1 5, (J) , AM f 1 0 1 , 1 J f 4 v , j „• t 1^, j f A c ( 1 0 1 , 15 , V ) 
COMMON/Fl 6 /M INF, PF 1 ,RF£,S TOT 


AA1 f 
AM 1 
TV 1 ( 
AM ( 1 
AIK 1 
A 11(1 
AIK1 
TV2 ( 

AM ( 1 
A M ( 1 
A M ( 1 
A M ( 1 
AO ( 1 
A 1 0 ( 
AA1 ( 
A 1 3( 
1 NXMil 
A N ( 1 
A0( 1 
IN'XMiJ 

A A? ( 
1) 

A13( 

1NXM4I 
AN(1 
AO ( I 
1NXM4 


3TZ3 ( 
AA2( 
A 1 3 ( 
1NXMU 


, 1 J N X M > s V T S M , ' • M X M V/ 0 F * V T S T ( 1 , 1 j -J X M ) 

1 jNXM)s*A 1 (1 ,’ 1 |NX“)*ni'( 1 , 1 JMXM )♦*},( ] , 1 1 NXM1 

, 1 > MXN ) eS 1 ( 1 , 1 » nxm ) *FMOA f 1 , * ; N*M) *Rn ( j , 1 j « AM J 

1 , 1 J N X m i - T V 5 ( ] , 1 j M X m 1 

1 ,?tNXMisTV1 ( 1 , 1 |NXM)*U( 1 , 1 

I , 3lNXM>srTtfl (1 , 1. |NX«)*v:r , 1 ; l 'X“) 

1, aiNXMi f Tvi(i,iiMXM)*(H(i,ir*iyM)-r'i,i r NXMvROCi,ifNXM)) 

, 1 »NXM ) aSl ( 1 , 1 »MM«1 f.QO( 1 , 1 5 N X • ) *U( 1 , J !'-X m ) 

1 , 1 |NXM)sTV? (1,1 jNXM) 

J, ?}NXMlaTV2(1,1?MXM)#')(1,llNX M )4ni(l,lrNXM5*p(l,li^XM) 

1 , 3fNXM)sTvS( 1 , 1 |NXM)*V( 1 , 1 U'XM) 

1 , U J N X M ) r T V ? ( 1 , 1 » N X M 1 * M ( 1 , ' [• V x M )' 

?|NXMa)sFMf'4(i,?|Mvy«>„7Hrsn,?;NXMa)/BT(l,2iNXM0) 

,2»NXMii)s$t(1,?|MXMU)*TMEC(1.2;NXM«)/BT(l,Z|NXMU) 

, 2|NXM4l)sEMnA(1 ,2jMXM4J}*»0(t,?:NXMU)*V(1,2?NXMU) 
,2|NXMU)s*B(1,?j k XM0)«(0S1 (1 ,2jNXMa)*AU(l,2,l|NXM0)+SSl(l,2| 

*AM(1 \ t NXM<!) 1-AA1 (1 ,2jMVM/j) 

2, l \ nxmu)st2<k i ,2f wx Mii >* a - 3( ’ , r t nx h«) 

2,l|NXMal«A9(1,?|NXM«l*AMfl,2,l,NXMa)+A10(l,2|NXM«)*AAl(1.2j 

»T23(l,2)NXM«)*A13(l,?|MXMil> 

,2tMXM01sAA1(1 ,2»NXM«)*U(1 ,2,mxm«).emoa( 1 ,2 f NXM0)*A(,(1 ,2jNXMa 
,2)NXMai*AB(i # ?|UXMa)*(0Sl(l,2;NXMa)*AU(l,2,2jNxMO'*0Sl(l l 2| 

*AM(1 ,2,?lNXMil) )-AA?(l ,2jNXMUl 

2,2|NXMtt)*TZa(l,?|MXM<|)*Ali(!,2,NXM<n 

2,2|NXM|J )SAQ(1 ,2tNXMij)*AM(l ,2,2|MXM4)«Alft(l ,2|NXM 0 )*AA2( 1 , ? ( 

+ TV?fl ,2|NXMfli*c"R1 (1 ,2|NXM<n*VM ,2,NX M 0)-S1 ( 1 . 2 t nymu )*CUR1 ( 

A M t a 3.KlvUfi\AiO/'i n - iu w u % . . 


in'Ar'iM/'Tivri i 1 cjimhu l w ^ I I n 1 * I 

21 ,2»Nxm«)*a#)( 1 ,2|NXM01.P( t ,2|NXM0)*S1 (1 ,2|NXMa)*AO( l ,2)NXM(>) ( 

1T71M \ aMIH 5.KIVMJI « 


,2jNXMui*a13(1,2|NXmu) 

,2)MXMU)aAAl (1 .2»NXMO)*V(l,2':NXM«)+EM0A(l,2|NXMO)*Pn , 2,NXMil) 
,2»NXM«)eAA( 1 ,2jMXM«)#(0Sl ( 1 , 2?NXM4>*AU( l ,2,3?NXM4)+SS1 (1,2? 
*AM(!,2,3»NXMtt))-AA2M,2jNXM«) 
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*M(i#2#S|NXt*<naTZflri,2|riXMa)#AlS(l#2»NXM«)' 

TVt (l#2|NX^«^BA«rt ,2»NX^a)*AMM ,2,3 10(1 *2»nX«a)*AA2(i,2f 
lnjXMa)»TV2O,?jNXM0)*cnRtn.2»N* Mt n*nn»2i^x M a) 

A0(1 ,?,3jnx*uisTvi f 1 , ?»NX M U> -P(l,2t^xM«)*(Ain(i #2 iNXmu ) *ehda ( 

31 ,2»NXM<n*Si (1 ,2;NXMUl*riJRl (l,?»NXMan«TZ3(l .2|NX^a)*Al3(l,2»NXHfl) 

AA?n*2|NX M ui3AAtn.2»^X MU )* H ('*^1 , ' JXMy 5-F M 0A(1.2|NXHa)*(VI8U»2l 
lNXMUl/Slr./RF* ( 1 •♦VIST ( 1 # 2>NXMU1*SIG/Sir.n*0SM(l ,2»NXM0) + 
?U(l»2!MX**ai*A6(1.2i w xMan 

'a13(1»2|I'jX m “)s*H( 1 ,?|MXMU)*(nsi C 1 , 2 jMXHUi*AU(l, 2 , 0 »NXMa)tSSl (1 ,2» 

lNXMa)*AMfl.?,ajNX M U)l»442n,2l l ‘JXMol 

AW(l,?,iJjMXM/ji3TZUMf2j N XMiji*Ai3(l,2jMXMm 

AO(i ,2,a»!VXM<l)sA9fl , 2 |NX*<n*A*M , ? , 0 , W X Ma ) ♦ A 1 0 ( 1 , 2 | N x«u ) * A A 2 (1 , 2 I 
lMXMUI-TZ3n,?;NX«U)*Al5(l#2l N,XMU ) 

RFTUKN 

END 



SUBROUTINE VF C2 
DJMFNSIOM A?0(I01) 

COMMOM/Ft /Nt , Ml <1,N50,N52, M 1 ,M1 1 , M 50 

COM^ON/FVGMI ,GMp,r,M3,GMS,r,M6,r.M7,G M «,KMq,B2,Rt »TW, GAMA, SIG,MM 
COMMON/F7/lJ(101 # 15),V(101,lS)#rnoi,lS).P(101,lS),ROUOt,15)f 
lVIS<101,15),MC101.1 , *),SMd01,15),Pw(l5),VTST(10l,lS) 

COMMON /Fft/BT (1 01 . ISJ.EMOAM M , 1 *5 J .ZNM01 . IS) ,R1 (1 71 , IS) 

cnMMt)N/FR/smoi,is),ssi (ioi , is) *osi M oi , is) ,cu«i f ioi , 15) 

COMMON /FI l/OUM 01 # 15)#0SM( 1 01 , m) , Afrf 1 ni . 15) , AB(1 01 # 15) 

COMMON /F 1 ?/T?2( 101 , 15) ,71 ( 1 , IS) , TZ3( 1 01 , 15) ,T7ii M 01 , 15) 

COMMON/F 15/ A 110(101 ,U), AMO( 1 01 , ?,U) , ANO( 101 ,«) , AQOdOl ,11) 

COMMON/F 1 fc/HJNF ,Rf 1 ,RE2,SIGT 
DO 10 0 I e 1 , 2 

A?0(l : M 1 )5S1 (1 , I* Ml )*R0(1 ,T»M1 ) 

AMOf I , I , ) .Ml )a2,*A?0(l JM1 )*IJ( 1 , T 1*1 ) 

AMOll , I,2jHllsA20(1 |M1)*<P(1, I j Ml )/RD(1 , I I Ml )+2,*IJ(l , !»M1)**2) 

amo(i, !,3jMn 8 AMn(i , i,i }Mn*vn , ij“n 
100 A m o ( 1 , ! , U I M 1 5 a A M 0 (1 , I , i * M : } « H f 1 , 1 ; M n 

A 2 Of 1 1 M 1 )sS« (1 , 1 ;M 1 )*S0C1 , 1 jMn*FMOA(l , 1 j'mi ) 

Allot 1 # 1 I M J )SA?0( 1 JMII 

AUP ( 1 » 2 j M 1 ) SA 20 ( 1 ) M 1 ) *l)( 1 , 1 t M t ) 

AMO (1,3 | Ml )SA20(1 )M1 1*V(1 , 1 (HI) 

Allot 1 ,«|M1 )SA 20(1 1*1 )*(H(1 , 1 |M|).C(1 . 1 fMI )/RO(l , 1 JMI) ) 

A?fl ( 1 I Ml )sfcMriA{ 1 , 1 |Mt )*R0( 1 , 1 > M 1 ) 

ANOf 1 , 1 |*M )SA?0( 1 I Ml )*v ( 1 , 1 |M1)«A8( 1, 1'fMl )*(0S1 (1 , 1 »M1 )*AU0(1 , 1 | H 1 

1 )»SS1 (1 , 1 1*1 )*AMH(1 , 1 , 1 t Mt ) ) 

•imp:'. ,?| M 1 ) e A 2 0 ( 1 } M 1 ) * ( U ( 1 f 1 fMi ) * V (1 ,1 jmi)«a<>( 1,1 I Ml )/RO(l , 1 » M 1 ))• 
1 AA( 1 » 1 I Ml ) *(0S1 (1 , 1 »M1 )*AUf1( 1,2 »M1 )*SS 1 ( 1 , l »M1 ) AA'lfU 1 , 1 ,2jM) ) ) 
ANO(1#3lMn=A20(1|Ml)*(V(l,t|Ml)*V(1,lj- < n+P(l,1|M1)/RO(l,ltMn)- 
1AB( 1, 1 } Ml )*(0S1 (1, 1 )M1 )*Alin(l,3»Ml)*SSl (1 , 1 »M1 )*AMO( 1, 1 , J|M1) ) 

A NO ( 1 # U|M* )SA?0(! S»1 )*(V( 1 , 1 |Mt )*Hf 1 , 1 fMi ).VIS(1 , 1 I’ll )*r>SH(l, 1 J*l ) 
J/STG/Rp/ROd , 1 JMI )*H(1 ,1 1M1)*A6(1, 1 » Ml) /ROT 1 , 1 1 M 1 ) ) -Aft ( 1 , 1 1 M 1 ) * ( OS 
21 (1 • 1 |M - )*AII0(1 , U,M1 ) ♦ SS 1 ( 1 , t J Ml )*AM0(1 , 1 ,<lin )) 

A303CU^1 (1 , 1 )*S1(1 , 1 ) 

AOOni1lM1)eA30*Hn(1,1|Mn*V(l,lVMl)i,ANO(l,ltM1)*TZ3Cl#l» fc ‘n 

AfJOt 1 ,?>M1 )o?.*A30*(R0 n , 1 f Ml ) *u c 1 , 1 |M1)*V( 1 , 1 |M1 )„A6( 1 , 1 |M1 ) ) ♦ 
1AM0(1,2jM1 ) *TZ3( 1 , 1 1*1) 

AQ0(1 , 3»M1)8A30*(^P(1 , ljMl)+Rf1(l, 1 j m 1)<.(V(1,1iM1)*V(1,1jM 1)« 
iu(i , i » m i )*u(i , i>Mnn*ANn(i,3jMi)*TZ3(i, i jmi) 

AOO(l,fl»Ml )KA30*(RU(1, 1 1 mi )**/( 1 , 1 JMI )«HI| , 1 JMD-Visn, 1 »M1 )*OSHC 1 , 
1 1 ;M1 )ZRE ZSIG-U(1 , 1 )M1 )*A6(1 ,1 |M1 )ltAN0(l,4|Ml >*TZ3(1, 1 )M1 ) 

AM0(1, 1 |M 1 )s«rza( 1 ,1 |MJ )*AN0(1 , 1 |M 1 ) 

. ANOd ,?»M1 )s-TZ«(1 , 1 »*M )*ANO( 1 ,2)M1 ) 

AN0d,3|Ml )s-TZ« ( 1 , 1 t«l )*ANOd ,3»M1 ) 

ANOd#OfMl)B«TZa(l , 1 |M1)*AN0(1 ,«|M1) 

RETURN 

END 
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* SUBROUTINE VFC3 

COMMON /F \ /U 1, NlU, NSO, NS?, MJ ,H||,M50 

COmmon/f?/n>:m i n*m),nxm'j,nxmo,nx^S,nxmb,nvm2 

COMMON/F3/CM1 ,5Mj,r,M3,r,M^,GM b .r,M7,RMH,G M 9,»2,RE, Tw, G4“*,SIGf MM 
COMMON/Fb/INTMSI , T IMP MSI # TEMP f 1 5 ) , Tg (t 5) > TE 1 (15),TE2(15),TE3(i?) 
I.TFU(lST r TF ,: :ciS),TFh(1S),TF7(1S) 

cnMM0N/F7/uri oi , isjfVnowisj.’npuisj.PMoi.isi.ROdoi^is), 
lVISn <, l,15 5,wnoi,JS),SM(liM,1S)#PW{tS),VISTftnt,iBJ 
COMMON /r "■ T ( 1 * 1 , 1 S ) , f'ip 4 ( \ 0 \ , 1 5 ) , ZN (i 0 ! , 1 5 ) , R i ( 1 0 \ , ':*■ ) 

COMMON /f 9/Sll t 01 , 15) .SSI ( 1 01 , IS) #nsi f 1 01 . 1 5) ,CURt n 3t , 15) 
COMMON/FIO/TMKCno? , 15),TMCS< 10) , 1 5 1 , T V U 1 0 1 , 1 5 ) , T V2 (1 0 l , « S ) 

COMMON 'E 1 3/ ill 1 (101 . ) 5, a > 

R0(?, 1 »nxm?) = aiji (2, 1 ,« ,nxm?)/si (2 1 JNXM21/EM0AC2, 1 , NXM2) 

U(?»2|NXM*j)s4IH(?,S t ?,NXM5)/ai|1(2,2,tJNXM5) 

V(2,! |NyM2)sAUl(?,t # 3»NXM2)/4iiu2, i , i jnxm2) 

TF(1 INI 1=IJBVC,4TH»( V ( 1 , 1» N XM1) , INTM • Vl) ,TF (1 JN1 n 
TE1 (1 »N1 ) S ‘3flvGATMR(V( 3, 1 ;NXM1> , INT ( 1 ,N1) , TEl 1 1 » ^ 1 ) > 

7E2(1 »N1 )s(3'.*TE(i |Nn + TEl (1 |Ml) )/«, 

V(2. * >NXMnsf38VSCATP (TF?( 1 iNn , IHTM JNI) , v(2 # t iNXMO) 

TV1 (2, 1 JNXM?1 = (C(2, 1 ;'!XM2)*ll(2, 1 J MX M ? ) ♦ v f 2 , 1 |MXM2)«V(2, 1 |NX w 2))/2, 

SH(2« 1 inxm?)s( AU 1 (?, l ,«}NXM?1 /aiji (2, 1 , 1 }NX«2)»tv1 (2, 1 >NXM?) )*GAMA 

M(2, 1 |NXM?)aSHf? # l JNXM2)4TVU2,I»NXM 2 1 

T(5»1 |NXm2)sSH(2 # l t mx m2 1/GM3 

P(2» 1 »MXM2)rGM9*(}0f ?, 1 » k * M 2)*SM(2, 1 jNXM?) 

TV?(2,l|NXM? 1 -vSGRT(T(?,i,Nx u 2)}T v ?C2«t|NXM2)) 

VIS (2# \ jMX^213T(2,1 >NxM 2)#TV2(2, 1 jNXM2)*GM6/(T(2, 1 |NX m 2)+GM5) 

RETURN 

END 



StIRROUTIKE. VP CS 

COMMON/Fl /N 1 ,Ht «,NS0 , N5?i M 1 , HI 1 ,M50 
CnMMON/F?/NVM,NXMt ,MX‘*3,HXMU # NXH5,NXMP,MXM2 

ronMnw/F3/r,Mt ,c.M2,G^i,fi«5,r,M6,r,M7,c;^8#GH<»,B2,Pf:#Tw,GAMA,sTG#MM 
COMMOM/Fh/I^T (15) #TIMP{ 15) .TE^Pt 15) ,TE fl5) #TE1 ( IS) »T£2(15) ,TE3(15) 
!,TFU(15)*TE5f 15),T£M15),TF7(1S) 

cnM«ON/F7/»(1ft1 »15),V(jni,i5).T(lPl f i5),P(ini # iS),RO(tOt,15), 
lVlsnoi,i5),H(ini,l5),SMMni,l5),PW(i5),viST(10t,15) 
CnMMON/F«/RTM01,1S),EMftA(101,l5),ZNM01,lS) # Rt (101,15) 
COmmON/Fo/SI (101 ,1S),SS1 (101 ,1S) ,nsi not ,15),CURt (im , 15) 

CPMM0M/Fl<!/THEC(lO1,lS)»THF5(i0t,15),TVl(101,15),TV2(101,15) 

COMMON /FI 3/*IH(t01, 15, «) 

(?0(?» 1 |NK M 5)sA<Jl (2,1 ,1 I M*M5) / si (2,l; k JXM5)/EMDA(2,l|NXM5) 
U(?#2|NXM8)saI.I1 (2,?,2|HVM«)/AII1 (2,2, 1 jNXMfl) 

V(2, IiNXMSIsA 01(2,1, 3» NXM5)/AU1(2, 1,1) M VMS) 

T F(1 l KJ 1 ) sO« VC; A THP ( V ( 1 ,1 JNKMJ ) , INT( 1 T N1 ) t TE ( 1 J M 1 ) ) 

Tt'1 (1 )NJ )sOflVRATHR(V(3, 1 tNXM1),INT(1tNl)|TEl(l|N|)) 

TF2( 1 |^1 )S( 3', *TE(1JN1WTE1(1|N1))/U ( 

V(2, 1 )NXM1)sOftVSC*TR(TF?(1 t M ),INT(1 ;N1 )','V(?, 1 ,NYM t ) ) 

TV1 (2,1 j NXM5)s (U(2, 1 jNXM5)*||(?,i • NXMtj) + v (2, 1 |NXM5)*V(2 # l|NXM5))/g 

SM(? # 1 fNXMSls ( AMI (2, t ,«|NXmS)/MH (?, 1 , 1 jN’XM5) w TVl(2, 1 jNXMS) )*GAMA 
H(2, 1 jNXm5)s.SH( 2, 1 |MXM5>*TV1 (2, ljNXMS) 

T (?, 1 | NXM513SH (?, 1 »MXM5) /0 M 3 

P(2, 1 jMXM5)sfi«Q*«0(? # 1 1 'IXMS)*SHf 2, 1 ; MX MR) 

U(2, Ml , m 1 1 ) ! M11V2’, AIUP, M50|M1 1 ) ♦,5*U(?, N52 jm) j ) 

V ( 2 , N 1 1 M 1 1 ) = ?,5*V(?,N1U;M1 ) )-?.*V(?,'50|Ml j ) + a 5*>v(2,N52;Mi 1 ) 
P(2, N 1'Mil'l = ^,5*P(p > M0jMli)-?.A p f2,N50jMl?) + ,5iH3(p t ,j5p.v!- ( ,* 

SM(2,Ni)M{i) s SH(?,NH,)?l)i)*p t 5 w 2.A5l J {?,N50|Mii) + > ^tSM(2iN5?iMll) 

H ( ? , M 1 1 M 1 1 )sSH( 2,N1 | Ml 1 ) ♦ (II (2, Ml > U 1 1 )*D(2,Nl f MH)«.v(2,Nl)M*i 1*V(2, 
1 M 1 } M 1 1 ) ) / 2 , 

T(?,AO fMH)5SM(?,N1,Mn)/GM3 

RrifP,Ni »mi 1)sO(2,ni ,m, j )/r.Mo/SH( 2 ,Ni ?mi n ~ 

TV2(2,ljNXM2)sV5fiRT(T(2,tjMxM?)jTv?(?,ljnXM2>) 

VISC2, 1 |NXM2)sT(2, J jNXM?)*tv2(2, 1 ;nxm2)*GM6/(T(2, 1 jNXM2) + cm5) 

RETURN 

END 
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APPENDIX C 


SAMPLE INPUT 

In this appendix, various input parameters are prescribed for the laminar 
flow past a sphere cone. The first sample problem investigates the flow at 
5° angle of attack by using the code given in appendix A. The second 
problem investigates the flow past a hyrerboloid at zero angle of attack by ns in 
the code given In appendix £. 

Input parameter NBODY is proscribed in the main program EBODY. The 
rest of the parameters are prescribed either in BOON or BHYP depending 
on whether NBCDY = 1 or 2. 


NBODY 


l 

CC 

= 

.001 

FM 

= 

10.3 

IF 

= 

U 6. 26 

PF 

= 

100. TT 

ALPHA 


.087266 (5°) 

THC 

= 

.7851 (15°) 

TW 

= 

330.6 

RN 

= 

.03175 

DY 

= 

.185 

BETA 

= 

1.1 

FDT 

= 

1.0 

THESH 

= 

.03^9 (2°) 

LMAX 

= 

20000 

LTURB 

=r 

25000 

GRIT 

— 

.01 


Input for Sample Problem 1 



APPENDIX C (continued) 


Input for Sample Problem 2 


NI30DY 

= 

2 

0 C 

= 

. 00 i 

FM 

= 

10.1 

TF 





PF 


100.77 

ALPHA 

■- 

0 

m r ] n 

i 1 L O 

- 

.78olt 



O’ -J o y 

1 '/V 


JoU . O 

PN 

- 

* 021.75 

* A 


• i.A5 

p f/rn 

= 

. 1 . . 1 

F [)■'[’ 

= 

i . 0 

TKH2H 

= 

.0:1^9 

] 


20000 

I.'TTJRB 

= 

25000 

CRT? 

nr 

.01 



APPENDIX D 


SAMPLE OUTPUT 


The sample outputs are given here for the two cases discussed in appendix 
The code initially prints the body geometry parameters and complete starting 
solution. The shock shape, heating rates, surface pressures, skin friction, 
etc. , are printed after every 500 time steps. The detailed flow quantities are 
printed after every 5000 time steps and after the solution is converged. Only 
a portion of the complete output is presented here which includes the body 
geometry parameters, a sample of gross flow quantities and a sample of detailed 
t low quantities. Detailed flow quantities are given here only at one body 
station although these quantities are printed at each body station in the" 
complete output. 
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APPEND IX 7) (continued.) 


DJ ST 


finer 
f 6DY 
bcur 
6DDV 
i tj u r 
u m 
t( ur 

fiUDY 

fct'fJT 

fc 01) Y 

f nor 
finer 
fDDY 
K:DT 

fc fiU r 


Output for Sample Problem 1 


Y- 

• 0 V *. 

Tlii* 

* * i •' 


Y • 

. * -ilL'J 

TH: - 

. . ■ 

‘ t « j 

r • 

• j 7 ; i... ^ 

Tm . - 

A • < j 

i r J 

r • 

• *> 6 i. V< ' .' 

T M - 

! .6 1 


Y ■ 

• 7906 

1 - 

♦ : ■) 

. f ^ 

Y- 

• 'V r ... « ( . 

Till » 

. 7 

1 *'j 

Y • 

l.lh j * i 

Ti* - 

• 7 

! 

Y- 

i.* 

7*0. 

* 7 

: * ,) 

Y ■ 

1 • 9 :o ■:< 

■ 

r 

. 7 

4 ■'* ■ 

Y* 

1.690: 

TH .. 

. / 

c M - 

Y • 

1 . d JG’f’ 

T Hr- 

, 7 

‘ 

Y- 

2 • 0 J ; i ( 

7 0* 

• 7 

1 * »J 

Y- 

UU K <*. 

T M » 

. / 

*0 

Y* 

2**00 1 

pH.* 

. 7 

■ « '■ . 

Y • 

2 * v , ( • 

7m- - 

.Vo 

■| s \ 

Y* 

n i 'j ! 1. j 

r H ■ « 

i , : 


Y- 

. 4 7.. ; 

llfp 



Y - 

« 7 t \ 1 L 

7-i . 

j . 


Y» 

. /90j 

IH * 



Y * 

* 7 2 : 

U-t * 

. 7 

H J 

v . 

dfu''l 

i- * 

» / 

9ij 

Y » 

1.2V vc v 

Tf ; > 

• 7 

< v , 

Y • 

1 *S . v , 


. 7 

to 

Y *■ 

1 . 7. 1 c '1 

T -U » 


* i> 

Y * 

J . V - - 

7 H . » 

. 7 


Y » 

2 • J 4 1 C 0 

T H r « 

. 7 > 

A ■ 1 

Y , 

2 . 2 . t. c 

r n * 


•i 0 

Y • 

‘ * i .) O . Y 

y l i m 





1 • y . 9 


* i T 

Y- 

2 « ! 7 c. c G 

IH ;■ 

. 7 

t V/ 


•s • 

.•• 'HO 

G J 5 - ■ 

i . f < \ o ..* 

► • 

• ! v f l 

U-'ip 

1*G( > 

r ■ 

. i 6 1 t £ 

C 12 • 

1 « t G cp C ’i 

► m 

• !' 2 r f < s 

G , ■ 

: .( aoi 

K * 

.f 7-fc ■ 

it. r • 

1 • C t v C : » 

f t 

. - o 3 .. 2 

C r- * 

. fif'C J9 

' ■ 

.9 3ot ? 

( n • 

. C,r i-»r,r» 

► ^ 

i.ooV*: 

l wP rf » 

. C u u 0 'j 

K ■ 

i . ' 'V . 2 1 

l. 1 2 • 

* ( 00 'j. 

*■■ ■* 

! . .< ;. , ' t 

C :J h - 

• V C yj G 3 

L m 

X * * ■. *> S \> 

l L K 

• toot e 

t ■ 

j • * . ^ 7 1 

cu * 

• C: C v 0 3 

f « 

i. v «? ; > i. 

0 J V “ 

.root o 

' * 

i . j 7 2 3 * 

C l ► * 

• f o o C 


1 . 9" i 1 * 

'■ v " 

*( C ■ ^ 7 
1 .< 0. 2.V 

V' . 

• 3 o i ti 2 

: k - 

• «uf 

* " 

. ■; .:** * ‘ 

l : v - 

1 . • 1 ■ 'jV 

^ » 

k' » 

„ | ^ v 

• u )’ 7 

P 

) . C L f 07t 

• i. 0 i <:■ .> 

t • 

• ! > .' Vi 

( 1. i"S • 

• L'\ C *> 

r ■» 

> .uf, /*.' 

(, « 

. t l- c ; ' 

K * 

i . ; ‘>.<2 6 

L L * * 

.urn 

K r 

I 4 f<‘ »t t. 

C ;.'► - 

.Cr, 

r * 

1 4 .<‘7 

C. L' Y * 

• f*o :t.c 

A' * 

1 .1 t ’ '7 1 

l. ^ * 

,r roc ’ 

K » 

:* 7’ i >2 

( * 

• CC v-G3 

2 * 

. * . <.39 

C • J t * 

»f 'jv 0 Y 

H * 

i . 9 4; 3 1 

cue* 

*( 'vide;* 


. along. 0 : 1 c y 

■ OJiA 

• 1 1 s o ( 

« i 7 v ( vp 

. ^ccc 

• 9 i s i 1* 
i . l icc n 
i . 2 *: no 
i . 9 tail c 

1 • 6 £ I G * 

1 . d : 0 U 

2 . 03166 
2 .22000 
2 • 9 o M f 
2* OOuO 

• 1 1 S (' f 
. 3 7 a. if. 
.SS6( (> 

. 79U 0 

•92S6r 

1.1160C 
1. 2V>f.'( 
1.9 60c v 
l.bGSOf 

1. tiSl' r* 

2. 33ion 

2 . 2 2 G 0 r. 

2 • 9 c 1- 0 0 

2 • 6960 ( 


M AMOS 

V - 

k/iJlOS 

M * 

K A L I G L 

R ■ 

k ,\ G 1 0 S 

K ■ 

K AJ1U5 

v . 

KaJI t j 

R * 

kAfUUS 

f - 

RADIOS 

R ■ 

RADIOS 

R • 

PAD I OS 

k • 

RADIOS 

2 - 

RADIUS 

V . 

RADIOS 

f . 

RADIOS 

k - 

RADIOS 

¥ • 


Of k s 1 e,L f f f 1 s [ . 

• 10 6 7 
. 2 9 )*5 

• ;99 c /c 
. j / 2 1 0 
. 1 C7 ( 2 
.Uviv 

• 2 ;< l 3 V 

./Git*. 

t ? /. j t ? 

• 271. I 
« 2 r 3 3 o 

. 

. lu-'t'i V 

.3:7 ‘f 

• j 3 2 7 1/ 

,lrn 

• , C)jj 
. 11^4 

• I - 17 9 

• 2 c ;: c ! 

• 2 3 6 9 2 
*2 j Vt 

• 2 9 i : 3 
.269! d 

• 2* 9A 9 
.I'u'-jl- 

• 3 y k * c 

• 3 1 j. 9 9 

.37*79 


4 0 jvOC 

PKt SSOR . 

.1^39: 

F’RfSSL'iV 

4 3ol 62 

f ; K! SSiVR ; 

4 J£( * <» 

RRl- *1 sur l 

• 6 792 9 

F R t SSURr 

• >01 JA 

F R t S ‘ J2 ^ 

. i 3 1 62 

PRI SSORc 

l.JUM 

FK-SSORt 

l.lCt2fc 

Pkw ssu.< 

U 329t > 

P R J S S U k c 

1.9 S 69 

F R . S S UR J 

A 4 J vC 7 J 

f'*U SS JR . 

1. 721S2 

P R : S S U R i. 

A • J ! c 3 A 

PRtSS0«‘. 

1. U ilt 

PRcS^OR* 


SM- c K : "V! 0 

M ANT IM rv ; . 

• fj'JO 'Y* 

♦ C 1 1 0 1 7 9 '* 

,0^( f j 

.0)*1 1 ^ <, 

• D-’t-t ? 

. C 1 2 r- 9 7 ; * 

• o i ’ s 

*"9 U7972 

• 0*/tci-Y 

. C07DG S96 

4 ■;'>*» l 9 

4 1. c 6 9 J r. ! ' 7 

4 D.U 1 , 

4 6 f 1 S v a 2 0 1 

• t-'. 1 \ 2 

4 S', 9 ■ 6c9 

-4 ODOvU 

. C 9.322* 

'-vt ■' i 

• 00 9 2*. 9 06 

-4 D v'-'. 9 

« 60 9 *9 VC 

-•C‘fiOC 3 

• Of. 3 9 3 9 9 9 

-• COo -j 3 

- <V3'9 i S 9 C 

• J c o ° 2: 

* S'. 3 27 2 VC 

. o ; c o i 

.70S 9 70 J A 

♦ 0C0C.9 

* C I v o 9 21 3 

4 U Yt* U 

.olUCrd, 


• C 1 ! 7^ **o 

• O'. / S 9 

• oi.coieo? 

4 nn'io 

.0077*231 

j 9 

*0267? 7 f* 7 

• oom i. 

. ( ( M 9 3 * 9 

” .CSi.f j 


-•00u<.3 

. JC962397 

-.I'C.S '2 

. ft: 1 91 23 3 

-*-7.n ■ 1 

• C0.S229 J7 

-.COf’fi D 

39 f >- 7 

-4 iv >1 

. ' C 9 Si f t; 6 

• 0v<( I 

4 G O 9 ciu 7 1 9 


DRAG* 

3 . R26S2 

OPi G» 

4.733-V 

9 R A G • 

1.793t7 

OP AG* 

1. ^372 

i./ R A C » 

i . 2 9 9 ii 

DRAG* 

1.27 V ? 

l/KAO 

i, 1 P j 

ORAC* 

1* 13S26 

1 i R A G ■ 

1 • LfV 22 

Drag* 

x • d 9 3 2 2 

DRAG. 

1 4 JC 1 29 

li k A G * 


Drag* 

1.37939 

Lk ao* 

1.07293 

Drag* 

1.0771S 


N FPICIIC'N 

• ALL P r . 

*00 2 1 1 1 / 

.91^26 

• r>- ?* v.9c 

• 6 9 6 C 9 

. C 1 1 1 ? ;: 1 

• 7 2 2 3 v 

*vl?v*fci?G 

• > f ? i /* 

. u i } S i 3 G 

•9 19c 9 

• Ou* j 926 

• 9C):7 

•00/3637 

• Cl 6 3 

•Yjtsjil 

• 3 v 9 f; vl 

• ) o ! C 1 G 9 

. *.1. 7 

• 00 9 9 30;. 

• 9 1 : :. 9 

•00912/6 

.92719 

•00396 6 3 

. 9 3 > 79 

4 Oul'S 3 V C* 

.991/1; 

* C 0 3 9 c ! * 

• 9 9 * 2 0 

« D'j 39 vCf 3 

.99.39 

4 0o’2t 2 2? 

• V 1 9 f . 7 

» t t 1 7 ^ !> r o 

•09913 

^011 6 397 

.71200 

• d 1 1 2 7 1 6 

• : 6 )? 7 

.00 7 t 7 i 1 

. 6 76 -/1 

*C o 60 1 63 

• S 7 t69 

*0099912 

• S 9 3 1 

• 0 0 7 S 1 6 6 

•6C 7d2 

•009966b 

p.MSH 

• 6099 i G^ 

•62321 

* 0O93S33 

4 6 22b9 

• U09<96‘t 

• 623t>: 

•0092976 

. .9 2 3 S 6 

.C09I6S6 

.62269 
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Zs 

• OOOOJ 
,C3J9 J 

• COC? 1 
•C3I23 
. C Q x 6 7 

• C J 2 I i 
« W 0 2 l ’ 7 
, O j 3 0 1 
.00362 

.05402 
. Ju453 
.00506 

.00615 
*05t72 
•0373 J 
.007 y 3 

• jot 5 i 
,00914 
. 0 U 4 7 7 
. Hi 1 0 4 5 
.041X3 

. 0 1 1 3 

• 01? 5 5 

. 3 » i 2 >5 
a U x 4 j 5 
. 0 1 4 W 
.01 1 GO 
» 0 ,i i> S i 

.01 724 
.0 A to 7 

.01(96 
.01 6 
. 0?<i77 
7 j 

. O 2 2 f o 
.023'ti 

• 0 ? 4 f 4 
.02 6t>7 
.92 6 72 

• C 2 7 7V 

• 0 2 1 8 7 
•G30J1 
•03115 
.03232 
•03352 

• 03** 74 " 

• 0 3 5 ■? J 
. 0 3 7 2 *j 

• 0 3 t .5 5 
, 0 3 9 ° 7 

• CM2 
. 0 4 2 5 7 

• 043V* 

. j 4 6 y 2 
. 0 4 t 7 
. u 4 t 3 5 
. 04 9 tJ o 

• 0 5 1 4 0 

• 0 5 2 v t> 
•05454 
.05616 
.05777 
,05946 
_• 06 1 X 5 _ 
,06267 
•064ol 
,0663a 
•ooei a 

• 0 7 0 0 J 
.07184 
.07371 
.07560 
.07751 
.07945 
.06141 
.08339 

• 0853V “ 
.06742 
,08946 
.09152 

• 09 36 'J 
.09570 
.09782 
.09995 
•10210 
•10426 
.10644 
.10863 
.11083"' 
.11305 
•11527 
♦11750 
•11974 
.12x99 
.12424 
•12650 

• 1 2 b 7 7 
♦13103 
.13330 
. 1365 7 


L 

V 

WP 

r 

. 0 OG 00 

. Cl OUC 

. 0 OC 

.9x82c 

. oL 799 

.3.7 

— « j 1 ; v 9 

» 4 * r 2 6 

. L is 79 

-.uCC 56 

- ♦ 0 1 *1 74 

.9 1.27 

• o 2 C e 

X 025 

-.U'l V 

. 4 1 e x 

• 0? 6 11 

- . J «j G 5 2 

L C i ' 

• 4 1 • x -v 

• C3wr.2 

- • t C t o 7 

- • - 3 1_ c 

. 9 i »■ 2 v 

•0316, 

XI 2c 

a 

• 4 x ■. :• j 

• C 2 c 11 

-.X ) 77 

-.ao r? 

. 9 j r 3 ; 

.t 4199 

-. J5 2 : i 

.41 Vr 

. 4 x t j x 

, 0 4 4 t 5 

-•ol?*: 

- , j 4 *» : r 

,918 21 

.44732 

-•■74 3 56 

— • i 4 7 j l 

. v n ; j i 

, 4 5 *« 5 

Oi 4i 5 

-.i’44 4 5 

* 4 * “■ ' 0 

• 6 ’ 2 r 

-•Cf 44 t 

- *c, 5 12 ^ 

. 4 . t 2 - 

*(5277 

Jl 5 75 

-.452;; 

* 'J 1 

• C 5 4 a 2 

>* t 54 

* L4 J2 

• 4 : c . •. 

. 0 5 6 0 5 

-.9c77i 

-. ,55 -5 

,MU3 

• v 5 6 e t 

~.C0 e 1 7 

-.0 5 5 c c 

. y 1" c J 

. < *i 65 4 

“.-•c9w. 

“•>16 54 

. 7 i c x 7 

,C57ut 

” • ’j . y >' 6 

- • o - 7 o c 

.9x513 

.1574^ 

“ • ul ^7, 

— • j 5 7 *i ? 

• V 1 r f '? 

.0 5 77y 

4i 1 54 

-.Ut ? /4 

.9 XC 1 i 

.-5K/4 

"ij’ui! 

-•'■5c J 4 

.9 WV- 

• j 5 *' t 4 

- , w 1 2 3 x 

5 < i 4 

. 4 x 7 . - 

• O 1 C. 4 0 

“.0142) 

- , . 5 C 4 6 

. 07.7 

. 0 > e 5 3 

~ . 1 j 1 5 1 1 


. t i 7 L- c 

• (- 51'. 4 

- • “’1 1 2 2 

- • 6 - c c 4 

.71773 

' t 74 

- . J 1 t c . 4 

*5 6 74 

. 4 1 7: ■, 

. t 5 c f. i 

“ . v 1 7 r t. 

- . 0 ^ t >3 

. 07 ; i 

. r 540? 

- . *• 1 t C 4 

-. a r rz 

. 9 : ? 4 <j 

.65951 

“ * v> 1 9c? 

- • C r V J 1 

. y 1 7:6 

• V ; 9 09 

«2i4 : 

-•-5904 

. vi 725 

• c : 9 i v 

“.Cx 1< 5 

‘ 4 X>- 

.>1713 

.05426 

>? ?' 9 

- , a 9 7 6 

. a 7c 

.05485 

“t 0 7 19 t 

-.-•5*35 

. > 4 c { o 

• 0 5 9 4 4 

“ , 4 2 5 .>6 

9 4 4 

. a 

• 0541 3 

-yV6i6 

- . ■:• 5 9 ;; 3 

. -'If 5 4 

• o ! 963 

X 7 -'3 

- ♦ 0 5 9 u 3 

. a t * (: 

• 4. 0 4 7? 

- . . 2 *: ! x 

-•'.<5 4 n 

• V 1 6 : . 

. 1- 0 4 L 2 

«,i il v 7 . 

- • 5 *. 4 o 2 

• 9 1 - 7 

. t 5 9 4 3 

- . V- 

- • < • : 9 r } 

ft V 1 L f D 

. C t i 3 

- . t 3 t ci 

-.*•6 0 

,91:: /: 

.or t- 14 

-i v? ^4 

— . 6 C* x 4 

.4x52^ 

. 0 C * 2 C 

- , J 3 4 . x 

-• 6 t u 2( 

.SiHl 

. ; t,c i- 

“ . * 2 6' ( 

— * w t ' ' i t 

,914 7 A 

• C t t 5 vi 

- , c 7 ; S 

la. 

. y „ 4 'x 

. ( but 3 

-.03^/5 

- . . f C > 3 

. 9 14 11 

. Cl L 7 h 

— . 0 U ? 4 

- . <v c C< 7 c 

.4x777 

• 5 U * 0 

-. *)4] i 7 

— . . f L 'il. 

.9.24. 

.06104 

-.0433i 

-. , 61 j4 

a a 3* 3 

. ( L 1 ] 9 

- » J4 s 4 2 

-.■>6x14 

• 4 1 2 t 2 

. V' 6 1 3 5 

- # x t <« s 

-.9613 - 

. V 1 ? 1 i 

.66X5 i 

- , U*t c i l 

-.o6 j 

. >1173 

.00167 

“.04475 

. - . 1 1 *J 7 

aa i a 

. 9< U * 

- • 0 ! 1 *. x 

- • •> t 15 5 

,>1.74 

• G 6 a y 3 

- » j 5 3 1 4 

— , 0 6 c \j 3 

■ » , u ■ ■ 

, c 2 2 * 
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